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In this review, we rederive the controversial influence functional approach of Golubev 
, and Zaikin (GZ) for interacting electrons in disordered metals in a way that allows us to 

^ ■ show its equivalence, before disorder averaging, to diagrammatic Keldysh perturbation 

ff^ ' theory. By representing a certain Pauli factor (5 — 2fP) occuring in GZ's efi^ective action 

' in the frequency domain (instead of the time domain, as GZ do), we also achieve a more 

l/^ , accurate treatment of recoil effects. With this change, GZ's approach reproduces, in a 

C D 1 remarkably simple way, the standard, generally accepted result for the decoherence rate. 

- The main text and appendices A.l to A. 3 of the present review are comparatively 
1/^ ' brief, and have been published previously; for convenience, they are included here again 

, (with minor revisions). The bulk of the review is contained in several additional, lengthy 

appendices containing the relevant technical details. 
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1. Introduction 

A few years ago, Golubev and Zaikin (GZ) developed an influen ce fianct ion al ap- 
proach for describing interacting fermions in a disordered 

conductojmm. 

Their 

key idea was as follows: to understand how the diffusive behavior of a given electron 
is affected by its interactions with other electrons in the system, which constitute 
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^ . its effective environment, the latter should be integrated out, leading to an influ- 

ence functional, denoted by e~Ti-'^'^^"-^^'\ in the path integral J'D'R describing its 
dynamics. To derive the effective action {iSn + Sj), GZ devised a strategy which, 
when implemented with sufficient care, properly incorporates the Pauli principle - 
this is essential, since both the particle and its environment originate from the same 
system of indistinghuishable fermions, a feature which makes the present problem 
conceptually interesting and sets it apart from all other applications of influence 
functionals that we are aware of. 

GZ used their new approach to calculate the electron decoherence rate 7,^ (T) in 
disordered conductors, as extracted from the magnetoconductance in the weak lo- 



I 



I 



April 10, 2008 10:57 WSPC/INSTRUCTION FILE ijmpb-review-vondelft 



2 Jan von Delft 

calization regime, and found it to be finite at zero temperaturJ^EEISEIZ!^ J^^{T 
0) — 7^'*^^, in apparent agreement with some experiment J^l. However, this re- 
suh contradicts the standard view, based on the work of Altshuler, Aronov and 
Khmelnitskii {AAK'P^, that 7^'^^(T ^ 0) = 0, a nd hence elicited a consid- 
erable controversj0. GZ's work was widely questioned,[IIIIlIIlIlIIIIIS!, with the 
most detai led and vigorous critique coming from Alei ner, Al tshuler and Gershen- 
zon (A Acjl^a nd Aleiner, Altshuler and Vavilov (AAvjlSHH^ but GZ rejected each 
critiquJMm 

with equal vigor. It is important to emphasize that the debate here 
was about a well-defined theoretical model, and not about experiments which do or 
do not support GZ's claim. 

The fact that GZ's final results for 7^^(r) have been questioned, however, does 
not imply that their influence functional approach, as such, is fundamentally flawed. 
To the contrary, we show in this review that it is sound in principle, and that the 
standard result J^^^{T) can be reproduced using GZ's method, provided that it is 
applied with slightly more care to correctly account for recoil effects (i.e. the fact 
that the energy of an electron changes when it absorbs or emits a photon). We 
believe that this finding conclusively resolves the controversy in favor of AAK and 
company; hopefully, it will also serve to revive appreciation for the merits of GZ's 
influence functional approach. 

The premise for understanding how 7^"^^ can be reproduced with GZ's methods 
was that we had carried out a painfully detailed analysis and rederivation GZ's ap- 
proach, as set forth by them in two lengthy papers from 1999 and 2000, henceforth 
referred to as GZ99^ and GZoP. Our aim was to establish to what extent their 
method is related to the standard Keldysh diagrammatic approach. As it turned 
out, the two methods are essentially equivalent, and GZ obtained unconventional 
results only because a certain "Pauli factor" {6 — 2fP) occuring in Sr was not 
treated sufficiently carefully, where fp is the single-particle density matrix. That 
their treatment of this Pauli factor was dubious had of course been understood and 
emphasized before: first and foremost it was correctly pointed out by AAgP^ that 
GZ's treatment of the Pauli factor caused their expression for 7^^ to aquire an ar- 
tificial ultraviolet divergence, which then produces the term 7^'^^, whereas no such 
divergence is present in diagrammatic calculations. GZ's treatment o f (5 — 2(P) was 
also criticized, in various related contexts, by several other 

authorjll|l^|15|16118]. 

However, none of these works (including our owrP^, which, in retrospect, missed 
the main point, namely recoil) had attempted to diagnose the nature of the Pauli 
factor problem with sufficient precision to allow a successful remedy to be devised 
within the influence functional framework. 

This will be done in the present review. Working in the time domain, GZ rep- 
resent {S - 2p°{t)) as 1 - 2nQ[ho{t)/2T], where uq is the Fermi function and ho{t) 
the free part of the electron energy. GZ assumed that ho{t) does not change during 
the diffusive motion, because scattering off impurities is elastic. Our diagnosis is 
that this assumption unintentionally neglects recoil effects (as first pointed out by 
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Eriksen and Hedeear cP) , because the energy of an electron actually does change 
at each interaction vertex, i.e. each time it emits or absorbs a photon. The remedy 
(not found by Eriksen and Hedegard) is to transform from the time to the frequency 
domain, in which {S — 2p°) is represented by 1 — 2nQ[h{e — uj)] — ta.nh[h{e — oj) /2T], 
where hCu is the energy change experienced by an electron with energy he at an 
interaction vertex. Remarkably, this simple change of representation from the time 
to the frequency domain is sufficient to recover 7^"^^. Moreover, the ensuing calcu- 
lation is free of ultraviolet or infrared divergencies, and no cut-offs of any kind have 
to be introduced by hand. 

The main text of the present review has two central aims: firstly, to concisely 
explain the nature of the Pauli factor problem and its remedy; and secondly, to 
present a transparent calculation of j^, using only a few lines of simple algebra. 
(Actually, we shall only present a "rough" version of the calculation here, which re- 
produces the qualitative behavior of 7^^^ (T) ; an improved version, which achieves 
quantitative agreement with AAK's result for the magnetoconductance [with an er- 
ror of at most 4% for quasi- 1-D wires], has been published in a separate analysis by 
Marquardt, von Delft, Smith and Ambegaokar^. The latter consists of two parts, 
referred to as MDSA-I and DMSA-II below, which use alternative routes to arrive 
at conclusions that fully confirm the analysis of this review.) 

We have made an effort to keep the main text reasonably short and to the 
point; once one accepts its starting point [Eqs. ([T]) to Eq. Q], the rest of the 
discussion can easily be followed step by step. Thus, as far as possible, the main text 
avoids technical details of interest only to the experts. These have been included 
in a set of five lengthy and very detailed appendices, B to F, in the belief that 
when dealing with a controversy, all relevant details should be publicly accessible 
to those interested in "the fine print" . For the benefit of those readers (presumably 
the majority) with no time or inclination to read lengthy appendices, a concise 
appendix A summarizes (without derivations) the main steps and approximations 
involved in obtaining the influence functional. 

The main text and appendices A.l to A. 3 have already been published 
previousljD], but for convenience are included here again (with minor revisions, 
and an extra sketch in Fig. [T]), filling the first 23 pages. The content of the re- 
maining appendices is as follows: In App. IA.4l we address GZ's claim that a strictly 
nonperturbative ap proach is needed for obtaining 7^, and explain why we disagree 
(as do many other^i^JiSLkyj). In App. B, we rederive the influence functional and 
effective action of GZ, following their general strategy in spirit, but introducing 
some improvements. The most important differences are: (i) instead of using the 
coordinate- momentum path integral / T>R J VP of GZ, we use a "coordinates-only" 
version J "D'R, since this enables the Pauli factor to be treated more accurately; and 
(ii), we are careful to perform thermal weigthing at an initial time to —00 (which 
GZ do not do), which is essential for obtaining properly energy-averaged expres- 
sions and for reproducing perturbative results: the standard diagrammatic Keldysh 
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perturbation expansion for the Cooperon in powers of the interaction propagator 
is generated if, before disorder averaging, the influence functional is expanded in 
powers of {iSii + Si)/h. In App. C we review how a general path integral expression 
derived for the conductivity in App. B can be rewritten in terms of the familiar 
Cooperon propagator, and thereby related to the standard relations familiar from 
diagrammatic perturbation theory. In particular, we review the Fourier transforms 
required to obtain a path integral P^g{T) properly depending on both the energy 
variable he relevant for thermal weighting and the propagation time r needed to 
traverse the closed paths governing weak localization. Appendix D gives an explicit 
time-slicing definition of the "coordinates-only" path integral / VR used in App. B. 
Finally, for reference purposes, we collect in Apps. E and F some standard mate- 
rial on the diagrammatic technique (although this is bread-and-butter knowledge 
for experts in diagrammatic methods and available elsewere, it is useful to have it 
summarized here in a notation consistent with the rest of our analysis). App. E 
summarizes the standard Keldysh approach in a way that emphasizes the anal- 
ogy to our influence functional approach, and App. F collects some standard and 
well-known results used for diagrammatic disorder averaging. Disorder averaging 
is discussed last for a good reason: one of the appealing features of the influence 
functional approach is that most of the analysis can be performed before disorder 
averaging, which, if at all, only has to be performed at the very end. 

2. Main Results of Influence Functional Approach 

We begin by summarizing the main result of GZ's influence functional approach. 
Our notations and also the content of some of our formulas are not identical to 
those of GZ, and in fact differ from their's in important respects. Nevertheless, we 
shall refer to them as "GZ's results", since we have (re)derived them (see App. B 
for details) in the spirit of GZ's approach. 

The Kubo formula represents the DC conductivity ctdc in terms of a retarded 
current-current correlator ([ji(l), j(2)]). This correlator can (within various approx- 
imations discussed in App. B.5.6, B.5.7, B.6.3 and IA.3|1 be expressed as follows in 
terms of a path integral P^g representing the propagation of a pair of electrons with 
average energy He, thermally averaged over energies: 

CTDC^^ Jdx2 jii'- j 22' J {de)[-n' {he)] J dTp^l'^gir), (la) 

The propagator P2v'lsi'^)^ defined for a given impurity configuration, is written 
in terms of a forward and backward path integral ^ ^ V'R between the specified 

initial and final coordinates and times. It gives the amplitude for a pair of electron 
trajectories, with average energy He, to propagate from r2' at time — to ri at 
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ir or from ry at time to r2 at — jT, respectively. [The sense in which both r 
and e can be specified at the same time is discussed in App. IA.3[ and in more detail 
in App. IC.4i Eqs. (|C.21|) to (|C.24[) ]. We shall call these the forward and backward 
paths, respectively, using an index a = F,B to distinghuish them. Sq = 
are the corresponding free actions, which determine which paths will dominate the 
path integral. The weak localization correction to the conductivity, ctqq , arises 
from the "Cooperon" contributions to ctdCj illustrated in Fig. [TJb), for which the 
coordinates ri, r^, and r'2 all lie close together, and which feature self-returning 
random walks through the disordered potential landscape for pairs of paths i?^/^, 
with path B being the time-reversed version of path F, i.e. RF {t^i) = R^{—t3) for 
^3 e (— ^T, ^r). The effect of the other electrons on this propagation is encoded 
in the influence functional ^-ii^SR+Si)/?! occuring in Eq. (jlbp . The effective action 
iSii + Si turns out to have the form [for a more explicit version, see Eq. (|A.7|) in 
App. A; or, for an equivalent but more compact representation, see Eqs. (|B.94p and 
(fBirll of Sec. EES]: 

rt'in ( ra' ^ 

(2) 

Here Sa stands for Sp/B — if; and the shorthand ^3^4'^ — ^[ts^ — t4^, , R°-{t3^) — 
(^4„/ )] describes, in the coordinate-time representation, an interaction propaga- 
tor linking two vertices on contours a and a' . It will be convenient below to Fourier 
transform to the momentum-freqency representation, where the propagators C 
and £° can be represented as follows [{duj){dq) = {duj dq) / {2tt)^]: 

^tK,= iduj)idq)e^^ ' J -"J£^{uj), (3a) 

^3.4„,-|[^^(^_2p0)]^^3^ if a'^B, ^''^ 
= J (dw)(dq)e*^"'(''"[''°^*^''^"''°'^*^"'^]~"''*^""*v))£^_'(^) . (3c) 

[Note the sign Sa' in the Fourier exponential in Eq. ([3c|) : it reflects the opposite 
order of indices in Eq. (j3bp . namely 34 for F vs. 43 for B.] Here is the Keldysh 
interaction propagator, while C^^^ , to be used when time ^4^, lies on the forward or 
backward contours, respectively, represent "effective" retarded or advanced prop- 
agators, modifled by a "Pauli factor" {S — 2p°) (involving a Dirac-delta Sij and 
single-particle density matrix p| in coordinate space), the precise meaning of which 
will be discussed below. £^'^''^(<D) denote the Fourier transforms of the standard 
Keldysh, retarded, or advanced interaction propatators. For the screened Coulomb 
interaction in the unitary limit, they are given by 

-R -A E^-iiD [^"(w)!"^ 
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Fig. 1. (a) Structure of vertices on the forward or backward contours of Keldysh perturbation 
theory. F: the combinations CI^^^^C^;^^ and G^^^^C^^^ occur if vertex 4 lies on the upper 
forward contour. B: the combinations ^G¥ and C¥ occur if vertex 4 lies on the 

lower contour. Arrows point from the second to first indices of propagators, (b) Sketch of a pair of 
time-reversed paths connecting the points at which the current operators jn'-j 22' s-^t [cf. Eq. Illal l] , 
decorated by several (wavy) interaction propagators C^^f^^ ^ (lo). In the Keldysh formalism, the 
electron lines represent the electron propagators G^^^{uj) or G^ (uj) = ta.nh{hLj/2T)[G^ — G^]{uj). 
The effective action defined in Eqs. ^ to Il4al l in effect neglects the frequency transfers uji in 
the arguments of all retarded and advanced electron Green's functions [G^^^{e — uii — . . .) 
G«/'*(e)], but, for every occurence of the combination C^/^{uJi)G^ {e — uJi), retains it in the 
factor tanh[/i{e — LOi)/fi\ of the accompanying G^ function. The latter prescription ensures that 
a crucial feature of the Keldysh approach is retained in the influence functional formalism, too, 
namely that all integrals J du)i over frequency transfer variables are limited to the range | fu^i \ < T 
[which is why the neglect of uJi in G^^^{£ — uJi — . . .) is justified]. In contrast, GZ also neglect 
the —UJi in tanh[?i(e — iL}i)/fi\ [see Sec.|4], which amounts to neglecting recoil. As as a result, their 
J du)i integrals are no longer limited to \hLOi\ < T, i.e. artificial ultraviolet divergencies occur, 
which produce GZ's temperature-independent contribution 7^^''"'^ to the decoherence rate [see 
Eq. {TTJ]. Thus, 7^^'*^^ is an artefact of GZ's neglect of recoil, as is their claimed "decoherence at 
zero temperature" . 



Cq{uj) = 2icoth{hOj/2T)lm[Cg{u))] , (4b) 

El = Oct , Eg = Oct + IH , (4d) 

where, for later reference, we have also listed the Fourier transforms of the bare 
diffuson V and Cooperon C (where 7// is the dephasing rate of the latter in the 
presence of a magnetic field, D the diffusion constant and v the density of states 
per spin). Finally, Cq {Co) in Eq. (j3c[) is defined as 

^(w) = tanh[ft(e - ub)/2T] cf^{u!) , (4e) 

where he is the same energy as that occuring in the thermal weighting factor 
[~n'{fi£)] in Eq. pa|) . 

Via the influence functional, the effective action ^ concisely incorporates the 
effects of interactions into the path integral approach. Si describes the classical 
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part of the effective environment, and if one would replace the factor coth(huj/2T) 
in {uj) by 2T/hw (as is possible for high temperatures) it corresponds to the 
contribution calculated by AAK?^. With Sr, GZ succeeded to additionally also 
include the quantum part of the environment, and in particular, via the Pauli factor 
{5 — 2p"), to properly account for the Pauli principle. 

Casual readers are asked to simply accept the above equations as starting point 
for the remainder of this review, and perhaps glance through App. A to get an idea 
of the main steps and approximations involved in deriving them. Those interested in 
a detailed derivation are referred to App. B (where Sr/j are obtained in Sec. B.5.8). 
It is also shown there [Sec. B.6] that the standard results of diagrammatic Keldysh 
perturbation theory can readily be reproduced from the above formalism by ex- 
panding the influence functional g-('^SR+Si)/h powers of {iSR + Si)/h. For present 
purposes, simply note that such an equivalence is entirely plausible in light of the 
fact that our effective action ^ is linear in the effective interaction propagators C, 
a structure typical for generating functionals for Feynman diagrams. 



3. Origin of the Pauli Factor 

The occurence of the Pauh factor {S - 2p°) in Sr was first found by GZ in precisely 
the form displayed in the position-time representation of the effective action used 
in Eq. However, their subsequent treatment of this factor differs from ours, in 
a way that will be described below. In particular, they did not represent this factor 
in the frequency representation, as in our Eq. (j4ep . and this is the most important 
difference between our analysis and theirs. 

The origin of the Pauli factor in the form given by our Eq. (j4e|) can easily be 
understood if one is familiar with the structure of Keldysh perturbation theory. 
[For a detailed discussion, see Sec. B.6. 2.] First recall two exact relations for the 
noninteraction Keldysh electron propagator: in the coordinate-time representation, 
it contains a Pauli factor, 

G^^ jdxu (G«-G^).,(^-2p%=|dxfe {5-2~p\,{G''-G^)u, 

(5a) 

which turns into a tanh in the coordinate frequency representation: 

G^{Cj) = tanh(to/2r) [cf (w) - G^(w)j . (5b) 

Now, in the Keldysh approach, retarded or advanced interaction propagators always 
occur [see Fig. [TJa)] together with Keldysh electron propagators, in the combina- 
tions G'j^p4pC!^4^ or L^^^G^^j^, where the indices denote coordinates and times. 
[Likewise, the Keldysh interaction propagators always come in the combinations 
^fpiF^^iF ^iB?'^^BjB-\ -'^-'^ ^"^^ momentum-frequency representation, the combi- 

^ RJ _ _ 

nations involving G^ therefore turn into (a)) \G^- — G"^] ^_-{e — u)) tanh[S(£ — 
Co)/2T]. Thus, in the frequency representation the Pauli factor is represented as 
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tanh[h{e ~ uj)/2T]. Here the variable he represents the energy of the electron line 
on the upper (or lower) Keldysh contour before it enters (or after it leaves) an in- 
teraction vertex at which its energy decreases (or increases) by to [see Fig. [Ija)]. 
The subtraction of Co in the argument of tanh thus reflects the physics of recoil: 
emitting or absorbing a photon causes the electron energy to change by fid), and it 
is this changed energy h{e — uj) that enters the Fermi functions for the relevant final 
or initial states. 

Of course, in Keldysh perturbation theory, he will have different values from one 
vertex to the next, reflecting the history of energy changes of an electron line as it 
proceeds through a Feynman diagram [as illustrated in Fig. [Tfb)]. It is possible to 
neglect this complication in the influence functional approach, if one so chooses, by 
always using one and the same energy in Eq. (|4ep . which then should be chosen to 
be the same as that occuring in the thermal weighting factor [—n'{he)], i.e. he = he. 
This approximation, which we shall henceforth adopt, is expected to work well if the 
relevant physics is dominated by low frequencies, at which energy transfers between 
the two contours are sufficiently small [h{e — e) ^ T, so that the electron "sees" 
essentially the same Fermi function throughout its motion. [For a detailed discussion 
of this point, see App. IB. 6. 21 ] 

Though the origin and neccessity of the Pauli factor is eminently clear when seen 
in conjunction with Keldysh perturbation theory, it is a rather nontrivial matter 
to derive it cleanly in the functional integral approach [indeed, this is the main 
reason for the length of our appendices!]. The fact that GZ got it completely right 
in the position-time representation of Eq. ^ is, in our opinion, a significant and 
important achievement. It is regrettable that they did not proceed to consider the 
frequency representation (|4ep . too, which in our opinion is more useful. 

4. Calculating a la GZ 

To calculate the decoherence rate jip — I/tj^, one has to find the long-time decay 
of the Cooperon contribution to the propagator P^gir) of Eq. ([T]). To do this, GZ 
proceeded as follows: using a saddle-point approximation for the path integral for the 
Cooperon, they replaced the sum over all pairs of self-returning paths R^^^ (t^j^^^ ) 

by just the contribution {e^i'-'^^"^^''''^^^)rw of the classical "random walk" paths 
Riw{t) picked out by the classical actions Sq, namely R^{t3j,) = Ri-vj{t3p) and 
R^{t3j^) = i2rw(— ^815)7 for which the paths on the forward and backward Keldysh 
contours are time-reversed partners. The subscript "rw" indicates that each such 
classical path is a self-returning random walk through the given disorder potential 
landscape, and ( )rw me ans averaging over all such paths. Next, in the spirit 
of Chakravarty and Schmid^, they replace the average of the exponent over all 
time-reversed pairs of self-returning random walks, by the exponent of the average, 
e--^(^), where F{t) = ^{iSn + S*/)™ (cf. Eq. (67) of GZOG^). This amounts to 
expanding the exponent to first order, then averaging, and then reexponentiating. 
The function F{t) thus defined increases with time, starting from F{0) — 0, and 
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the decoherence time can be defined as the time at which it becomes of order 
one, i.e. F{t^) w 1. 

To evahiate {iSu + S'/)i.w, GZ Fourier transform the functions £3^4' = 
£[^34, R'^{t3) — R°-' (ti)] occuring in Sji/j, and average the Fourier exponents usinj^ 
the distribution function for diffusive motion, which gives probabihty that a random 
walk that passes point R^witi) at time will pass point R,-^{t3) at time ^3, i.e. 
that it covers a distance R = Ri-^lt^) — Ri-^lt^) in time |t34|: 

v^l*34iy 

=e-g'^|t34| c0(|t34|) = e-^4l*^*l . (6) 

(Here t^^ = — t^.) The arrow in the second line makes explicit that if we also ac- 
count for the fact that such time-reversed pairs of paths are dephased by a magnetic 
field, by adding a factor e""*""'*^*', the result is simply equal to the bare Cooperon 
in the momentum-time representation. 

Actually, the above way of averaging is somewhat inaccurate, as was pointed 
out to us by Florian Marquardt: it neglects the fact that the diffusive trajectories 
between ^3 and ^4 are part of a larger, self-returning trajectory, starting and ending 
at ri ~ r2 at times ^\t. It is actually not difficult to include this fact, see MDSA- 



i^QI, and this turns out to quantitatively improve the numerical prefactor for r^p (e.g. 
in Eq. p8|) below). However, for the sake of simplicity, we shall here be content with 
using Eq. ([6]), as GZ did. 

Finally, GZ also assumed that the Pauli factor [5 — 2fP) in Su remains un- 
changed throughout the diffusive motion: they use a coordinate-momentum path 
integral jT^RfVP [instead of our coordinates-only version JVR], in which 
{S — 2ff) is replaced by [1 — 2no{ho)] = tanh(/io/2T), and the free-electron en- 
ergy ho[R{ta), Pita)] is argued to be unchanged throughout the diffusive motion, 
since impurity scattering is elastic [cf. p. 9205 of "n depends only on the 

energy and not on time because the energy is conserved along the classical path"]. 
Indeed, this is true between the two interaction events at times ^3 and ^4, so that 
the averaging of Eq. ^ is permissible. However, as emphasized above, the full tra- 
jectory stretches from — to ^4 to ^3 to ^r, and the electron energy does change, 
by ±ftd), at the interaction vertices at t^ and t^. Thus, GZ's assumption of a time- 
independent Pauli factor neglects recoil effects. As argued in the previous section, 
these can straightforwardly taken into account using Eq. (j4e]), which we shall use 
below. In contrast, GZ's assumption of time- independent n amounts dropping the 
—Hu) in our tanh[/i(e — lu)/2T] function. 

If one uses GZ's assumptions to average Eq. ([2]), but uses the proper tanh[ft(e — 
w)/2T] function, one readily arrives at 




= 2Rc 



iz /(rfc^)(dg)rj^''H [r"-/H(r) 



-2« 



(7) 
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where /^^^'^ — f^"^ are the first and second terms of the double time integral 



corresponding to self-energy (a — a' — F) and vertex (a ^ a' = F) contributions, 
and the 2Re[ ] in Eq. ([7]) comes from adding the contributions of a' = F and B. 
Performing the integrals in Eq. ([8]), we find 



Cli-L0)r + \cli- 



- 1 



(9a) 
(9b) 



Of all terms in Eqs. ([9]), the first term of f^°^^, which is linear in t, clearly grows most 
rapidly, and hence dominates the leading long-time behavior. Denoting the associ- 
ated contribution to Eq. © by ^(i^fl/S'/)[°^'^'"s,soit ^ ^^«/A>5<=if^ corresponding 

rates 7,^^^''^'''^ obtained from Eqs. ^ and 



{duj){dq) tanh 



h{e — uj) 



2T 



are: 



2Re 



ILU) 



= — / {duj){dq) coth 



2T 



2Re 



Let us compare these results to those of GZ, henceforth using 



(10a) 

(10b) 
0. Firstly, both 

our 7;^'"''" and 7^'*'*''^ are nonzero. In contrast, in their analysis GZ concluded that 
{SrItw — 0. The reason for the latter result is, evidently, their neglect of recoil 
effects: indeed, if we drop the —hoj from the tanh-factor of Eq. ()10ap . we would find 
7^ = and thereby recover GZ's result, since the real part of the factor in square 
brackets is odd in uj. 

Secondly and as expected, we note that Eq. (jlObp for 7^^''''''^ agrees with that of 
GZ, as given by their equation (71) of GZ993 for 1/t^, i.e. 7^''^°'^ = 7^^. [To see 
the equivalence explicitly, use Eq. (IA.9P .] Noting that the Jdw-integral in Eq. (jlObp 
evidently diverges for large Cu, GZ cut off this divergence at I/t^i (arguing that 
the diffusive approximation only holds for time-scales longer than Td, the elastic 
scattering time). For example, for quasi- 1-dimensional wires, for which J{dq) = 



'J dq/{2-K) can be used (a^ being the cross section, so that ai = 0^0-^™'^° 



IS 



the conductivity per unit length, with (t^q 
GZ9P): 



Drudc 



2e^vD), they obtain (cf. (76) of 



1 

rGZ 



'2D 



hai 



{duj) 



coth 



2T 




1 



(11) 



[The use of a self-consistently-determincd lower frequency cut-off is explained in 
Sec. [B]. Thus, they obtained a temperature-independent contribution 7^'*^^ from 
the +1 term, which is the result that ingited the controversy. 
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However, we thirdly observe that, due to the special form of the retarded inter- 
action propagator in the unitary limit, the real parts of the last factors in square 
brackets of Eqs. (jlOap and (jlOb|l are actually equal (for •jh — 0). Thus, the ultravio- 
let divergence of 7^'''°'^ is cancelled by a similar divergence of 7^'*^°'^ Consequently, 
the total decoherence rate coming from self-energy terms, 7^®'^ = 7;^''^'''^ + 7,^'^°'^, is 
free of ultraviolet divergencies. Thus we conclude that the contribution 7"'*^^ found 
by GZ is an artefact of their neglect of recoil, as is their claimed "decoherence at 
zero temperature" . 



5. Dyson Equation and Cooperon Self Energy 

The above results for 7^'''°'^ -I- 7^'''°'^ turn out to agree completely with those of a 
standard calculation of the Cooperon self energy S using diagrammatic impurity 
averaging [details of which are summarized in Appendix F] . We shall now summarize 
how this comes about. 

Calculating S is an elementary excercise within diagrammatic perturbation the- 
ory, first performed by Fukuyama and Abrahami^. However, to facilitate compar- 
ison with the influence functional results derived above, we proceed differently: We 
have derived [Sec. B.6.1] a general expressiorPSl, before impurity averaging, for the 
Cooperon self-energy of the form E = 12aa' ^aa' + ^aa' j which keeps track of 
which terms originate from iSji or Sj, and which contours a, a' = F/ B the vertices 
sit on. This expression agrees, as expected, with that of Keldysh perturbation the- 
ory, before disorder averaging; it is given by Eq. (|A.10|) and illustrated by Fig. lAll 
in App. A. We then disorder average using standard diagrammatic techniques. For 
reference purposes, some details of this straightforward excercise are collected in 
Appendix F.2. 

For present purposes, we shall consider only the "self-energy contributions" 
{a = a') to the Cooperon self energy, and neglect the "vertex contributions" (a a'), 

R/I 

since above we likewise extracted 71^ from the self-energy contributions to the 
effective action, {Sh/jY^^'^™^'^'^^^. After impurity averaging, the Cooperon then sat- 

self 

isfies a Dyson equation of standard form, Cg (uj) = Cq(w) +Cg{u;) 
with standard solution: 

1 



(w) 



(12) 



where S ^}_^^ S„„ , with ^^ pp (w) 



"— i?y/,scif 
q,BB 



and 



■"q.FF 



{duj)(dq) coth 



2r 



Im[Zj(w)] C^_q(t- 



(13a) 



^fl.self 1 



{duj){dq)<. tanh 



2T 



iCq (w) 



(13b) 
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In Eq. (jl3bp . the terms proportional to (l?")^ [(C*^) ^ + (2?") J stem from the 
so-called Hikami contributions, for which an electron line changes from G^/^ to 
qA/r |.q qR/a ^jjg interaction vertices. As correctly emphasized by AAcff^ 
and AA\flS!, such terms are missed by GZ's approach of averaging only over time- 
reversed pairs of paths, since they stem from paths that are not time-reversed pairs. 

Now, the standard way to define a decoherence rate for a Cooperon of the form 
(fT2|) is as the "mass" term that survives in the denominator when uj = Eg = Q, 



I.e. 7; 



self 
V 



'(0) 



-2Re 



-/-|-fi,self , 



In this limit the contribution of the 



^0,FF (0) 

Hikami terms vanishes identically, as is easily seen by using the last of Eqs. ()4a|) . and 
noting that Re[i(P )~^] = R.e[i] = 0. (The realization of this fact came 

to us as a surprise, since A AG and AAV had argued that GZ's main mistake was 
their neglect of Hikami termJI^II^, thereby implying that the contribution of these 
terms is not zero, but essential.) The remaining (non-Hikami) terms of Eq. (jl3bp 
agree with the result for E of AAV^^ and reproduce Eqs. (fTO)) given above, in other 
words: 



1. 



—self 



1 



= hEo (0)]= {^SR + SI 



leading, sell 
rw 



(14) 



Thus, the Cooperon mass term — Sq (0) agrees identically with the coefficient of 
T in the leading terms of the averaged effective action of the influence functional. 
This is no coincidence: it simply reflects the fact that averaging in the exponent 
amounts to reexponentiating the average of the first order term of an expansion 
of the exponential, while in calculating the self energy one of course also averages 
the first order term of the Dyson equation. It is noteworthy, though, that for the 
problem at hand, where the unitary limit of the interaction propagator is considered, 
it suffices to perform this average exclusively over pairs of time-reversed paths — 
more complicated paths ar e ev idently not needed, in contrast to the expectations 
voiced by AAcfl^and AA\tIHl. 

The latter expectations do apply, however, if one consideres forms of the in- 
teraction propagator (lj) more general than the unitary limit of (|4ap (i.e. not 

proportional to [P°(ci;)]~^). Then, the Hikami contribution to 7^°'^ — I]q°'^(0) in- 
deed does not vanish; instead, by noting that for uj = q = jh ~ the second line 
of Eq. (jl3bp can always be written as 2Re[P°(a})] , we obtain 



„self 



- J {du)){dq) |coth 
xIm[£^((D)^ 



q 



2T 
24 



tanh 



2T 



(15) 



which is the form given by AA\iISl. 
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6. Vertex Contributions 

Eq. (|10bp for has the deficiency that its frequency integral is infrared diver- 

gent (for ill ^ 0) for the quasi-1 and 2-dimensional cases, as becomes explicit once 
its q-integral has been performed [as in Eq. (fTT|) ]. This problem is often dealt with 
by arguing that small-frequency environmental fluctuations that are slower than the 
typical time scale of the diffusive trajectories are, from the point of view of the dif- 
fusing electron, indistuingishable from a static field and hence cannot contribute to 
decoherence. Thus, a low-frequency cutoff j^p is inserted by hand into Eqs. pil)) [i.e. 
did — > dio] , and j^p determined selfconsistently. This procedure was motivated 

in quite some detail by AAG"'^''', and also adopted by GZ in [see Eq. (|lip 

above]. However, as emphasized by GZ in a subsequent paper, GZOC-, it has the 
serious drawback that it does not necessarily reproduce the correct functional form 
for the Cooperon in the time domain; e.g., in d = 1 dimensions, the Cooperon 
is knowrff^ to decay as e-'^i^/^^r^ i.e. with a nontrivial power in the exponent, 
whereas a "Cooperon mass" would simply give e"'^/'^^. 

A cheap fix for this problem would be to take the above idea of a self-consistent 
infrared cutoff one step further, arguing that the Cooperon will decay as e~'^'^f ^'^\ 
where 7^°'^(t) is a time- dependent decoherence rate, whose time-dependence enters 
via a time-dependent infrared cutoff. Concretely, using Eqs. ^3]) and pO]) . one would 
write 



l/r 







■jcoth 






2T 



i s tanh 



s=± 



h{e — suj) 
2T 



{dq) l_ 



(16) 



It is straightforward to check [using steps analogous to those used below to obtain 
Eq. (fT8|) ] that in J = 1 dimensions, the leading long-time dependence is 7^°'^(t) cx 
r^/^, so that this cheap fix does indeed produce the desired e~°-^'^^'^-f^ ' behavior. 

The merits of this admittedly rather ad hoc cheap fix can be checked by doing a 
better calculation: It is well-known that the proper way to cure the infrared prob- 
lems is to include "vertex contributions" , having interactions vertices on opposite 
contours. In fact, the original calculation of AAKP^in effect did just that. Likewise, 
although GZ neglected vertex contributions in they subsequently included 

them in exploiting the fact that in the influence functional approach this 

is as straightforward as calculating the self-energy terms: one simply has to in- 
clude the contributions to {iSR/Si)rw of the vertex function — J™''* in Eq. (O, 
too. The leading contribution comes from the first term in Eq. (|9bp . to be called 
(i5fl,/S'7);.°^^'^'"e'™''*, which gives a contribution identical to {iSn/ SiY^^'^'''^''"'^^ , but 
multiplied by an extra factor of —^^^^7^ under the integral. Thus, if we collect all 
contributions to Eq. ([7]) that have been termed "leading" , our final result for the 
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averaged effective action is j-{iSii + ■5'/)™'*'"^ = F^{t), with 



Fdi^) = / {duo) uo |coth 
{dq) 1 



hio 

2T 



tanh 



h{e — uj) 



2T 



^ sin(ciJT) 



hiy {Dq 



2^2 



(17) 



This is our main result: an expression for the decoherence function F^{t) that is 
both ultraviolet and infrared convergent (as will be checked below), due to the 
(coth + tanh) and (1 — sin)-combinations, respectively. Comparing this to Eqs. (jl6p . 
we note that F^{t) has precisely the same form as T7^°'^(r), except that the infrared 
cutoff now occurs in the / {dui) integrals through the (1 — sin) combination. Thus, 
the result of including vertex contributions fully confirms the validity of using the 
cheap fix replacement J^{dLd) — > J-^^^{did), the only difference being that the cutoff 
function is smooth instead of sharp (which will somewhat change the numerical 
prefactor of r^). 

It turns out to be possible to also obtain Eq. p7|) [and in addition all the 
"subleading" terms of Eq. ((T])] by purely diagrammatic means: to this end, one has 
to set up and solve a Bethe-Salpeter equation. This is a Dyson-type equation, but 
with interaction lines transferring energies between the upper and lower contours, so 
that a more general Cooperon Cq(rii, ^2)1 with three frequency variables, is needed. 
Such an analysis will be published in DMSA-li2QI. 

To wrap up our rederivation of standard results, let us perform the integrals 
in Eq. (|17p for F^{t) for the quasi- 1-dimensional case d = 1. The /(dq)-integral 
yields uo^'^l'^ \J D j^j {a\hl To do the frequency integral, we note that since 
the (coth + tanh)-combination constrains the relevant frequencies to be jtoj ;< T, 
the integral is dominated by the small-frequency limit of the integrand, in which 
coth(/MD/2T) ~ 2T/to, whereas tanh, making a subleading contribution, can be 
neglected. The frequency integral then readily yields 

F^{t) = —=—4=-^—=[-] , (18) 



so that we correctly obtain the known e^°-i'^/'^v)'^' decay for the Cooperon. Here 
g^{L) = {h/e'^)a^L'^~'^ represents the dimensionless conductance, which is 3> 1 
for good conductors. The second equality in Eq. (|18p defines t^, where we have 
exploited the fact that the dependence of Fi on t is a simple r'^/^ power law, which 
we made dimensionless by introducing the decoherence time . [Following AAgP^, 
we purposefully arranged numerical prefactors such that none occur in the final 
Eq. (fT9|) for below.] Setting t = in Eq. (fT8|) we obtain the self-consistency 
relation and solution (cf. Eq. (2.38a) of AAqIS): 



J_ _ T/h _ / h^ai ^ 



(19) 
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The second relation is the celebrated result of AAK, which diverges for T ^ 0. This 
completes our recalculation of 7^^^ using GZ's influence functional approach. 
Eq. (fT8|) can be used to calculate the magnetoconductance for d = 1 via 

Diudc roc 

<^'Sh{H)--"-^ dr(7^o(r)e-^^M. (20) 

TTVh Jq 

(Here, of course, we have to use 7// 7^ in C^^q{t). Comparing the result to AAK's 
result for the magnetoconductance (featuring an Ai' function for d — 1), one finds 
qualitatively correct behavior, but deviations of up to 20% for small magnetic fields 
H. The reason is that our calculation was not sufficiently accurate to obtain the cor- 
rect numerical prefactor in Eq. (llSp. GZ did not attempt to calculate it accurately, 
either]. It turns out (see MDSA-I^^ that if the averaging over random walks of 
Eq. ^ is done more accurately, following Marquardt's suggestion of ensuring that 
the random walks are self-returning, the prefactor changes in such a way that the 
magnetoconductance agrees with that of AAK to within an error of at most 4%. An- 
other improvement that occurs for this more accurate calculation is that the results 
are well-behaved also for finite 'Jh, which is not the case for our present Eq. (|10aP : 
for 7// 0, the real part of the square brackets contains a term proportional to 
1h/E^) which contains an infrared divergence as q ^ 0. This problem disappears 
if the averaging over paths is performed more accurately, see MDSA-i2Q!. 



7. Discussion and Summary 

We have shown [in Apps. B to D, as summarized in App. A] that GZ's influence 
functional approach to interacting fermions is sound in principle, and that standard 
results from Keldysh diagrammatic perturbation theory can be extracted from it, 
such as the Feynman rules, the first order terms of a perturbation expansion in the 
interaction, and the Cooperon self energy. 

Having established the equivalence between the two aproaches in general terms, 
we were able to identify precisely why GZ's treatment of the Pauli factor {5 — 2p°) 
occuring Sr was problematic: representing it in the time domain as tanh[/io(t)/2T], 
they assumed it not to change during diffusive motion along time-reversed paths. 
However, they thereby neglected the physics of recoil, i.e. energy changes of the dif- 
fusing electrons by emission or absorption of photons. As a result, GZ's calculation 
yielded the result (jS'7j^)rw = 0. The ultraviolet divergence in (S'p^)rw, which in 
diagrammatic approaches is cancelled by terms involving a tanh function, was thus 
left uncancelled, and instead was cut off at w ~ I/toi, leading to the conclusion that 
I'^^iT is finite. 

In this review, we have shown that the physics of recoil can be included very sim- 
ply by passing from the time to the frequency representation, in which [5 — 2/5*^) is 
represented by tanh[fi,(e — uj)/2T]. Then {iSii)^vj is found not to equal to zero; 
instead, it cancels the ultraviolet divergence of {Si)rw, so that the total rate 
7v = 7ip + 7ip reproduces the classical result 7^^^, which goes to zero for T ^ 0. 
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Interestingly, to obtain this result it was sufRcient to average only over pairs of 
time-reversed paths; more complicated paths, such as represented by Hikami terms, 
are evidently not needed. (However, this simplification is somewhat fortuitous, since 
it occurs only when considering the unitary limit of the interaction propagator; for 
more general forms of the lat ter, the contribution of Hikami terms is essential, as 
emphasized by AAG and 

The fact that the standard result for "f^ can be reproduced from the influence 
functional approach is satisfying, since this approach is appealingly clear and simple, 
not only conceptually, but also for calculating j^. Indeed, once the form of the 
influence functional ([2]) has been properly derived (wherein lies the hard work), the 
calculation of {iSii + S'/)rw requires little more than knowledge of the distribution 
function for a random walk and can be presented in just a few lines [Sec|4]; indeed, 
the algebra needed for the key steps [evaluating Eq. ^ to get the first terms of ([9|) , 
then flnding (fTO|) and (fTT]) ] involves just a couple of pages. 

We expect that the approach should be similarly useful for the calculation of 
other physical quantities governed by the long-time, low-frequency behavior of the 
Cooperon, provided that one can establish unambiguously that it suffices to include 
the contributions of time-reversed paths only — because Hikami-like terms, though 
derivable from the influence functional approach too, can not easily be evaluated in 
it; for the latter task, diagrammatic impurity averaging still seems to be the only 
reliable tool. 
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Appendix A. Outline of GZ's Influence Fiinctional Approach 

Without dwelling on details of derivations, we outline in this appendix how the 
influence functional presented in Sec. O is derived. (A similar summary is contained 
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in a previous paper by this authoJ^^^^S however, it is incomplete, in that we have 
introduced important improvements since.) Before we start, let us point out the two 
main differences between our formulation and that of GZ: 

(i) GZ formulated the Coopcron propagator in terms of a coordinate-momentum 
path integral J jTyP, in which {6 — 2ff') is represented as [1 — 2no{ho)] — 
tanh(/io/2T), where the free-electron energy ho[R{ta), P{ta)] depends on position 
and momentum. This formulation makes it difficult to treat the Pauli factor with 
sufficient accuracy to include recoil. In contrast, we achieve the latter by using a 
coordinates-only version JVR^ in which exact relations between noninteracting 
Green's functions make an accurate treatment of the Pauli factor possible, upon 
Fourier-transforming the effective action to the frequency domain. 

(ii) GZ effectively performed thermal weighting at an initial time to that is not sent 
to — oo, but (in the notation of the main text) is set to to — ^■'"/2; with the latter 
choice, it is impossible to correctly reproduce the first (or higher) order terms of a 
perturbation expansion. GZ's claim in that they have reproduced these is 
incorrect (see end of App. C.3), since their time integrals have — t/2 as the lower 
limit, whereas in the Keldysh approach they run from — oo to +oo. We have found 
that with some (but not much) extra effort it is possible to properly take the limit 
to — !■ — oo, to correctly recover the first order perturbation terms [App. C.3] and 
to express the conductivity in a form containing thermal weighting in the energy 
domain explicitly in the form of a factor J (de)[— riQ(fe)]P^, where is an energy- 
dependent path integral, obtained by suitable Fourier transformation [App. C.4]. 

A.l. Outline of Derivation of Influence Functional 

Consider a disordered system of interacting fermions, with Hamiltonian H — Hq + 
He 

Ho = J dxip^x)ho{x)il){x) , (A.la) 

H, = ^Jdxidx2 ■4Hxi)i^{xi):Vi^' :^/'t(x2)^(a:2): (A.lb) 

Here J dx = J2a I '^'"'^ ^'^'^ = '^{I'l is the electron field operator for creating 
a spin-cr electron at position r, with the following expansion in terms of the exact 
eigenfunctions il)\{x) of ho{x) = + V.^^^{r) - pt: 

^P{x) = ^a(x)ca, [ho{x) - ^x\i:x{x) = 0. (A.2) 

A 

The interaction potential V^j* = V""*(|'''i — ''"2I) acts between the normal-ordered 
densities at ri and r2. The Kubo formula for the DC conductivity of a c?-dimensional 
conductor gives 



ctdc = -Re 



lim -j^Y] ( dX2jlV ■ 322' Jiv .22['^o) 
^o->o dix>o ^ J 



(A.3a) 
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Jii',22'(c^o) = / dh2e'^"''^e{h2)qn',22'] , (A.3b) 

J —oc 

C[ll',22'] = ^{[i'Htl, Xv)^p{ti, Xi),'ip'^ {t2, X2')4>{t2, X2)]) H, (A.3c) 

where jui = -^^(Vi — Vi') and a uniform applied electric field E{ijJq) was rep- 
resented using a uniform, time-dependent vector potential, E{ll!q) = iuJoA{uJo)- A 
path integral representation for C[ii'.22'] can be derived using the following strategy, 
adapted from (1) introduce a source term into the Hamiltonian, in which an 

artificial source field W2'2 couples to 'ip^ {t2, X2')il^{t2, X2), and write C[ii/_22'] as the 
linear response to the source field ■i;22' of the single-particle density matrix pu' = 
{tp^ {ti,xi')iip{ti, xi))h- (2) Decouple the interaction using a Hubbard-Stratonovitch 
transformation, thereby introducing a functional integral (. . .)y over real scalar 
fields Vp/B, the so-called "interaction fields", defined on the forward and backward 
Keldysh contours, respectively; these then constitute a dynamic, dissipative envi- 
ronment with which the electrons interact. (3) Derive an equation of motion for 
pXi' 1 the single-particle density matrix for a given, fixed configuration of the fields 
Vp/B, and linearize it in 'C2'2, to obtain an equation of motion for the linear re- 
sponse 5pXy{t) to the source field. (4) Formally integrate this equation of motion 
by introducing a path integral J 'D'{R) over the coordinates of the single degree of 
freedom associated with the single-particle density matrix SpYy- (5) Use the RPA- 
approximation to bring the effective action Sy that governs the dynamics of the 
fields Vp/B into a quadratic form. (6) Neglect the effect of the interaction on the 
single-particle density matrix whereever it occurs in the exponents occuring under 
the path integral J V R, i.e. replace there by the free single-particle density 
matrix 

p^. = {i^Hx,)4'{x^.))o = ^^I(x,>A(x.)no(a) , (A.4) 



where thermal averaging is performed using (6)0 — Tre[^^^°0]/Tr[e^'^^"]. (7) Per- 
form the functional integral {. ■ ■)v (which steps (5) and (6) have rendered Gaussian) 
over the fields Vp/s', the environment is thereby integrated out, and its effects on 
the dynamics of the single particle are encoded in an influence functional of the 
form e~'''^^+'^^). The final result of this strategy is that J22' • jiv C\^\v .22'\ can be 
written as [cf. (11.49)] 

/■If /-Ik 



dx2 322' ■3ivCyxv.22']^ Idxo.roBPLdsi f ^'(^) (A-5) 

j(i2j-J-(i2jlj-(ti)e-'^^«+^^l^*^'*"^/' 



4' 
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where ^ ^ is used as a shorthand for the following forward and backward 

path integral between the specified initial and final coordinates and times: 



P'H^(tf)e-'^^o"(*f.*f)M... (A.6) 



R^[tf)= 

The complex weighting functional e^'^'^o^""^'?-' occuring therein involves the action 
for a single, free electron. Expression (jA.sp has a simple interpretation: thermal 
averaging with at time io (for which we take the limit — oo) is followed by 
propagation in the presence of interactions (described by e~[*'^"+'^^l) from time 
up to time ii, with insertions of current vertices j{t2a) time ^2 on either the 
upper or lower Keldysh contour, and j{ti) at the final time ti. 

For the purpose of calculating the Cooperon propagator, we now make the follow- 
ing approximation in Eq. (jA.Sp [referred to as "approximation (ii)" in App. B]: For 
the first or second terms, for which the current vertex occurs at time on contour 
a — F or B respectively, we neglect all interaction vertices that occur on the same 
contour a at earlier times or ^4- e [^01^25]; however, for the opposite contour 
containing no current vertex, we include interaction vertices for all times € [^Ojii], 
with to —00. [This turns out to be essential to obtain, after Fourier transforming, 
the proper thermal weighting factor [— no(fe)] occuring in Eq. (fTa|) . see App. C.4.] 
The rationale for this approximation is that, in diagrammatic language, this ap- 
proximation retains only those diagrams for which both current vertices J22' and 
jii' are always sandwiched between a G^- and a G'^-function; these are the ones 
relevant for the Cooperon. The contributions thereby neglected correspond to the 
so-called "interaction corrections" . [If one so chooses, they latter can be kept track 
of, though.] 

This approximation (ii) is much weaker than the one used by GZ at a similar 
point in their calculation: to simplify the thermal weighting factor describing the 
initial distribution of electrons, namely to obtain the explicit factor po in Eq. (49) 
of GZOS"^, they set to — > t2 (their t' corresponds to our t2), and thereby perform 
thermal weighting at time t2, instead of at —00. As a consequence, in their analysis 
all time integrals have t2 as lower limit, which means that (contrary to their claims 
in GZod^ they did not correctly reproduce the Keldysh first order perturbation 
expansion for C[ii'_22']5 in which all time integrals run to —00. A detailed discussion 
of this matter is given at the end of App. C.3. [Contrary to our initial expectations, 
but in agreement with those of GZ, it turns out, though, that the choice of to 
does not have any implications for the calculation of t^, which does not depend on 
whether one chooses to — t2 or sends it to — cxd.] 

Having made the above approximation (ii), the effective action {iSu + Sj) oc- 
curing in Eq. (jA.5|) is found to have the following form (we use the notation iSfiJ Sj 
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to write two equations with similar structure in one line, and upper or lower terms 
in curly brackets refer to the first or second case): 

[iSR/Si]{h,to) 
{^L^/L')3,4s 

Here ~Su = and = -tj„,,<(izj-r|'(tj^,)) 

are the standard retarded, advanced and Keldysh interaction propagators. For each 
occurrence in Eqs. (|A.8|1 of a pair of indices, one without bar, one with, e.g. 4^ and 
4a, the corresponding coordinates x1 and are both associated with the same time 
^4, and integrated over, J dx^dx'^, in the path integral J V^R). (This somewhat 
unusual aspect of the "coordinates-only" path integral used in our approach is 
discussed in explicit detail in App. D.4; it is needed to account for the fact that the 
density-matrix pP is non-local in space, and arises upon explicitly performing the 
J "DP momentum path integral in GZ's formulation.) The Su functions on the right 
hand side of Eqs. (|A.8p will kill one of these double coordinate integrations at time 
U. 

Eqs. (|X7| and fO]) are the main result of our rederivation of the influence 
functional approach. They are identical in structure (including signs and prefactors) 
to the corresponding expressions derived by GZ (Eqs. (68) and (69) of GZ9jl^, as 
can be verified by using the relations 

to relate our interaction propagators Cy to the functions Rij and ly used by GZ. 
However, whereas Eqs. (68) and (69) of GZ9P are written in a mixed coordinate- 
momentum representation in which it is difficult to treat the Pauli factors {S — 2p°) 
sufficiently accurately, our expressions (jA.8[) are formulated in a coordinates-only 
version. Formally, the two representations are fully equivalent. The key advantage of 
the latter, though, is that passing to a coordinate- frequency representation (which 
can be done before disorder averaging, allows us to sort out the fate of {S — 2p°), as 
discussed in Sec. AI3] [and extensively in App. B.6.2]. 



ti 

to J to 



J2 * / *(ii''/i')3.4., 



;;0 



[S - 2p 
'^3M ~ ^^^^ 



3b4:f 3b3b ' 



'^4i54f 



T2« (^34 03^3^ 1-4^3^ 



(A.7) 
(A.8a) 
(A.8b) 
(A.8c) 
(A.8d) 
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A. 2. Cooperon Self Energy before Disorder Averaging 
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Fig. Al. First order contributions to the irreducible self energy of the Cooperon, illustrating 
Eqs. l lA.lOll . The arrows associated with each factor Gij or Cij in Eqs. l lA.lOl l are drawn to point 
from the second index to the first {j to i). Filled double dots denote the occurence of a factor 
{& — 'ip)4,p4,p on the upper contour or (5 — 2p)4^4^ on the lower contour. Bars on filled dots are 
used to indicate the barred indices to which the interaction lines depicting C^j are connected. Both 
filled and open single dots indicate a delta function 5; the open dots stand for delta functions that 
have been inserted to exhaust dummy integrations, as discussed after Eqs. l|A.8| l [and, in more 
detail, in Sec. 6.1]. The diagrams in (b) and (c) coincide precisely with those obtained by standard 
Keldvs h di agrammatic perturbation theory for the Cooperon self energy, as depicted, e.g., in Fig. 2 
of Ref. liSI, (There, impurity lines needed for impurity averaging are also depicted; in the present 
figure, impurity averaging has not yet been performed.) 



From the formalism outlined above, it is possible to recover the standard results 
of diagrammatic Keldysh perturbation theory, before disorder averaging, by expand- 
ing the path integral (jA.5[) in powers of the effective action i [iSn+Si). For example, 
using Eqs. (jA.Sp [and being sufficiently careful with signs, see App. B.6.1] one read- 
ily obtains the following expressions for Cooperon self energy E^/^ = Saa' ^^a'^' 
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summarized diagrammatically in Fig. lAll 



[t'^/jY^'^ = (G^/^)3-5-G43j£^7£^)3-^- , (A.lOa) 



= (G^/^f ^^^Gt3.(^^/^^^)3/^ (A.lOb) 
(S^/;)'^'^ = G«-^-^.((5KM)_^^^(^~A/i^~K) 3.^ (A.lOc) 

\ / 4:BOB ^ 

(S^//)'^'^ = G«-3-^-(G^M)4.3.(^-^/^^)4.3. . (A.lOd) 

To obtain this, we exploited the fact that every vertex occuring in the effective action 
is sandwidched between retarded propagators if it sits on the upper contour, and 
advanced ones on the lower contour. The Keldysh functions arise from using some 
exact identities, valid (before impurity averaging) in the coordinate-time represen- 
tation: depending on whether a vertex at time t/^i sits on the forward (time-ordered) 
or backward (anti-time-ordered) contour (a! = F/ B)^ the factor [6 — 2fP)L^/^ oc- 
curing in L^^, is sandwidched as follows (on the left hand sides below, a coordinate 
integration J dx4^, over the un-barred variable at vertex 4 is implied) : 



G«4.(<5-2p°)4,4. /:f4,Gf^^.^- G«4^(£--(i)£^i^(<i)Gi;^.^(£-), (A.Ua) 



i^-mis^sGt.u ^-Gti^e)Ct^,{u)Gf^,^{e-u) (A.llb) 



The left- and right-hand sides are written in the time and frequency domains, re- 
spectively. To obtain Keldysh functions from the left-hand side expressions, we 
exploit the fact that the upper or lower contours are time- or anti-time-ordered to 
add an extra —G^^^ = 0, and then exploited Eq. ([5a|) to obtain a factor ±G^ (see 
Sec. B.6.2). 

A. 3. Thermal Averaging 

It remains to figure out how the thermal weighting in Eq. (|la[) can be derived from 
our general path integral expression Eq. (|A.5[) . This is a standard, if nontrivial, 
excercis e in Fourier transformation, carried out (along the lines of a similar analysis 
by kk\^) in App. C.4. The result is an equation for the conductivity similar to 
but more general than Eq. |T]), with J^dr P21, cffiT) replaced by J^dTi2 '^('''12), 
involving a slightly more complicated path integral [Eq. (jC.2ip ]. defined as 

P^lr'{T)= dfe"'W M I?'(i2)e-['^«+'^^l/'' . (A.12) 

Note that the duration of the forward and backward paths differs by a time f, in 
contrast to the path integral (jlbp used in the main text. The combination J de J df 
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of integrals from Eqs. (fTa|) and (jA.12p have the effect of fixin J2l] ^j^g average energy 
of the forward and backward trajectories to be close to the Fermi energy, with 
energy spread of roughly zLT (see App. IC.4l for a detailed discussion). This energy e 
is the same as the one that in perturbative calculations shows up in the tanh[?i(£ — 
[D)/2T]-factors of the Keldysh electron Green's functions {e — cu), which play a 
role in determining the phase space available for electrons to get scattered upon 
absorbing or emitting a noise quantum. In App. [C^ we argue that the simplest way 
to keep track of this in the influence functional approach is to replace Eq. (jA.12[) 
by Eq. (|lb[) . which mimicks the effect of the former's integral / dfe^^'^ by using (i) 
forward and backward paths of equal duration r and (ii) an effective action whose 
time integration boundaries are fixed at ±t/2, but which depends explicitly on the 
average propagation energy e [via Eqs. ([2), (|4ep . or equivalently Eqs. (|B.94p . (IB.97P ]. 

Note that GZ's approach in effect employs the same simplification, since they 
likewise have no J dfe'^"^ integral and use forward and backward paths of equal 
duration r. Their effective action depends on the average energy e, too, via the 
tanh[?ie/2T]-factor in their Sr. However, lacking the —Co recoil shift, their tanh- 
terms turn out to yield zero after averaging over random walks, so that (iS^^^)rw — 
0. 

A. 4. Perturbative vs. Nonperturbative Methods 

We conclude this overview-style appendix with some general comments on whether 
it is sufficient to calculate perturbatively, as we contend (in agreement with 
otherjn™) 

, or whether a truly nonperturbative approach is needed, as GZ have 
argued in We have made an effort to keep the discussion as nontechnical 

as possible and accessible to casual readers that have not studied App. [B] in detail, 
although we will on occasion refer to results from the latter. 

In GZ's influence functional approach, the decoherence time is deflned as the 
scale at which the function F{t) = j^{iSR + S'/)rw, which in their theory is linear 
in the interaction propagators R/I, becomes of order one. This means that is 
the crossover scale between the regimes where perturbation theory is rigorously 
valid or breaks down, F{t) ^ 1 or ^ 1, respectively. To determine this scale, 
we contend that it is sufficient to calculate F(t) perturbatively (assuming, strictly 
speaking, F{t) <C 1), and then to enquire for what time the perturbative result 
so obtained ceases to be small, setting F{Tip) ~ 1. (This is analogous to the fact 
that the crossover scales Tk or Tc, the Kondo temperature in the Kondo problem 
or the critical temperature in the theory of superconductivity, can be calculated 
perturbatively as the scales where perturbation theory breaks down.) An accurate 
knowledge of F{t) for t > would be needed only if we desired to accurately 
include exponentially small (e^^*-^-* ^ 1) contributions to weak localization, which 
is usually deemed not worth the effort. (In contrast, for the Kondo problem or 
superconductivity, nonperturbative treatments are worth the effort, because the 
phenomena of interest become strong in the nonperturbative regimes.) 
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GZ have argued in GZOG^ that a perturbative treatment of weak locahza- 
tion is insufhcient, because according to them it fails to disentangle the effects 
of preexponent and exponent in an Ansatz for the Cooperon of the general form 
C(t) — A(r)e~^^'^-' : when this is expanded in powers of the interaction, both A and 
F contribute to the first-order term C*^^-' . The influence functional approach avoids 
this problem by very naturally generating a general expression for the function F 
in the exponent - which in GZ's approach turns out to be linear in the interac- 
tion propagator [Eq. or Eq. (|B.84p ]. However, the problem of disentangling the 
exponent from the preexponent is easily avoided in the diagrammatic approach, 
too, by calculating not the Cooperon itself, but its self energy, to linear order in 
the interaction; Fourier transforming the resulting Cooperon C{uj) into the time 
domain, this automatically yields an expression of the form ^(r)e~^*-^\ again with 
F linear in the interaction propaga tor . [The prefactor arises from wave-function 
renormalization effects, see DMSA-l]2S, Eq. (14a).] Since both the influence func- 
tional and diagrammatic strategies yield results for which the exponent F is linear 
in the interaction (and contains contributions with a similar coth -I- tanh structure) , 
it is reasonable to expect that if both approaches are implemented with sufRcient 
care, their answers for F should agree completely. 

They do agree, in fact, if the recoil-incorporating effective action proposed in this 
work and featuring tanh[^(£ =F p)/2r] -factors is used. (This agreement is demon- 
strated explicitly in DMSA-Il^^.) But they differ if GZ's procedure is followed with- 
out modification, leading to their no-recoil tanh[fe/2r]-factors. It is important and 
instructive, therefore, to identify at which point of the derivation of the influence 
functional approach the need for a modiflcation of GZ's approach flrst manifests 
itself. We shall now argue that this point is reached when the order in which two 
distinct averaging procedures are performed, over paths R"" and fields V, is tacitly 
interchanged, an aspect that has not been emphasized in the preceding sections. 

To be concrete, let us focus on an intermediate stage of GZ's first principles 
calculation of the weak localization contribution (t^q to the conducitivity. Following 
the enumeration of steps used in App. lA.ll p.fTSl this stage is reached after steps (1) 
to (6) [or according to the enumeration of App. IB.41 p. [36l after steps (A) to (G)], 
resulting in the following expressions [the first of which corresponds to Eq. (jB.55a[l ]: 



CTDC.roal^^^ J dX2 jll' ■ j 22' J'u' ,21' i^) ^ (A. 13a) 

J[2'.2i'io) = J[3':2i'm ((e^(^^^-^v)/'') ) , (A.iab) 

\ \ / cqp/ V 

Sy"" ^ dhhl"'{h,R^/^{H)) , (A.13C) 

^(2'?2i'(0) {de)[-n'{he)] dfe^^ (A.lSd) 



cqp 
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The correlator J^^i 21' (0) originates from C[ii/.22'] of Eq. (jA.3c[) . It has here been 
expressed [starting from Eqs. (|B.49p to (jB.52p . and using the results of Eqs. (|C.14[) . 
(|C.19a|) and (|C.21[) ] as a double average ((. . ■)cqp)y over a pair of phase factors 
gH-Sv-Sv)^ which describe the influence of interactions, represented by fluctuating 
flelds Vp/B(t3,r3), on a pair of closed quantum-mechanical paths (cqp). [The de- 
tailed form (|A.13c[) of the phase factors follow from Eq. (|B.56p . with hy given by 
(|B.36bP : see also Eq. (|B.58|) . and the discussion thereafter]. Eq. (jA.lBbp instructs us 
to first pick out a specific configuration of the fields Vf/s(^3j '"3): then to calculate 
the average (. . .)cqp of this phase factor over all closed quantum-mechanical paths 
with boundary conditions specified in Eq. (jA.13dp [as obtained from Eq. (|C.21[) ]. 
and to evaluate the average over all field configurations in the end. Thus, for a given 
V, the set of paths making the dominant contribution will depend on V. 

Now, the next step of GZ's strategy [step (7) according to App. lA.li or step 
(H) according to App. IB.4j . is to perform the average (. . .)y over the interaction 
fields. To carry out this step, GZ (tacitly) interchange the order of averages [as do 
we in App. IB. 575] . in effect replacing Eq. (jA.13bp by 

\ \ / V / cqp \ / cqp 

SMR"] = ^{(Sv - S^)iSv - Sv))v = + Si ■ (A.14b) 

Eq. (|A.14a|) instructs us to first pick out a specific pair of paths R^^^{t:}), and then 
to calculate the influence functional (^e^'^^v-s^/)/h^^ which describes how the chosen 
pair of paths are effected on average by interactions. Within the RPA approxima- 
tion, the (. . .)y average can now be done exactly, yielding an effective action 5off 
that is linear in interaction correlators R/I [Eq. (|B.84[) ]. The sum over all closed 
quantum paths is to be performed at the end. 

Now, this seemingly innocuous change in the order of averages is without con- 
sequence only if both averages are performed exactly, as is possible in an order-for- 
order perturbation expansion (or, to all orders, for exactly solvable models such as 
the Caldeira-Leggett model). However, this is not the case in GZ's theory (or our 
version thereof), which proceeds to use the semiclassical approximation of replacing 
the sum over all closed quantum paths by a sum over only the saddle point paths 
that extremize the action. In principle, this can be done in at least two different 
ways, which we indicate schematically as follows: 

/ cqp \ / bare 

-pe«Kl\ 'dcally / ^^^^^^^ 
I cqp \ / dressed 

Here the subscripts bare/dressed indicate that the sums over the paths on the right 
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side of Eq. (|A.15|) are taken only over those paths, with boundary conditions as 
specified in Eq. (jA.13dp . which extremize the bare action 5*0 = Sq — Sq ("bare" 
paths, Eq. (|A.15ap . used by GZ), or the full action ^tot = + iSeS ("dressed" 
paths, Eq. l|A.15bp . discussed below). On the far right of Eq. (|A.15a|) we indicated a 
further (uncontroversial) approximation, used by GZ and others when an exponen- 
tial is to be averaged over bare paths, namely to lift the average into the exponent. 
In practice, [e.g. Sec. [4j, the averages (. . .)barc on the r.h.s. of Eq. (jA.lSaP are re- 
placed by ((. . .)bare)dis ^ (• ■ ■)rw, whcrc the latter average is over diffusive random 
walks (rw) with appropriate boundary conditions. In other words, (. . .)barc is ap- 
proximated by considering only semiclassical trajectories in a disordered potential 
landscape (while fluctuations about these semiclassical paths are neglected), and, 
after after implicit disorder averaging, these semiclassical trajectories are treated as 
random walks. 

Ideally, one would of course prefer to average over dressed paths [Eq. (jA.lSbp j , 
which "know" about the effects of interactions due to the role that iScS plays in de- 
termining the saddle point trajectories. (Even more ideally, one would also take into 
account fluctuations about these dressed paths). In such a calculation, the iSj term 
in iScS would cause the dressed paths -RJjrcssod acquire an imaginary component 
(we thank Igor Gornyi for alerting us to this fact), implying that the contributions 
of the two terms in {iSr + '5^/)[-Rdrossod] partially cancel, even though both Sr 
and Si are purely real functionals of their arguments (GZ overlooked the possibility 
of such a partial cancellation, because they considered only bare paths, see below. 
Marquardlff^ has illustrated how such a partial cancellation occurs in the Caldeira- 
Leggett model). Note that such a dressed-path procedure would require only that 
5*101 [-Rdrcsscd]/^ > 1, and would not require 'S'cff [-Rdrcssodl *o ^e small. Indeed, 
its results would be nonpcrturbativc in the interaction correlator R/I, since ScS, 
though linear in i?//, is a non-linear functional of -RJ^i-egs^d ' which itself is nonlin- 
ear m R/I (as illustrated explicitly in the Caldeira-Leggett model, where all these 
nonlinear functions can be evaluated explicitly). 

However, in the present theory, using fully dressed paths is not technically fea- 
sible. Therefore, GZ made the standard and seemingly natural choice of averaging 
purely over hare paths. Indeed, they write (just before Eq. (61) of GZ9{P|): "In 
the zero-order approximation one can neglect the terms Sr and Si describing the 
effect of Coulomb interaction" so that "the path integral is dominated by the saddle- 
point trajectories for the action 5*0" . In other words, bare paths don't "know" about 
the interactions at all. Consequently, GZ used an effective action S'^^ obtained by 
treating the Pauli factor 5—2pm. S^s as time- independent (arguing that the energy 
argument of the corresponding Fermi function is conserved during propagation), 
essentially replacing it by tanh[fe/2T]. [See p. 9205 of "n depends only 

on the energy and not on time because the energy is conserved along the classical 
path".] 

Once the approximation of using purely bare paths has been made, the ef- 
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fective action S^g[R^^^J is linear in the interaction propagators (since -Rtarc 
now independent of the interaction). This imphes in our view that GZ's results 
are purely perturbative. GZ dispute this characterization, calling their approach 
"nonperturbative" because in their view it does not require S^e'[Rbaic\ ^ only 
S^siRharcl ^ ^o[Rha.rc\- ^c disagree, contending that GZ do need the former con- 
dition, because without it, their use of purely bare paths would not be be justified: 
a semiclassical treatment requires the evaluation of the action S'tot[-R"] to be ac- 
curate to within ft, implying that the effects of SeS on the semiclassical paths can 
be neglected only if S^^iR^^y-J ^ ^- [Note, also, that an approach that rehably 
evaluates (Scs) in the regime where the result is < ft would yield the function F{t) 
in the regime where it is < 1, which is entirely sufficient to reliably extract r^, as 
argued in the second paragraph of this subsection.] 

While it is a matter of somewhat empty semantics whether an approach using 
purely bare paths can be called nonperturbative or not, the validity of such an 
approach can be subjected to a hard test: does the result which this approach 
produces for F(t) after the average (. . .)baro has been performed agree, in the regime 
F{t) <C 1, with that obtained from Keldysh perturbation theory? (GZ's claim in 
GZoP that their approach agrees with Keldysh perturbation theory, is true only if 
the perturbation expansion is performed before averaging over paths, see Secs. lB.6.1l 
and lC.3( l. The answer is no: The perturbation expansion obtained by expanding the 
first-principles expression (|A.13p in powers of {Sy — Sy)/h shows unambiguously 
that the paths arising in the perturbation expansion, do know about the interactions 
(in contrast to bare paths): energy conservation induces recoil at each interaction 
vertex, so that the electron frequencies incident and leaving a vertex differ from 
each other, e vs. e =F w, in a way relevant for the Pauli factors, which depend on 
e =F w. This effect is negligible for retarded and advanced electron propagators, 
which depend only on the combination h{e =F cD) — ± ih/Tc\, with = P'^ /2m 
[Eq. (|F.2p ]. since there energy shifts by ft|a)| < ft/rd are negligible. However, it is 
not negligible for the Keldysh propagator, which contains fermion functions of the 
form [e^'''-(^^'^)-^p]/'^ + 1]-^ [Eq. (jB.47bp ]. in which the largeness of fie is cancelled by 
that of e^^. Thus, interaction events with recoil energies of order fi|(D| > T strongly 
change the value of the Fermi function which specifies the phase space available 
for a given transition. Since these recoil effects are present in the original order 
((• • ■)cqp)v of doing the average but absent in the switched order (e~'^<=ff ^'')barc if 
GZ's version of the effective action is used, something is clearly amiss in the latter 
approach. 

The main assertion of our own work is that for the purpose of describing deco- 
herence in weak localization, recoil effects can be taken into account in the influence 
functional approach provided GZ's use of bare paths is supplemented by the use of 
an effective action that keeps track of recoil (rec): 



JvD 



3 h "-^cff [^barc] 



(A.16) 
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[In practi ce, w e perform the averages over bare paths on the r.h.s. the same way 
as GZ do,l^^ i.e. using an average over semiclassical diffusive random walks (and 
neglecting fluctuations about these), ((. . .)barc)dis (• • •)rw] Here S^"^ is the effec- 
tive action obtained by representing {S — 2p) by tanh[/i(e =F aj)/2T]-factors [as made 
explicit in Eq. jEMJ]. The result for [-Rb^re] so obtained [Eq. or Eq. (|R94l) 
with (jB.97p ] is linear in the interaction propagators R/I, just as GZ's effective ac- 
tion is, but in contrast to the latter, an expansion of e~'i^'^<=ff ^•">''= to first order in the 
exponent yields results consistent with the Keldysh perturbation expansion also if 
the average over paths in (S'^g'^)barc is performed explicitly first and the exponen- 
tial expanded only thereafter. Moreover, the results for F{t) so obtained agree fully 
with those from a diagrammatic Bethe-Salpeter calculation of the Cooperon (see 
DMSA-II^'-'). A crucial ingredient for ensuring this agreement is that the ultravio- 
let divergencies arising in each of the two terms in {{iS]^^ + 5}'''^)[-Rbare])barc cancel 
each other [Sec. H] . This cancellation is possible because the functional S'^'^ [R^^j.^] , 
despite using only bare paths, is not purely real, thereby capturing an essential 
feature of [^Srcssod] that is not present in ^^^[-Rbarc]- 

Since GZ contend that their approach is nonperturbative, they reject arguments 
based on perturbation theory, defending their use of purely bare paths by evoking 
only the standard semiclassical approximation. But the need to keep track of recoil 
arises within the latter framework, too, in a way very similar to that described 
above: The standard condition for the validity of the semiclassical approximation 
is that the propagation energies and momenta of the quantum particle that is to 
be described semiclassically should be much larger than the typical frequencies and 
wave numbers characterizing the potential landscape which it is moving in, so that 
the latter appears "smooth". If one were to consider a single noninteracting elec- 
tron propagating with energy he ep through a disordered potential landscape, 
this implies the conditions ep 3> h/rd (or kp 3> l/Zci), which certainly are satisfied 
in the regime of weak localization. However, GZ's theory for an electron propagating 
through and interacting with a Fermi sea of other electrons shows that the propa- 
gation energy enters not only in the free part of the action, but also in the Fermi 
function [e^'''^^^^)/^ + 1]^^ arising from the Pauli factor S — 2pm Sr. To ensure that 
this factor is treated accurately, the standard semiclassical condition 3> h/rd 
evoked by GZ is not sufficient, since inside the Fermi function the largeness of he 
is cancelled by the largeness of ep- Thus, interaction-induced changes in e of order 
Lu < 1/toi will produce strong changes in the Fermi function between ~ and ~ 1. 
These changes need to be kept track of. As argued above, this can be accomplished 
by the recoil contributions ^lo in our tanh[/i(e =F a;)/2T]-factors. 

Appendix B. Derivation of Influence Functional Approach 

In this appendix, we rederive the influence functional approach of GZ, with the 
aim of establishing clearly (i) how far it can be taken without any approximations, 
and (ii) what the consequences are of the approximations that they eventually do 
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make. We generally follow the strategy they have chosen to use, but the details of 
our notations and derivations deviate from GZ's whenever we believe that greater 
compactness, clarity or generality can thereby be achieved. The most important 
difference is that instead of using the coordinate-momentum path integral / 'D{RP) 
of GZ, we use a coordinate-only version / T^'R, since this enables the Pauli factor 
to be treated more accurately. 

The outline of this appendix is as follows. After a summary of our notational 
conventions. Sees. IB. II to IB. 31 define the model and decouple the interaction using 
a Hubbard-Stratonovich transformation within a Keldysh framework. In Sec. IB. 41 
we summarize GZ's procedure for deriving their influence functional approach, and 
in Sec. IB. 51 repeat their steps in explicit detail, though with some changes. Finally, 
Sec. IB. 61 establishes a link between the influence functional so derived and Keldysh 
perturbation theory, and discusses the fate of the Pauli factor. 

Notational Conventions 

We begin by summarizing, for ease of reference, some notational conventions to be 
used throughout: We shall use the shorthands x = (r, cr) for electron position and 
spin, and J dx = J2a J^*"- Operators will generally carry hats (e.g. Hq), and the 
subscripts 5, H and I will distinguish operators in the Schrodinger, Heisenberg or 
interaction pictures, respectively. For c- number flelds, the shorthand Vi = Vi{ti) = 
V{ti^ ri) will often be used, i.e. the time argument, when not displayed explicitly, will 
be understood to be ti. c- number functions depending on two different coordinates, 
i.e. coordinate-space matrices, will generally carry tildes [e.g. pij — p{xi,Xj)], and 
their Fourier transforms w.r.t. — rj will carry bars, e.g. 

p{Ry,p) = j dnje-'P""'' p{Ry + \ry,Ry - \ry) , (B.l) 

where Rij and rij will generally denote center-of-mass and relative coordinates, 

R^j = [r, + r.j)/2 , = r, - . (B.2) 

We do not set h — 1 but display it explicitly throughout. Hence, the variable p in 
Eq. (|B.1[) (and likewise fe, q below) denotes a wave-number, with units of inverse 
length, not a momentum; the corresponding momenta will always be denoted by 
capital letters: 

P ^hp K ^hk, Q^hq. (B.3) 

For correlation functions, the shorthand Gy = Gij{tij) = G{tij;Xi,Xj) will often 
be used, i.e. the time argument, when not displayed explicitly, will be understood 
to be tij — U — tj. [For the step function, we use By = 0{tij).] The corresponding 
frequency Fourier transform w.r.t. ty will be denoted by 

G,(%) = J {doj)e-^-'^^Gy{u) , j {dLu) = / £ , (B.4) 
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where uj has units of inverse time. If coordinate-space subscripts are not dis- 
played explicitly, they are understood to be summed over, e.g. [G{t)G{t')]ij = 
J dxkGik(t)Gkj{t'). We distinguish forward and backward parts of the Keldysh con- 
tour by an index a = F, B [GZ use a — 1,2 instead], and use boldface for matrices 
in Keldysh space, e.g. Gij. 

A pair of indices such as li', appearing once without prime and once with, will 
denote independent coordinates Xi and Xi' referring to the same time (i.e. U = tii 
is to be understood), which are, however, to be set equal at the very end of the 
calculation, after being differentiated upon, i.e. 

(V, - V,)p.,e ^ [(V. - V.rP^^'U=^, ■ (B.5) 

We shall often encounter double summations over coordinates referring to the same 
time. For such coordinates we shall use the index pair ii, one without bar and one 
with, take it to be understood that = ti, and denote the double summation by 

j dxij = ^ j dridrj . (B.6) 

When taking the limit of infinite volume, we shall use the shorthand notation 
(dp) = f = hm — y , (B.7a) 

^ ^' J (2tt¥ VoWoo Vol ^ ' ^ ' 

S'-'^\p^p)= lim Vp'VoI, (B.7b) 

Vol — *oo 

so that J{dp)S{p) = 1, i.e., S{p) equals {2tt)'^ times a d-dimensional Dirac delta 
function. If the integrand under J (dp) depends only on the energy 
and if it decays at least as fast as 1/^p for ^ oo, we shall use J (dp) v j d£^p. 
Here v denotes the density of states per spin at the Fermi surface, which in d = 3 
or 2 dimensions is given by 

d-2 , I 1 



27rn2 \ ] 2£p 



(B.8 



where (n) = /^^^^(dfc) = ^ (^) is the average electron density per spin. The 
'purely 1-dimensional case d = \ will not be considered here; nevertheless, d — 
or 2 of course include the case that a sample is quasi 1-dimcnsional, in the sense 
that only one of its dimensions is larger than the phasebreaking length, L > Dt^, 
where D = UpTd/d is the diffusion constant. 

For quasi-d-dimensional diffusion, the actually measured (bare) DC conductivity 
cTj; is related to the Drude conductivity o'd™"^'^ ~ 2e^vD by an extra factor a''"'^, 
which accounts for the sample's transverse directions along which motion is not 
diffusive (d = 3 or 2 is the actual dimension of the sample, d— 3, 2 or 1 the effective 
dimension for diffusive motion). Hence, it is customary to define [cf. AAl412J, after 
Eq. (5)] cTj — (Tp^'^'^a'^^'' as the conductivity per unit length and unit area of 
a 3D sample (for d ~ '6), or the inverse square resistance of a film (for d = 2), 
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or the inverse resistance per unit length of a wire (for d = 1). Likewise, the weak 
localization correction to the conductivity is often expressed in terms of these actual 
conductivities by defining = c^c a'*^''. 

The fact that the weak localization correction is small compared to the Drude 
term is often made explicit by writing the prefactor of the Cooperon term as 
crQ™'^°/gj(Lif) [see Eq. (|C.9p ]. where Lh = \^Dth is the magnetic length and 

g3{L)^{h/e')a^L'^-\ (B.9) 

is the so-called dimensionless conductance, defined as the conductance, in units of 
e^/fi, of a rectangular (d-dimensional) block with volume a'^^'^L'^, measured along 
one of the "long" directions (of length L). 

For good conductors, gj(L) {tt^-'' /d){akF)'^^^{kikFy-\L/lci)^-^ is large 
whenever L is large: we may assume Icikp 3> 1 and akp ^ 1 throughout, thus 
for d = 3 or 2, any length L > l^j implies gg{L) ^ 1; for d = 1 the function 
gi{L) likewise starts out being ^ 1 for L ~ Id, but decreases with increasing L; 
nevertheless, it reaches gi{L^) ^ 1 only when L exceeds the very large length scale 

= (aA;F)''"^/ci > 

B.l. The Model and Kubo Formula 

Following GZ, we consider a disordered system of interacting fermions, described 
by the Hamiltonian H ~ Hq + Hi , where 

Ho = J dxi,lix)hoix)iJ;six) , ho{x) = + V;,„p(r) - /i , (B.IO) 

r 

H^^- dxidx2 -.hiis- Vlf ■.h22S-, (B.ll) 

■ n^Js■■ = fiijs - {ntjs)Q , "jjs = V's(a^j)V's(^») ' (B-12) 

(0)o = Tr{6po}, Po^e-^"°/{TTe->'"°}. (B.13) 

^mp('") is the disorder potential. Wc shall assume that the interaction potential 
^12* — V^™*(|ri — r2\), which guarantees that its Fourier transform in d effective 
dimensions, Vj'^*(q) — Vj"*(|<z|) = a^"'^/ d'^r e^^^"^V"^^{r), is real. For example, for 
the Coulomb interaction, they are given by V72* — \ri-r2\ ' ^^'^ 

Vr^ = %, V^'' = '-^, ^«=a^ln(,-V). (B.14) 

g \q\ 

Eq. (|B.12p corresponds to a normal-ordering prescription which subtracts the ex- 
pectation value w.r.t. the free density matrix po- The second-quantized electron field 
il^g{x) = "05 (t", a) (in the Schrodinger picture) destroys a spin-a electron at position 
r, and can be expanded as follows in terms of the exact eigenfunctions ^\{x) of 
hQ{x), with eigenvalues 

[ho{x)^i)]^x{x) = i). (B.15) 

A 
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The current density operator has the form 

^2 



Jiv 



(B.16) 



where A is the vector potential, jiy = (— ie?i/2m)(Vi — Vi'), and the convention 
of Eq. (jB.Sp was used for the indices ii'. An external electric field, E{t2,r) — 
— V14xt(t2, r) — dt2A{t2,r), switched on at time I'q, is described by the perturbatior0 



^ext(i2) = 0{t2 - Q j dx2 hf\t2,X2) , (B.17a) 

hT{t2,X2) = /i°22' ^22'S > /»22' = e^xtfe, - ^(^2, Ta) • ^22' • (B.17b) 

According to the Kubo formula, the linear response of the current density to this 
perturbation is 



-—Mti^'ri){n^^g) -i / dt2 / rfa;2 ^22' Jii' C[ii',22': 



(B.18) 



The first term of Eq. (jB.18[) is the diamagnetic contribution, (n) = 21/6^ / d being the 
average electron density per spin [cf. Eq. (|B.8|) ] . In the second term, the correlator 



C[ii 



(B.19) 



is to be evaluated with i?oxt set to zero, where Bnit) — e^-^*'^^Bse^^'^*'^^ describes 
time evolution in the Heisenberg picture, and thermal averaging is defined by (O) — 
Tt{0 ph}, where pn = e~^^ / {Tre~^^} is the full equilibrium density matrix. 

The DC conductivity is defined via the low-frequency limit of the current re- 
sponse to a spatially homogeneous applied AC field E{t2) = J {dojo)e~^'^°*^ E{(Zio)- 
For a d-dimensional isotropic sample it can be written as 



CDC 



lim — - 



(B.20) 



where E{llJq) can be represented by either of the choices (related by a gauge trans- 
formation) 



A{ujo) 



E{u;o) 



-r2 ■ E{uJo) 



lujo 



(B.21a) 
(B.21b) 



GZ use choice (jB.21ap (but note our footnote [aj), A AG use choice (|B.21bp . Taking 
the limit —00, one then readily finds from Eqs. (jB.20[) and (jB.lSp that ctdc 



We use e < for the electron charge, as do AAcffZI, whereas GZ use — e < 0, hence our 
potentials are related to GZ's by a minus sign: eV^^" = —eV^^, and likewise eVj!'^''" = —eV^^, 
gyhcrc _ _^yGZ f^j. ^j^g potentials introduced in Eq. l lB.28dl l below. 
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can be written in either of the following forms, depending on whether the electric 
field is represented using a scalar or a vector potential: 



--,y^ldx2jii'-r2 lini Ji2' ,2i'{<^o) ■ 
a ^ — 'J ujo^o 



1 



crDC= lim — 

CTi 



dX23\vj21' Jl2',21'(wo) 



m 



(B.22a) 
(B.22b) 

(B.23) 

Sometimes it is covenient to average the coordinate ri over the volume Vol, in which 
case one should replace J^ai i'^ (|B-22p by / 



where we have introduced the retarded correlator [with 612 = ^(^12) 

Jl2',21'(^0) = / dtl2e*"°*i^ 012 C[ir^22'] ■ 



B.2. Keldysh Approach with Source Fields 

We now use the Keldysh real-time approach to rewrite C[iy^22'] in terms of correla- 
tors whose dynamical and statistical properties are governed entirely by Hq: First, 
thermal weighting is done in the infinite past using the free density matrix po, and 
then the interaction is turned on adiabatically. For arbitrary operators Ah and Bh , 
this amounts to the replacement 



{AH{ti)BH(t2)) {Aniti - to)BH{t2 - to))o , 



(B.24) 



where the initial time to is sent to —00 so that all disturbances associated with 
switching on the interactions have decayed in the infinite past (the limit to —00, 
will be understood but not displayed below). Second, the time evolution of all 
operators is expressed in the interaction representation, using the familiar operator 
identity 



AiiU) 



Te 



JHoit,~to)/h 



Ase 



-iHo(ti-ta)/h 



(B.25a) 
(B.25b) 
(B.25c) 



where T is the time-ordering operator (the anti-time-ordering operator, needed for 
tj\, will be denoted by T). The correlator C[ii'.22'] then becomes 

C[iv,22'] = ^{Uj{h,to)hn'i{ti)Ui{ti,t2)n22'i{t2)Ui{t2,to))o (B.26) 



-^{Uj{t2,to)h22'i{t2)Uj{tut2)nn'i{ti)Ui{h,to))o 



This expression can be recovered via functional derivates from the following con- 
struction: 



l',22'' 



^Pll'(^l,^o) 



Sij2'2 



(B.27a) 



v=0 
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Uj^{h,to) = Te-i ■'■'0 MHu(t,)+vj(H)] ^ 27c) 

Vi{h) = j dX3,3 ^33(^3) "33/(^l) ■ (B.27d) 

The index a = F, B will be used to distinguish propagation associated with Uj or 
J7| in Eq. (|B.26[) . i.e. with the forward or backward parts of the Keldysh contour, 
respectively. Since uj^Uip = 1, the denominator in Eq. (jB.27bp . included for later 
convenience, in fact equals unity, p^^, (ti, to) — p{ti, tQ;xi, xy) is the reduced single- 
particle density matrix. We call it "reduced", since the thermal average ( )o in 
Eq. (jB.27bp traces out all electron degrees of freedom but one, to be called the 
"singlcd-out electron", for which the others constitute an effective environment. 
Note that we have defined p^^^^, (ti, to) in the presence of a source terrrlB to generalize 
this to vits), defined by Eq. ()B.27dp [which uses the conventions of Eq. (|B.6|) ] 
on the interval ^3 G [toj^i] in terms of a real c-number "source field" £'33(^3) = 
v(t3;x-^,X3) that couples to the (not normal-ordered) operator ^357(^3). The source 
field is devoid of physical meaning, and is introduced merely as a mathematical 
device to generate C[ii/^22'] via functional differentiation. For v = 0, our reduced 
density matrix Pi2{t, to) corresponds to p{t; ri, r2) of (GZ-II.20) of GZ, who simply 
call it "density matrix" . 

In the usual Keldysh approach, all time integrals involving the interaction extend 
from —00 to 00. This can also be achieved in the present approach, if desired, by 
inserting a factor of 1 = Ulg{too,ti)Ujp{too,ti) just to the left or right of nii'i{ti) 
in the first or second lines of Eq. (jB.26p . respectively, and taking the limit too — > 00, 
^0 — > —00. However, the actual value chosen for t^o does not matter, and in the 
present approach, it is actually simplest to use too = ti. 

B.3. Hubbard- Stratonovich Transformation 

Following GZ, we now decouple the interaction term Hi in Uja using a Hubbard- 
Stratonovich transformation that introduces a path-integral over a further pair of 
real, spin-independent c-number fields Vait^^r^): 

U it , ^ /^^afe,r3)^.(ti,to)e-^r(*^-*»)A 

''''^^'''''^ - lVV^{t,^r,)e^^l^^^^Min ' ^^'^^^^ 

S^y^t.M) = /t^Ss^^^^^^ , (B.28b) 



to 



(27r)3 2F'"*(g) 



''For our purposes it turns out to be sufficient to use the same source term v and source field C33 
on the forward and backward contour; to calculate correlators more general than p^^, one would 
introduce a separate source term i)" and corresponding source fields Cgg for each of the forward 
and backward contours, a = F/B. The corresponding generalizations below are straightforward. 
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Ua{ti,to) - Te-i^'o 'it3[K(t3)+«Kt3)] ^ (B 28c) 

Va{t3) = / dx3eVait3,r3) -.nssi: (ts) • (B.28d) 



The fields Va{t3, r^) and their Fourier transforms 14(^3, q) are defined on the interval 
^3 S [^Oi^i] on the upper or lower Keldysh contour for a = F or B, respectively 
(i.e. the time argument of Va is understood to carry an implicit index a). By using 
Eqs. (jB.28a|) to rewrite all Uia in Eq. (|B.27bp in terms oiUa, the reduced density 
matrix can be expressed as 



Pi,,{ti,to) = {p^,,{ti,to))y , (B.29a) 

Pii'(ti,tQ) = y/, , X , (B.29bj 

Ziti,to) = {ul{ti,ta)UF{tiM))a , (B.29c) 

iS^^\ti,to) = i{S^^ - S^^) + hlnZ . (B.29e) 

In Eq. (IB.29ap . the reduced density matrix p^y is expressed as a functional aver- 
age, over all configurations of the fields Va, of the functional pj^y{ti, to). The latter, 
defined in Eq. (jB.29bp (and called pv by GZ), is the reduced density matrix cor- 
responding to a particular configuration of the fields Va- For any such functional 
JF[V"a], the functional average is defined by the functional integral (IB.29dp . with 
an effective action Sy^ given by (IB.29ep . Note that Sy^, via its dependence on Z, 
depends on the source field v. 



B.4. Roadmap for GZ's Strategy 

If, in Eq. (|B.29bp for p-^-^,(ii,to), the evolution operators are expanded in powers 
of the Va's, the standard Keldysh perturbation expansion for t hese correlators would 
result (as recapitulated in App. IE|). The approach of AAgP^ amounts to doing just 
such an expansion to order V^- However, such a perturbation expansion has infrared 
divergencies which are cured only when the leading divergencies are summed to all 
orders (or by introducing an infrared cut-off by hand, such as an external magnetic 
field, as done by AAG). At present, no exa. ct way of summing the entire perturbation 
series is known. Already in 1982, AAk[I^ were able to perform a summation of the 
leading infrared divergencies by treating Va as a classical field; this indeed cured 
the infrared problems, but neglects the quantum nature of Va, hence corrections to 
AAK's calculation are to be expected at sufficiently low temperatures. 

GZ attempted to proceed both beyond AAK's calculation (by including quan- 
tum corrections) and beyond perturbation theory (by summing an infinite subset of 
the perturbation series). The essence of their idea was to integrate out all electron 
degrees of freedom but one, the "singled-out electron", thereby deriving an influ- 
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ence functional describing the effect of the other electrons (an effective dissipative 
environment) on the diffusive motion of the singled-out electron. To this end, they 
adopted the following strategy, which we shall repeat below in our own notation: 

(A) An exact equation of motion is derived for to) [(GZ-II.24), our (jB.33p ]. 

(B) From this, another exact equation of motion is derived for the linear response 
Spii'{ti, to) to the source field v [(GZ-II.39), our (jB.35p ]. together with the form of 
the effective Hamihonian Hfj [(GZ-II.40), our (jB.36P ] which governs the dynamics 

of (5pii'(ii,io)- 

(C) This second equation of motion is integrated exactly [(GZ-II.41), our (IB.38|) ] 
in terms of effective evolution matrix functions U^{t, t') [(GZ-II.42), our (IB.40|) ]. 

(D) A functional derivative of J/Sj^j^, to) w.r.t. to the source field v is taken to 
obtain an expression for the conductivity [(GZ-II.49), our ((B33)) or ((R55)) ]. which 
involves a functional average of the form {U^ p^U^)v over the fields Va [Eqs. (|B.52p ]. 
where = pVa=o is the (initial) density matrix in the absence of interactions. The 
purpose of the subsequent steps (E) to (G) is to facilitate the evaluation of this 
functional average. 

(E) The evolution fmictions introduced in (C) are represented as path in- 
tegrals over the degrees of freedom of a single electron, whose Hamiltonian de- 
pends on the fields Va- We shall use a coordinate-space-only path integral / 'D'{K) 
[Eq. (jB.56p ]. thereby deviating somewhat from GZ, who use position-momentum 
space integrals JVE^VP" [(GZ-IL44), our Eq. dPTall ]. The relation between GZ's 
position-momentum and our coordinates-only path integrals is explained in time- 
slicing detail in App. |D] 

(F) The action 5y* (more specifically, the term \nZ) that governs the weights of 
different configurations of the fields Va in the functional average {U^ p°lJ^)v, is 
expanded to second order in Va, corresponding to the standard RPA approximation 
[(GZ-II.30), our Eqs. (fRBl) and (IrM)) ]. 

(G) The density matrix Pjj/ (ti, to): whereever it still occurs, is approximated by its 
noninteracting {Va ~ 0) version p^^, . [GZ make this approximation twice: (i) in the 
propagators U^, to obtain (GZ-II.43), and (ii) in the initial-time thermal averaging, 
to obtain (GZ-II.49); we use the analogue of (i) [Sec. IB.S.T) . but do not need (ii).] 

(H) The functional average ( )y, which through the approximations (F) and (G) 
has been reduced to a Gaussian functional integral, is performed to yield an effective 
action iSR + Si [(GZ-IL54), (GZ-II.55), our (lR79)l . ([BJS])]. This effective action 
determines the influence functional of the environment (the other electrons) on the 
singled-out electron. 

In GZ's paper the above steps are presented in a somewhat different order: ap- 
proximation (F) is discussed already after (B), and approximation (G) is made 
directly after (C). We prefer to carry out the steps in the order stipulated above, 
because this allows us to postpone each approximation to the latest possible stage. 

The results derived by the above steps are used in Sees. IB.6I and IC.3I to make 
contact with diagrammatic perturbation theory, and in the main text [Sec. g to 
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calculate the decoherence time. For the latter, we continue to follow GZ's approach 
in spirit, but use a more careful treatment of a certain "Pauli factor"; remarkably, 
this turns out to lead to AAK's result for instead of GZ's. Although the details 
of this calculation are presented in the main text, we shall now summarize them 
here, too, in order that the present brief overview of GZ's strategy be complete. 
(I) The term iSr in the effective action turns out to depend on a certain "Pauli 
factor" (S — 2ff), which we treat differently from GZ: In their position- momentum 
path integral it is represented as [1 — 2no{hQ)], where hQ{R{t), P{t)) is the single- 
particle energy of the singled-out electron, which GZ assume to remain constant 
during the diffusive motion. In our opinion, this assumption neglects recoil effects 
associated with electron-electron interactions [see Sec. IB. 6.2] . Therefore, we instead 
use a Fourier representation of the effective action, in which the Pauli factor is 
represented as [1 — 2nQ{h{e — Cj))] [Eq. (jB.91|l ]. where he corresponds to GZ's Hq, 
and u) is the frequency transfer upon emission or absorption of a photon. 
(J) The path integrals J T>R' for the singled-out electron are performed in the 
saddle point approximation, meaning that only the contributions of pairs of time- 
reversed diffusive (or "random walk") paths are retained. 

(K) The average of the influence functional over all such random walk paths, namely 
(^^~{iSr+Si)^^^^ jg approximated by the exponentiating the average of the effective 
action, e-<*^«+^^>™ [(GZ-IL67), our Seed]. 

(M) The exponent F{t) — {iSb. + <S'/)rw, a growing function of time, is evaluated 
by Fourier transforming the effective action into the frequency-momentum domain 
and averaging the Fourier exponents, using (e*'? [-'*(*3^^-'^(*'*^l)i.w — e^^ ■Ci|t34| j^^j. 
Eq. dSl)]. 

(L) The resulting function F{t) is used to identify the decoherence time as the 
time for which F{t^) becomes of order unity [(GZ-II.67), (GZ-II.70), or (GZ-IIL6) 
(GZ-III.22), or our Eq. 

B.5. Repeating GZ's Strategy in Detail 

The remainder of this appendix is devoted to a detailed discussion of steps (A) to 
(I), using our own notation. 

B.5.1. Exact Equation of Motion for Piii{t,t()) 

To derive an exact equation of motion for Pjj/(i,io)j we start from the simple rela- 
tions 

ihdtipjit^x) — ho{x)i()j{t,x) , (B.30a) 
ihdtKit) = [Va{t) + viit)] USt) ■ (B.30b) 

Since all functions in Eqs. (|B.30[) are evaluated at the same time i, as are all other 
functions needed below up to Eq. (|B.37|) , we shall suppress the time argument below 
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and use the shorthand notation 

From Eqs. (jB.30[) . we then readily find 

7-1 



ihdtZ 



ihdt{U^g n-j, Uf)o - Pu'idtZ 



Eq. ()B.32a|l can be brought into the form 



I-Oi' 



Hi' 



dxk 



P^ki<^Vl 



Fk 



eV, 



BkfPki' ^ Pik^ki' ^ikPki' 



by using the identities 



il4'-4^k 



\k^k^' 



Z \lll^l^li:kAllF)o = Pu'Pkk - hkhr 



(B.31a) 
(B.31b) 

(B.32a) 
(B.32b) 



(B.32c) 



(B.33) 



(B.34) 



The last of these can be checked by expanding both sides in powers of V^, and 
evaluating each term in the expansion using Wick's theorem. Since Va is quadratic 
in ■0's, one readily finds that the combinatorial factor for each topologically distinct 
diagram is just equal to 1, and that the left- and right-hand sides of Eq. (jB.34[l 
generate precisely the same set of topologically distinct diagrams. 

Eq. (jB.33P is the desired equation of motion for p-^, . [It reduces to (GZ-II.24) 
upon setting the source fields to zero, v — Q and recalling our footnote m] The 
term on the right-hand side of Eq. (|B.33P that contains a term quadratic in p, 
coupling to e{VF — Vb), will be seen below to be responsible for enforcing the Pauli 
principle. Note that Eq. (|B.33P contains only c- number functions (no hats occur). 
Hence the order of factors in products does not matter as long as their subscripts 
are displayed explicitly (the derivatives contained in the functional operator h^ii 
should be understood to act on index i' of p^^, even if we write them in the order 
Pii'hoi'). Nevertheless, the subscripts do imply that the products have the structure 
of matrix multiplication in coordinate space; we hence chose to write the factors in 
an order that is suggestive of this matrix multiplication. [This order conforms to 
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that used by GZ, in whose notation the coordinate indices are not displayed, but 
are implicit.] 

B.5.2. Equation of Motion for Spu' (i, to) 

Next, we expand the reduced density matrix to linear order in the source field 
(which is sufficient for a linear respone calculation of the conductivity) by writing 
Pii' — PiT'' ~^^Pii> ' where the superscript (ns) denotes "no sources" and Sp-^, is linear 
in V. It satisfies the following equation of motion, found by expanding Eq. (jB.33|) . 

ihdtSp^.^, ^ Dw + JdxiH:^6p^^, -Jdxr Sp^^.H^i,, (B.35a) 

A^' = j dx^v.^jp[p - J dxv , (B.35b) 

where the effective Hamiltonians and are defined as follows 

= ho, 4 + , Hi = 4 ho, + h^,, , (B.36a) 

(B.36b) 

a=± 

w^^eSa, wtr^SahiSu-^p^^). (B.36C) 



[Eqs. (|B.36p correspond to (GZ-II.39,40); their —eVx{ri) corresponds to our vu.] In 
Eqs. (|B.35[) . the combination of indices a or iV, one without bar, one with, will 
always refer to two independent position indices associated with the same time (i.e. 
ti = ti). The Hamiltonians and are associated with propagation along the 
upper and lower Keldysh contours, which is why in Eq. (|B.35ap they are contracted 
from the left or right with the first or second index of Sp^^,. In Eq. (jB.36c[l for 
the vertices w^" and elsewhere below, the symbol Sa stands for "sign of a", with 
sp/B = ±1- The fields Vai = Vai{ti,ri) (with a = ±) occuring in Eq. (jB.36bp are 
defined as symmetric and antisymmetric linear combinations of the fields Vai (i.e. 
the time and coordinate arguments of V+j and V-i on the upper and lower Keldysh 
contours are both equal to (ij, rj)]: 



V+A ^ /1/2 1/2 \ fVp, 



(B.37) 



Note that H^, like hoii is a c- number functional operator - the derivatives contained in hiuS^^ 
get "transferred" onto the function it multiplies: 
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Since both and iJ^ depend, through Vm, on both Vpj and Vbz, crossterms 
will occur below that link the forward and backward Keldysh contours. Note that 
the field Vm, which shall always carry a "barred" index below, is contracted with 
the second index of w^" in or the first index of w?^" in h^^, respectively. Thus 
V- and 'w°'~ "do not commute" , which will be important below. The factor {6 — 2p) 
in w°'~ will be seen below to account for the Pauli principle. 

All functions occuring in Eqs. (jB.36|l depend on the same time argument t, which 
we henceforth display explicitly again. It is worth emphasizing that, through their 
dependence on p[j^\t,to), the expressions Sp-^,, 1)^^, Hf^ and C/| [defined below in 
Eqs. (jB.40|l ] all explicitly depend on the initial time to, too, although, for notational 
brevity, this ^o-dependence will be displayed below only for pff{t, to). 



B.5.3. Exact Expression Jar 5pu' (t, to) 

The formally exact solution of Eqs. (IB.35P can be written in the form 
W - fdt'fdx.rk Uf,it,t')D,-,it')U-^^,{t\t), 



(B.38) 



where the functions {t,t') and Ujl{t',t) are defined by the requirements that 



M{t,t)^ufs,t) = s^^ , 



ihdtM{t,t') 



ihdtUfS'.t) 



dx,H^,{t)U(At,t') 



dx-M {t',t)H-f^{t) 



Eqs. (|B.39[) are fulfilled by time-ordered exponentials 



1 



(5^. - - I dt^Hf^its) - I dt3 I dti I dxkHf^{t:i)H^^{ti) 



= 5,, + - 



' Jt, 



(B.39a) 
(B.39b) 

(B.39c) 
(B.40a) 

(B.40b) 



di3 Hf^ih) 



1 



dti / dta / dxk Hfk{t:i)H^Si) + 



where we always take t > t' , and where each "internal" product of two factors 
^ik^kj that arises when expanding the exponential involves a further coordinate 
integral J dx^ ■ [Below, we shall often suppress time arguments and use the short- 
hand Ufj = U^{ti,tj) and likewise f/| = U^{ti,tj).] Note that the time-ordered 
exponentials (|B.40aP and ()B.40bp for and U^^ are defined in terms of the power 
series expansions indicated above; the same is true for all path integral representa- 
tions of and U-^ to be used below. Note also that U°j is spin-diagonal, since this 
is the case for Hfj{t) = 6a,a,H°-{t, ri,r^) and pa{t, to) = 8^.„.^p{t, to; n, n). 
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Eq. (|B.38P corresponds to GZ's exact Eq. (GZ-II.41). Note that the procedure 
by which it was obtained, namely, first differentiating p and then integrating dp 
w.r.t. t, has produced a result in which the reduced density matrix p^'^^^ appears in 
the exponent, via its occurence in and U^. Accordingly, the effective action to 
be derived below will likewise depend on . 

Let us now also derive an equation of motion for the time evolution of the density 
matrix in the absence of source fields, p[ji''\ since we need it in Eq. (jB.SSp . where 
it enters via the D of Eq. (|B.35bp . [This point is not discussed by GZ, who simply 
replace p''^,^^ in Eq. (jB.38P by p^,, as discussed in Section [B.5.71 ] Evidently, the 

desired equation of motion for p'^j,^'' is the v — version of that of py', namely 
Eq. (jB.BSp .-i-n. which can, in analogy to Eq. (jB.35a|) (without its first term), be 
rewritten as 

^n^,fP = jdx-M'iPP -Jdxr f^H^, . (B.41) 

Here the primed Hamiltonians H'^^ are defined by equations identical to Eqs. (|B.36p 
for the unprimed ones, except that the vertices wf^ of Eq. (|B.36cP are replaced by 
primed vertices w'^'^ that are defined as foUowsEI 

w'^+ ^e~5u. w[r ^ sa\e{5.^ - y^2p^'^) . (B.42) 

Here the y^^^ £ [0, 1] are (arbitrary) real numbers, with + = 1. It will 
turn out below to be convenient to let the choice of values for depend on which 
contour the current vertex at time t2^ is located: if it is on contour F/B, we shall 
choose y^^^ =0=1 — y^^^ (compactly: = for j(^2a) on contour a; Fig. IBll 
below shows an example with d = F). The solution of Eq. (|B.41|) can be expressed 
as 

pP{t,to) ^U'i{tMplkU§{to,t) , with yP + y^^l, (B.43a) 



K (tM _ ^^/l-,U-^,ito,t) if = 0, = 1 



uKit,to)pl-,u;^,ito,t) 



(B.43b) 



if = 0, = 1 



The primed propagators [7|° are defined analogously to Eqs. (|B.40p . but with H 
H' everywhere. In Eqs. (jB.43[) we have implemented the standard initial condition 



'^The reason for the extra j;" in front of p]'^^ for w'^ , which is the only difference compared 
to w^-~ of Eq. l lB.36cl l. is as follows: The linear response equation of motion for Sp -, contains 

two different contributions that are quadratic in p, namely p'tJ'''6^_fci5p^^, and (5p^j,ey_j.p^"^' , 
which in Eq. I IB.35al l were grouped with the first and second terms of respectively. In contrast, 
the equation of motion for p^'^i'^ turns out to contain on the right-hand side just one type of term 
quadratic in p, namely p^5J''^eV_j,p^'Jf', with total weight 1. By using + = 1 in Eq. I IB.421 1. 
we have distributed this term with weights y^ and y^ among the two terms on the right-hand 
side of Eq. llRiTl l. 
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for the Keldysh approach, namely that at time t — t^, the density matrix was free, 
i.e. p-^f^ {to, to) = p^, . Below [cf. Eq. (|B.51bjl ]. we shall insert Eqs. (|B.43bp . with to ^ 
— oo , into Eq. (|B.38|) , where it enters via the D oi Eq. (jB.35bp , to ensure that thermal 
averaging is done in the infinite past. This is an important improvement relative to 
an approximation used by GZ, who simply replace p'^y,^\t, to) in Eq. (jB.SSP by p^, ; 
they thereby effectively perform thermal averaging with a nonequilibrium initial 
density matrix, as discussed in Section [B. 5. 71 

The way in which t/J^, and p\^^'' differ from their free versions is evidently 
through their dependence on the fields Va and the density matrix pij in Eqs. (IB.36j) . 
Let us now briefly discuss their free versions. First, in the absence of all interactions 
the expression for the reduced density matrix p'^^^^ reduces to the form 



p° = ^^(^^) noi^x) = ip%r , (B.44a) 

A 



where no{(,x) = [e^^^"^ + 1] is the Fermi function, and ipx{xi) are the exact single- 
particle eigenfunctions of hoi, with eigenvalues [cf. Eq. (|B.15|1 ]. Next, let U^'^ 
denote the propagator to which reduces in the absence of interactions, i.e. for 
Vai = in Eqs. (|B.36P [so that H°j = hoiSij]. Its explicit form is easily found by 
constructing an object satisfying the defining Eqs. (jB.39|l for Vai — 0; the result is 
independent of whether a = F or B, and given by: 

Ul" ^ Ul = ^^(^^) e-'^^*'^/" = tHG^ - Gf) , (B.45) 

A 

where G^^^ = ±0{±tij){G^ — G^) are the standard free retarded and advanced 
electron Green's functions, with 

[{^pJ{t„x^)^p}{tj,Xj)}o) X 

It follows that for a given time order, as occurs under a time-ordered integral, U^"" 
is equal to either a retarded or an advanced Green's function; e.g., for ti > tj, we 
have U^^ = ihG^ and Ujf = —ihG^^. Nevertheless, it will be useful to generally 
retain both terms in Eq. (jB.45|) . because that allows expressions involving the free 
reduced density matrix to be simplified by Fourier transforming from the time to 
the frequency domain: for example, denoting the frequency Fourier transform of 
Uij{t) by C/j^(w), we immediately find the exceedingly useful relations: 



dx,f^,U-^^iLo) - noihio) Uliio) - -ihG<iLo) , (B.47a) 
dx, [S - 2p% m[Lo) = [1 - 2no(fi^/2T)] U^Alo) = ihG^itu), (B.47b) 



to the frequency represantion, the Pauli factor {6 — 2fP) in Eq. (jB.47bp gets mapped 



where G^ — G^ + G^ is the Keldysh function. Note, in particular, that by passing 
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onto [1 - 2no{hLj/2T)] = ta,nh{hij/2T), a fact that will be very useful in Sec. IB. 6. 21 
below. 

For future reference, we note that when the matrix propagators and 
[with t > t'] are expanded in powers of hy [i.e. in powers of the fields Va, see 
Eqs. (|B.36bp ]. they take the fo rm of time-ordered or anti-time-ordered power series, 
respectively: 




I I 'J '-'il ■ ■ ■ "'VNn'^ Nj 

\(,A\NfTOBlB frOBlB frOB 

[[+^) ^ jN'^VNN ■■ -^21 '^Vll^li 



dxi idx2_2 ■ ■ ■ dx]\[ ]\[ 

(B.48) 



These expansions [illustrated in App. D by the third row of Fig. lDlj are alternative 
but equivalent to those of Eqs. (jB.40p . and, just as the latter, can be regarded as 
formal definitions of U^, and of all path-integral representations thereof to be used 
below. Note that for each occurrence of a "vertex" hfl, t or h^,^ , , the vertex 
coordinates xi^ and xj^ are both associated with the same time U, and both are 
integrated over in / dxi j [cf. Eq. (|B.6|) ]. This need for a double position integral 
at each vertex is a direct consequence of the fact that the effective Hamiltonians 
of Eqs. (|B.36p are nonlocal in space. Since the integrals in Eq. (|B.48p are time- 
ordered, each [/"^ occuring in can be replaced by ihG^, and each U'^^ in 
by —ihG^ [see Eq. ([D.llcP ]. Indeed, the latter replacements are, in effect, used in 
the path integral representation of i7° to be introduced below [Eq. (|B.56p ] . We have 
nevertheless chosen to write Eq. (jB.48P in terms of ?7°° functions, as a reminder 
that the density matrices occuring in the interaction vertices hy can be converted 
to Fermi functions using Eqs. (jR4Z|. 

B.5.4- Exact Expression for Conductivity 

The density-density commutator C[iy,22'] [Eqs. ((RTQ)) . ((R27)) . (IB.29ap and (|R38)) ] 
can now be obtained by taking the functional derivative of Spi2{t) with respect to 
the source field v [occuring in ([R38l) via D of (|B.35bp ]. Henceforth writing t = ti 
and t' = t2, the result can be written as 



C[ll',22'](^l - ^2) — ^ >^12''',21'(^l'*2;io) + C[^r*22V (B.49) 

a.=F,B 
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where C^^!''22''] ^ contribution irrelevant for weak localization, which will be 
dropped henceforth^ The J^^^'s are defined in terms of the correlator 



"^12',22',21'(*l' ^o) = -^Uu' {tl,'t2)p''^,\t2,t())ij2ii (t2,ti) (B.51a) 
^]:J dXoroU[2'itl,t2)U!2^{t2,to)ploUci{to,t2)Uiy{t2,h) , (B.51b) 

[the second line follows via Eq. (jB.43a[) ] by: 

'^i2'^2i'(^i'^2;io) = /'^2;2 (ji^, 22 2i'(*i'^2;io)) (B.52a) 

' J \ ' ' ' V.ns 

1 / ~ p ~ IP ~ g \V^=0 

dXofi ^ \^t^l2' U20 Poo '^01' )y^^ ' 

j[^,\^,iti,t2;to) = - Idx2 bX2 22' 21'itl,t2;to)) (B.52b) 
' J ^ ' ' V,ns 

^2;o,0^ yC^lO Poo ^02' ^21'/^^^_^ • 

J^^^ [illustrated in Fig. IBlj and J*^^^ denote correlators that have a current vertex 
inserted on the forward or backward Keldysh contours, respectively. As a notational 
reminder, the indices 2, 2', and 2 here all refer to the same time, t2 in this case, 
and after performing the derivatives implicit in jui and h'^} , we have to set 2 = 2'. 
However, X2 in Eqs. (jB.52p is an independent integration variable. The subscript ns 
(for "no sources") in ( )v,ns indicates that, following the prescription of Eq. (jB.27ap . 
all remaining v-dependencies are to be dropped henceforth by setting v — 0. The 
second lines of Eqs. (|B.52a[) and (jB.52bp . in which we set = for the corre- 
lator J^°''> containing the current vertex on contour a, follow from Eq. (jB.Slbp for 
^^12' 22' 21' using the first or second line of Eq. (jB.43bp for U'^ p° U'^ , respectively 



The term Cm"'22°1 '"^ 



^hV%2'] = «(Pl2(*l'*0)7r J--^)v,nB - i(Pi2{*l:*0)>V,ns (7^ J-T}v,ns, 

^^^^ ' ^ {Ul(tlM)Up(H,t2)h22>l(t2)Up{t2,tQ))o 



{Ulit2,to)h22'l{t2)Uliti,t2)Up{t2,to))o 



and arises since the effective action * of Eq. l|B.29e| l in the functional average p.29d[l depends, 
via InZ, on v too. CjJjJ'*22/'] corresponds to (GZ-II.47) and is neglected by GZ [see the discussion 
after (GZ-II.47)], because in the absence of interactions, it vanishes entirely, and hence does not 
contribute to the weak localization correction to the conductivity (in other words, C^JJ^'22°] 
irrelevant to the question how this correction is affected by interactions). We shall not consider it 
further either, since in diagrammatic terms it corresponds to Hartree contributions to the electron 
Green's fun ctions, which merely renormalize the magnitude of the conductivity (and were neglected 
by AAdm, too). 
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[thereby conveniently avoiding primed propagators U' under the a;2-integrals on the 
"other" contour a' ^ a, which thus have the form J dx2U^2 ~ ^ij ' latter 
composition rule follows from Eq. (jB.45|) and the completeness of the waveftmctions 
ip\{xi) occuring therein.] 




Fig. Bl. A pair of backward (B) and forward (F) paths contributing to jjj/ (*i > *2 ; to)i with 
ti > t2 > to- There are two ways to view this figure: (i) Ignore the wavy interaction Unes, double 
dot vertices and frequency assignments; then this figure illustrates the second line of Eq. IIB.52al l. 
and the solid or dashed lines represent the full unprimed [Uy ) or primed {U'y ) propagators, 
respectively, (ii) Imagine the propagators U^/^ and tj'^ to have been expanded in powers of the 
interaction [as in Eq. | |B.48| |]. This generates a forward and backward backbone of free propagators 
W '^^ j"^ (represented by either solid or dashed lines, which now have identical meanings), which 
are respectively decorated by the vertices hy^:^ and h^-^ (represented by a pair of dots; both dots 
are associated with the same time, but the one drawn on the side of earlier times is distinguished 
by a bar; the origin of this convention is explained in App. |D] Fie. IDll l. The vertices generate, 
after averaging over the fields Vai, the wavy interaction lines Ciaj , , connected to the barred 
dots. [The interaction lines are labelled according to Eq. (IB.84t below: C^a' stands for C^^^j , , 

, or , if generated by (V+i^V+j ,)v, {V+'ia^-J i)v, or {V-ia^+J respectively, 

cf. Eq. l lB.75bl l.] For both cases (i) and (ii), arrows are drawn to point from the second index 
to the first index of each of U^j, and C^Jj^f^. Thus, they point from later to earlier times 
along the backward Keldysh contour, and from earlier to later times along the forward Keldysh 
contour (i.e. they form a continous loop, starting on the backward contour from ti backwards 
to to = — oo, then continuing on the forward contour from to = — oo forwards via t2 to ti). 
Finally, the frequencies label the interaction correlators Caa'i'^) ^md Green's functions G^^^{uj) 
and G'^iuj) = ta.nh{hjj /2T)[G'^ - G^]{uj) that arise (before disorder averaging) upon Fourier 
transforming from the time to the frequency domain, as for Eqs. I IB.86I 1 or I IC.16t below. The 
effective action defined in Eqs. JJ} to H4a| l of the main text neglects the frequency transfers uii in the 
arguments of all retarded and advanced electron Green's functions [G^^^{e—LUi — ...) G^/'*(e)], 
but, for every ^(u}i)G^ {e — uii), retains it in the factor tanh[?i(£ — uji)/fi\ of the accompanying 
G^ function. [As discussed in Sec. [3l or IB. 6. 21 this is justified by the fact that all integrals over 
frequency transfer variables are limited by Fermi factors to the range | tkUi | ^ T] . 



Inserting Eq. (|R49)l for C[ii/.22'] into Eq. (|B.22a|) . the expression for ctdc that 
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results upon representing the applied field in terms of a scalar potential, and then 
relabelling X2 <^ 2^2 the term containing J^^\ we find: 



CDC 



2m d 



J2f'dt2f dx2-2 (B.53) 



[rz-rg] • (Vi - ViO(Ji^,,22,21'(*l'*2;io))y,„s ■ 

Eq. (jB.53p for the DC conductivity is analogous (but, as discussed below, not 
identical) to (GZ-II.49) [the factor Ul^,{ti,t2)Ug,{t2,ti) which occurs m our 
'^i2',22,2i'(*i; ^2; ^o) IS the analogue of the function J{ti,t2;ri,r[;r2,r2) occuring 
in Eqs. (GZ-II.49) and (GZ-II.50)]. In deriving Eq. (|B.53[) . no approximations have 
been made, apart from not displaying the Hartree terms [cf. footnote |ej. 

Instead of (jB.22a|) and (|B.53[) , it will be more convenient for our purposes to use 
Eq. (jB.22bp as alternative expression for ctdci derived by representing the applied 
external field via a vector potential. The correlator Ji2',21'('j-'o) occuring therein can 
[via Eq. (jB.49[) ] be expressed as: 

/C30 
dii2 e*'^«*i^ lim V euJi^l 2i'iti,t2;to) ■ (B.54) 

Since Ji2',2i' stems from the commutator C\iy ^22'] [Eq. (|B.19p ]. whose terms satisfy 



y'-lV H^22' HI ~ y'-22' H^lV HI 



the correlators J^"' satisfy 



Jl2'',2vi^"L,t2;h) — — ^j/^|2'l(*l' *2; ^o) , >^1^^21'('^o) — ~ >^1^2^2'1 (""^o) 

The first of these [which implies the second] is manifestly obeyed by Eqs. (|B.52[) . 
Taylor-expanding Eq. (|B.22bp using J{loq) = J(0) -I- ujqJ'{Q) + . . and separating 
CDC — CDC. real + icTDc.imag into real and imaginary parts, we obtain 



CTDCroal J dX2 ju' ■ j 22' J[2',21'i^) ' (B.55a) 

1 \ie'^{niiH) 1 /" , • . 7 

«0-DC,imag = hm } h 3 /rfa;2 Jll'-j22' Jl2'.21'(0) 

0)0^0 ujo m a J 



Since we have taken the DC-limit loq — > 0, the imaginary part CTDCimag must be 
strictly equal to (to all orders in the interaction), which is a useful consistency 
check. 

In App. [C] we show how Eqs. (|B.55p can be massaged into more familiar ex- 
pressions for (Tdc, both in the absence and presence of interactions [cf. Eqs. (jC.Sp . 
dCJi), dCM ]. 

B.5.5. Coordinate- Space Path-Integral Representation for 

In this subsection we shall derive path integral expressions for the objects in terms 
of which the conductivity is expressed in Eq. (|B.55ap , namely the propagators C/| 
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[Eq. (jB.40p ] and the interaction-averaged correlators (J^°^)y [Eqs. (|B.52|) ]. We devi- 
ate from GZ's approach, who used a path integral / T>R J P over both coordinate 
and momentum space, in that we shall use coordinate-space-only path integrals 
/ T>'{R), because that makes possible a more accurate treatment of the crucial 
nonlocal Pauli factors {6 — 2p) in the effective Hamiltonian H'^ of Eq. (|B.36p . 

We begin from the power series expansions (IB.48|) of the evolution matrix func- 
tions U°j{t,t') of Eqs. ()B.40|) in powers of hy, and introduce, as a shorthand for 
these expansions, the following coordinate space path integrals: 



dt 



V3f3f 



(B.56) 



Here Sa stands for sf/s = ±i and the index value a — F or B should be used for the 
upper or lower term in the curly bracket, respectively. The coordinate-space path 
integral is over all paths ^^{tz) that begin at time t' at point Tj and end at time t 
at point r^; the time that is used to parametrize this path R'^it^) is understood 
to refer to the upper or lower Keldysh contour for a ~ F oi B, respectively [in this 
sense, an index a on ^3 is implicit, as in i?°(t3^)]. The objects 5*0 and hy in the 
exponential factors in Eq. (jB.56|) are both functionals of the path R'^{t3): Sq is the 
standard action for a noninteracting electron in a disorder potential, 

SSit,t')[R''{t3)] = J^' dh - 14mp(fi"(i3))] , (B.57) 

whereas in the second exponential, we used the following shorthand notation: 
^y3.3. = E<\^"3. = 5]^^"[t3,r3^(i3),r3^(i3)]V^a[i3,r-3^(i3)] , (B.58a) 

a— ± Q— ± 

~K3b3b = Y.^^^b'^I:3b = E^"[*3,r-3^(i3)]*^"[i3,r-3^(i3),r-3^M .(B.58b) 

a— ± a— ± 

In App. [D] we give an explicit definition of the path integral Eq. (|B.56p by time- 
slicing the time interval [i', t] [Sec. ID. 4] . and a detailed demonstration that it satisfies 
the defining Eqs. (|B.39p [Sec. ID. 2] . The explicit derivation given there shows that, 
when writing down the path integral (|B.56P , the following points are to be implicitly 
understood [see also Fig. IDll of App. D]: (i) The path integral (|B.56p is simply a 
short-hand for the time-ordered power series expansion (|B.48p . with {—i/K)U^^j^ re- 
placed by Gf^j^ and {i/K)!)^^^^ by Gf^^^ [cf. Eq. (jP.llcp ]. (ii) For each occurrence 
of a "vertex" ^^3^3^ or ft-^g^g^, the vertex coordinates r^{H) and ^^(^3) are both 
associated with the same time ^3 , and both are assumed to be integrated over in the 
path integral [as in Eq. (|B.48P ]. thereby taking into account the nonlocal nature of 
the Hamiltonians hy^^. (iii) The associated integrations / dx^^^ are understood to be 
included in the measure / VR"- (the prime serves as reminder of this fact), in addi- 
tion to the integrations associated with propagators between vertices, (iv) Vertices 
are connected by propagators of the form G^'^j^ or Gj^^^ on the forward or back- 
ward Keldysh contours, respectively. However, since these propagators occur under 
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time-ordered integrals anyway, they can equally well also be written as {—'i/fi)U^pjp 
or (j/^)^j'jfii3 1 is convenient [in order to exploit Eq. (|B.47p ] whenever they are 
contracted with a density matrix p^^j^ or p^^i^- 

Now use the path integral representation (|B.56p (twice) in Eq. (jB.51bp for 



{J^)vas^ and interchange the order of averages, ^ / V^R) . . .J — > J 'D'{R) 
[The latter step could have been postponed until the beginning of Section IB. 5. 
but is used already here, since it simplifies subsequent expressions. Its use, sooner 
or later, is a crucial ingredient in GZ's approach. Its far-reaching consequences are 
discussed in detail in App. IA.4I ] We obtain 



xlf ' if/ V\R) / " i " V\R) T^t^,to)m , (B.59) 

J2'p J2b JOf JOb 

where ij" ^ T>\R) is used as a shorthand for the following forward and backward 
path integral between the specified initial and final coordinates and times: 

V'{R)...= / P'i2^(t3)e<(*-*f)/'' 

X / P'i2^(<3)e-'^o"(*?'*f)/^.. (B.60) 

jRB{tf)=rf 

The influence functional T(^ti.to)[^°'] i^^ Eq. (|B.59[) is defined by the following func- 
tional integral over all configurations of the fields V±3 = V±(i3, r^) of Eqs. (|B.37p . 
with t3 e [to,<i]: 



St?'-8-V (ti,to) 



B-V{h,to)= ['dt3 [ drg ^ 6^(^3,^3)14(^3, r-3) , (B.61b) 
Baih, r,) ^ SFWl%S{r3 ~ rg J + SBWl%S{r3 - rg J , (B.61c) 



• 33 = ^sz '«33 + t'23 



W^^ = 032 W^? + 023 W'^3 ■ (B.61d) 



Here S'^'(ii,io) is given by Eq. (|B.29e|) . and Eq. (|B.61cp . which defines the field 
Ba3 ~ Ba{t3,r3), foUows from using Eqs. (jB.SSp or a primed version thereof, for 
ts > t2 or t3 < respectively. The distinction between the two time orderings, 
which is reflected in the definition (jB.61dp of the vertices (and not noted by 

GZ, since they set — ^2), is necessary, since Eq. (IB.51bl) correspondingly fea- 
tures unprimed or primed propagators Uf^ or U'^" , respectively, which have different 
vertices [compare Eqs. (|B.36cP and (jB.42p ]. Note that Ba3 is itself a functional of 
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both the paths R^^t^) and R^lt^). The influence functional JF(t^ jjj-) [il"] describes 
the eflFect of all other electrons on a pair (forward and backward) of singled-out 
electron trajectories R^{t3) and R^{t3) between the initial time to and final time 
ti. Importantly, this influence functional incorporates the Pauli principle, via the 
presence of the Pauli factor [S — 2p)in w°'^ . 

B.5.6. RPA Approximation 

To evaluate the influence functional ^(ti,to) explicity, our next task is to perform 
the functional integrals J Wa stipulated in Eq. (jB.61aP . As a first (standard) step 
toward making these integrals Gaussian, i.e. doable, we apply (following GZ) the 
RPA approximation: we approximate the effective action Sy^ of Eq. (jB.29eP by the 
part quadratic in the fields V , say 

z4'\ii,to) =z(5°^-5°^) + ftZ(2) = -i [V-A-V]ih,to) (B.62a) 
= -l f 'dt3 f 'dU fdrs dri ^ V^aA^fv^,^ , (B.62b) 

^Jto J to J 

so that Eq. (|B.61aP becomes 



e 2R 



To find we have to find an explicit expression for the term ^iZ^^^ in Eq. (jB.62a[) . 
which arises from expanding the factor \nZ = Z^^^ + 0{y^) to second order in T4, 
using 

Z{t^M) = 1 + Z^^^ + + Oiyl) , (B.64a) 

+ {'^VB{t3)VBiU)) -2{VBit3)VFiU))o} , (B.64b) 

and noting that Z'^^ vanishes, since Va is normal-ordered [cf. Eq. (|B.28dp ]. Express- 
ing Eq. (|B.64bp through the fields (a = ±1) of Eq. ((537| . we find 



= - 



f'hf'duf drs dr4(ll/_3X34V^+4 + Vl3%4l^_4) , (B-65) 

tJ to ^ tQ J 

where Xij (the charge susceptibility) and fjij (characterizing charge fluctuations) are 
defined as 

X^-^-z^0(%)([:n,,Kti):,:n„-/(t,):])o - 46^ ft Im [GfG< ] , (B.66a) 



' The term "influence functional" is used here in precisely the sense in which Feynman used it: 
Our T23 [-R"] is analogous to the quantity F\q(t), q' {t)] of Eq. (12-90) of R. P. Feynman and A. R. 
Hibbs, "Quantum Mechanics and Path Integrals", McGraw-Hill (1965). 
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2h 



{{■.h,u{U):r-n,ji{t,):})o - -e'hRe[G>G<^ 



(B.66b) 



with equal spins, (Ji = aj (for (7^ ^ aj, both these quantities vanish). [Eqs. (|B.66[) 
correspond to (GZ-II.31) and (GZ-II.32).] The right-most equalities were obtained 
by using Wick's theorem to rewrite the correlators in terms of the single-particle 
Green's functions G^'^ [Eqs. (|B.46p ]. The Fourier transforms of and ijij satisfy 
X*-ki^^) ~ Xki^^) — X-fc('^) 8.nd ^_fc(— a;) ~ fjkioj) (thus the latter is real), and 
are related by the fluctuation dissipation theorem [(GZ-II.33)]: 

fjkicj) ^ -^coth{huj/2T) Imxkiuj) ■ (B.67) 

Now, if we write the second-order contribution iSy {t,to) in the form of 
Eq. (jB.62bp , and Fourier transform|3 we obtain from Eqs. (|B.28bp and ([RBS]) : 



i2n) 
I 



olol' 







2if]k{u)) 



(B.68a) 



(B.68b) 



/ act' 



where efc(ijj) = 1 — ^"^'(fc) Xfc(^) is the dielectric susceptibility. [The latter relation 
is a generalized version of (GZ-II.35); in (GZ-II.36), GZ added to efc(w) an electron- 
phonon contribution, which is not important for the present discussion and neglected 
by GZ themselves later on, after (GZ-II.75).] 

Having found A, let us now also find and discuss some useful properties of its 
inverse, A^^ [it will be needed in the next section after evaluating the functional 
integral Eq. (jB.63P ]. Using e_fe(— cj) = e*f.{ijj), we find 



with matrix elements given by 



_/fc(cj) iRk{uj) 
iR-k{-uj) 



Ik{uj) 



efc(t^) 

277fc(cc>)|F-*(fc)|2 



coth(ftw/2T) Imi?fc(w) 



(B.69) 

(B.70a) 
(B.70b) 



IskicoW 

where the last equality in Eq. (jB.70bp follows from Eq. (jB.67p . Note that the as- 
sumptions [stated before Eq. (IB.14|) ] that V^'"'(fc) is real and symmetric, imply that 

R*_k{~Lj) ^ Rk{uj) ^ R^k{uj) , r_k{-cu) ^ Ikiuj) ^ i-k{cj) ^ Ik{-io) , (B.71) 



s strictly speaking the Fourier transform I IB.68al l is an exact representation of V • • V only if the 
time integrals in Eq. l lB.62bl l arc unbounded, e.g. for tg = —oo and ti = oo. In our formalism, this 
indeed is the case, since we do take the limit to — ♦ — oo, and may also take ti — > +oo (because the 
ti -dependence drops out, anyway). 
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SO that the functions Ry and ly are both purely real: Rg — i?*- and lij = /* — 
Iji. For reference purposes, we note also that their frequency Fourier transforms, 
denoted by Ryiuj) and lij^w), satisfy the relations Rij{u!) = Rji{uj) — Rlj{—U!) and 

iy{uj) = i,,{oj) = ^i-Lu) = = - coth{nuj/2T)lm[R,j{uj)] . (B.72) 

Furthermore, Rij{oj) is analytic in the upper half plane, implying that Rij is pro- 
portional to 9(tij). In contrast, lij is symmetric in its indices and thus nonzero for 
both tij > and < 0. 

The components of are of course related to field correlation functions of 
the type {VaiVa'j)v,ns', as follows from a simple exercise in Gaussian integration: 
Introducing the generating functional 

QIC] = {e-^^-^)v,ns , (B.73a) 



[C-V](ii,to) = 



r di3 / drs J2 C"3 ^3 , (B.73b) 



where (^ai = CaiU, fi), with a = ±1, are two real source fields, we find 
QIC] /^^/^^ e-Av-4 v e-K V 

^^^^ /W+/W_e-^^-^-^ ■ ^ ' 

The field correlators are then easily found to have the form 

hi iR'. 



-{yaiya'j)v,ns — — n-. 



n oQaiOQa'j 

i 



^1'4 o) (^-^^^^ 

where Eq. (|R69)) has been used, and the functions [cf. (GZ-II.56) and (GZ-II.57)] 
{R/i)y^ /(dfc)(du;)e-^'^(*'-*')+*('-'-'-^)(i?//)fc(u;) (B.76) 



are defined via their Fourier transforms, given by Eqs. (jB.70a[) and (jB.70bp above. 
Eq. (|B.75bp expresses the general fact [reviewed in App. IE.2| that the field correla- 
tors can also be written in terms of the standard retarded, advanced and Keldysh 
components of the interaction propagator, C^, and , implying that these 
are proportional to Rij, Rji and ly [cf. (GZ-III.A14)]. This implies that Ry is a 
retarded propagator and thereby confirms that it is is proportional to 9{ty ) [as had 
already been concluded above from the analytic properties of Ry{uj)\. 

To obtain explicit expressions for Rk{uj), one needs £k{'^) and hence Xfc(cj), 
for which one has to calculate a polarization bubble [see App. F, Fig. IFlT e)]. If 
= Air/k'^ represents the unscreened Coulomb interaction [Eq. (|B.14[) with 
Ao — 0] and, as is usually the case in the presence of disorder, only small frequencies 
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and wave numbers are of interest, a standard calculation yields [of. Eq. (jF.5e[) and 
(GZ-II.36)]: 



5.5.7. Approximating p^'^ P% 

Even after having made the RPA approximation, the functional integral in 
Eq. (jB.63P over all field configurations of Va is not yet Gaussian. The reason is 
that the term ;8' • V in the exponent depends, via on the full, interacting 

density matrix pfj^\t'), which depends on the fields Va too, in a highly nontriv- 
ial way. To make further progress, we shall ultimately have to neglect the effect 
of interactions on the single-particle density matrix, by replacing p^j^^\t') by its 
noninteracting (and hence time-independent) version : 

Pij (O — ' Pij = {%s)o ■ (B.78) 

GZ use this approximation at two points in their calculation [see the comment after 
(GZ-II.43)]: (GZi) to simplify the propagators U^, namely when passing from (GZ- 
11.40) to (GZ-II.43); and (GZii) to simplify the thermal weighting factor describing 
the initial distribution of electrons, namely to obtain the explicit factor po in (GZ- 
II. 49). In our formalism, (GZii) would correspond to setting to t2, i.e. making the 
replacement p^^\t2,to) p'^^, in Eq. (|B.51ap for J^, 21' (^i: ^2; ^o) and inserting 
the result into Eq. (|B.53p . since this would reproduce (GZ-II.49). 

We shall use similar but weaker approximations, and proceed in two separate 
steps: 

(i) We "linearize" the exponential factor B ■ V in Eq. (|B.63|1 by making the re- 
placement B[p\^^^^] — > ;S[/5|], so that the functional integral (jB.63|l becomes truly 
Gaussian in V and can readily be performed [see Sec. IB.5.8| . We thereby neglect 
the effect of interactions on all occurences (via w"" in hy) of p'"''-' in the propaga- 
tors U^, the rationale being that in order to calculate the decoherence rate, we are 
interested in how the interaction affects the time-evolved propagation of electrons 
along time-reversed paths, and not how it modifies equal-time objects like pij. Dia- 
grammatically, this corresponds to neglecting diagrams which modify the Keldysh 
Green's function without affecting the retarded or advanced ones, i.e. which modify 
only the tanh factor, but not the propagator Uij in Eq. (jB.47bp . 

(ii) For the propagator J^°'\ which is defined as the sum of all terms for which 
the current vertex j22' occurs on contour d at time t2^ , we neglect all interaction 
vertices that occur on the same contour d at earlier times t^- or t^- € [to,t2^]. 
Thus, in the second lines of Eqs. (jB.52a[) and (|B.52b[) . we make the replacements 

U20 and [/g^, — + Uq^,. However, for the opposite contour containing no cur- 
rent vertex, we include interaction vertices for all times G [to,ii]- The rationale for 
this is that, in diagrammatic language, this approximation retains only those dia- 
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grams for which both current vertices J22' and jn' are always sandwiched between a 
G^- and a G'^-function, i.e. G^ j G^. These are the ones relevant for the Cooperon; 
the contributions thereby neglected correspond to the so-called "interaction correc- 
tions" , which feature at least on current vertex sandwidched between two retarded 
or advanced functions, i.e. G^ j G^ or G^j G^. 

Note that this approximation (ii) is much weaker than (GZii): we do not re- 
place p"^} by ill ^l- (!B.51ap (i.e. we do not set t^ ^2), but instead use 
p'^\t2,to) ~ U2oPqqU^, [Eq. (jB.43a[) ] and send ta —00. Also, we wish to em- 
phasize that "interaction correction" terms can be calculated from our formalism if 
one so chooses, by avoiding our approximation (ii) altogether and keeping track of 
all interaction insertions on the entire interval [to,*i] of both contours [Eqs. (|C.10c[) 
and (jC.lOdp give examples of such contributions] . For the sake of greater generality, 
we shall thus for the moment use only approximation (i), and postpone the use of 
(ii) to Sec. EH 

B.5.8. Integrating out the Fields Va to obtain iSn + Si 

The approximations discussed in the previous two subsections render the functional 
integral (|B.63|1 for !F^ti,to)[^'^] Gaussian. In fact, Eq. (jB.63p is just of the form 
(|B.74ap . with ( • V replaced by 6 • V of Eq. (|B.61cp . so that we get 

:F^t^,to)[R''] = e-^^--^"-^ = e-[^s„+s.](ti.to)/n (B.79) 

The exponent {iSu + Si)[R"'] = ■ A^^ ■ B, which is a functional of the paths 
can be regarded as an "effective action" that describes the effect of interactions on 
the "singled-out" electron traveling along the paths i?" . The indices R,i are meant 
to distinguish terms depending on the interaction propagators R and /. Before 
working out the effective action explicitly form, however, let us first collect results 
to obtain path integral expressions for the corr elators j[2'2i' of Eqs. ([R52|) . These 
contain the correlators (J^, 2i')v,(ns)j for which we use Eq. (jB.59|) . with J^{ti,to) 
given by Eq. (|B.79[) . and / dxj integrals, which we perform in the same way as for 
the second equalities of Eqs. (jB.52[) : 



J'l'2/,2i'iM;to) = ±- / dxo,-,, p'o.ou (B-80) 



X < 



ni-F p-^F n^B _ 

F F i V'{K 

J2'p J Of J Ob 
Jof J2r Jop 



y ^-[iSR + S,](ti.to)/h 



f/B=0 



'2'p JOf JOe 
"If /"I'n r2'a 

V'{R) 

/Of J2b JOb 

Combined with the current vertex insertions J dx2j22' • jiv of Eq. (jB.55ap . we 
obtain 

fdx2j22'-jiv J^3,2v{tiM-M) = fdx,^^^^pl^-^J 'i^V'{R) (B.81) 

J a—FB J Of ^ Ob 



I 



I 



April 10, 2008 10:57 WSPC/INSTRUCTION FILE ijmpb-review-vondelft 



54 Jan von Delft 



This expression, which is the first central result of our formalism, has a simple 
interpretation: thermal averaging with pjjg at time (— > —oo) is followed by prop- 
agation, in the presence of interactions (described by e~^^^"-^^'^), along all possible 
paths from time to up to time ti , with insertions of current vertices j {t2^ ) at time 
^2 on either the upper or lower Keldysh contour, and another current vertex j{ti) 
at the final time. 

Let us now determine the effective action explicitly, by using Eq. (|B.61cp for B 
to evaluate ■ ■ B: 



^ „ „./ Jto Jto 



(B.82) 



^FWiX 'BWfX + 2^34 ssWf:,^ (A 



sbW-! 



(^-1)^3", the inte- 



There are now two somewhat different routes to proceed, which lead to two some- 
what different (but equivalent) representations for the effective action. The first, 
followed in the present section, exploits symmetries under 3 4 write the effec- 
tive action in terms of as few terms as possible, leading to expressions [ (jB.83[) . 
(|B.84p . or (|A.7p . (|A.8|) ] useful for recovering the Keldysh diagrammatic results 
for the Cooperon self energy [ (jB.89p or (jA.lOp j. The second, summarized in Sec- 
tion IB.6.3[ does not combine similar- looking terms, and is useful for establishing 
contact with other, more standard influence-functional approaches. 

Let us proceed with the first route. Since {A'^^)'j^" 
grand in Eq. (|B.82p for B ■ A^^ ■ B is symmetric under the exchange of variables 
X]q / dt^dr^ ^ J2a' I dtidr^. We have exploited this fact to insert a factor of 
2^34 into the first and last terms of Eq. l|B.82p . which both individually have this 
symmetry, to obtain time-ordered integrals for these, which has the advantage of 
reducing the number of terms in subsequent expressions. (We could similarly have 
inserted 26*34 into the second and third terms of Eq. (|B.82p . too, but since only their 
sum has the above-mentioned symmetry, this turns out to be inconvenient.) 

More explicit expressions for Sr/i can be found with the help of Eqs. (jB.61d|l 
for W|", Eq. (|B.75P for A"^ and recalling that 6*34^43 = 0. Using the shorthand 
(iSii/Si) to present two similar equations in one line, and writing [iR/I)^^^ , = 
{iR/i)[tij,r^{ti) — r?'(tj)], where tij = ti — tj [and hkewise for C^^'^l^], we find: 



[tSR/Si]it,,to)[R''] 



E 

aa' 



dt:, 



dti^, (zL^/Z^)3„4„, [HI 



(B.83) 
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- (zL«/L^)3b4b 



The 6u functions in the second Unes of Eq. (jB.84[) wiU remove one of the coordinate 
integrations J dx^^i that are contained in the path integral / 2?'(i?). The second and 
third terms of Eq. (|B.82[) are equal, as can be seen by setting 3^4 and a a' in 
the third and recalling that {A~^)%^'^ = {A~^)'^2 ; we exploited this property above 
to combine those contributions from these terms that are proportional to -Rs^j^ [or 
-^3f4b] together into Eq. (jB.84bp [or Eq. (|B.84cp ]. hence the factors of 2R in these 
equations. 

Note that if we make approximation (ii) of Sec. IB. 5. 71 a useful simplification 
occurs [which was exploited in App. A to obtain Eqs. (jA.8|) from Eqs. (jB.84[) ]: all 
the factors (6*4^,2 + y"^ ^24^') above then reduce^ to 1, because y° ^ 1 was needed 
only to deal with interaction vertices occuring at times t^^, earlier than a current 
vertex on the same contour a', and these are precisely the ones that are dropped 
under approximation (ii). 

Eqs. (|B.84p for the effective action (iSi + Sr) constitute the second central result 
of this section. It should be emphasized that in the path integral (|B.8ip . the Pauli 



" To see this explicitly, we argue as follows, discussing in parallel the cases of J(S-=F/B)^ having 
a current vertex on the upper/lower contour and for which we have decided to use y<^=^/s = 0: 
if an interaction vertex lies on the same contour as the current vertex, i.e. for j(''=^/^) on the 
upper/lower contour at time t^^^^^ (hence a' = F/B), approximation (ii) says that it must lie at 

greater times than the current vertex, t^p^ ^ > t2p^g , implying that (84^^ g2+y^ ^ ^ 824^,^ = 1. If 
instead the interaction vertex lies on the opposite contour than the current vertex, i.e. for J(<^=^/s) 
on the lower/upper contour at time ^4^/^? (hence a' = B/F), the fact that + = 1 (always) 
and that we chose y"-^^/^ = 0, implying y^^^ = 1, also gives (64^yp2 + y^ ^ ^ 824^ ^ p) = 1, 
independent of the value of t4^^p . 



^034 spsp W[X i^R/Ih.i. (B-84a) 

1 . . ([S~ {042 + y^e24)2p%pip\ ~R/K 

\ SBSF (2^i^//)3B4. ^f^t^ (B-84b) 

1^. /[5- (042 +y^e24)2p°]4,4,\ ~,RIK J_ 
\^4p\p 



^3B4i. ^3B3i; 



= "i SFSB <^t. (2*^/^~)3.4b <b4b (B-84c) 

_ , 1 . ^ .A/K {[5 - (042 + y^024)2p"]4B4. 

- ±2* "SfSf ^4b3f 1 1? 

I- 2''4b4b 

= -034 SBSB {iR/ih^A^ W-l\^ W-ll^ (B.84d) 

1 ■ rA/K , - (042 + y^024)2p°]4B4, 

- T2^ ^iB3B '^3b3b I 
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principle is fully accounted for by the Pauli factors {6 — 2p) in S^-. The ability to 
incorporate the Pauli principle into an influence functional for interacting electrons 
may be regarded as one of the main achievements of the formalism developed so 
far. 

This concludes our rederivation of GZ's influence functional. In the remaining 
section B.6, where we show how it is related to diagrammatic Keldysh perturbation 
theory, and in the main text, where we use it to calculated the decoherence rate 7^, 
our analysis differs significantly from GZ's, since we come to different conclusions. 

Let us here just mention one such difference: According to the first lines of 
Eqs. (jB.84[) . iSn and Si are, respectively, purely imaginary or purely real functionals 
of the paths i2°, since W^"", Rij and ly are all purely real functions. GZ have used 
this fact to argue that after averaging e~(*'^«+'^^) over all paths [as required by 
the path integrals in Eq. (jB.80[) ]. e'^^'l^ will produce an exponentially decaying 
function of time and thereby determine the interaction-induced decoherence rate, 
whereas e"*"^^ will just produce an oscillating time dependence, and hence, quite 
generally, cannot contribute to decoherence; in particular, they argued that "iSn 
can never cancel any contribution from Sj" [discussion before (GZ-III.22)]. This 
general argument would work if the measure used in the path integral were real; 
however, it does not apply to the present case of Eq. (jB.80p where the measure 
giiSfi '' /h jg complex, since the average of a purely oscillatory function, using a 
complex measure, can well contain a decaying component, too. Indeed, it is shown 
in the main text [end of Sec. [4] that contributions from iSr and Si do partially 
cancel each other. 

B.6. Influence Functional vs. Keldysh Diagrammatics 

To check the general formalism developed above, it is important and instructive to 
verify that it can reproduce the standard results of diagrammatic Keldysh pertur- 
bation theory, before disorder averaging. We shall do this by expanding the path 
integrals (jB.80p in powers of the effective action {iSn + Si): 




(B.85) 



Now and henceforth using approximation (ii) of Sec. IB.5.7[ we shall use this ex- 
pansion to reproduce the Keldysh expressions for the conductivity in first order 
perturbation theory [Eqs. (|B.86[) ]. and to obtain general expressions for the first or- 
der contributions to the Cooperon before disorder averaging [Eqs. (jB.89[) ]. thereby 
reproducing the familiar Keldysh diagrams for the Cooperon self energy [Fig. lAll of 
App. A]. 



I 



I 



April 10, 2008 10:57 WSPC/INSTRUCTION FILE ijmpb-review-vondelft 



Influence functional calculation of decoherence in weak localization 57 




Fig. B2. Feynman diagrams for the first-order correlators j'^^^''^ of Eq. I|B.86| |. and for the 
vertices of Eqs. I|B.92|| . 



B.6.1. First Order Terms and Cooperon Self Energy E^/^ 

The iV = 1 terms of Eq. (|B.85p can be used to obtain the first order contribu- 
tions, j'^2''2V^^)^ correlators needed for crDc,rcai [Eq. (|B.55a|) ]. This straight- 
forward, if tedious, excercise is discussed in App. IC.3I The result can be written as 
•/i2'l2i'(0) = /(rfe)4'l2i'- where [see Fig.ll] 



?/(!)£ 

•^12', 21' ^ 


aa' 


(B.86a) 


j'WI) _ 




(B.86b) 


f(R/i) 


Gf/^(e - Co) GUe) Gi{e) G^Ae - Co){C^' /\C'')^,{Co) , 


(B.86c) 


f'iR/i) 


Gf3(e - u) Gl,{e) G^,{e) cf/^ie - j ^^{u^) , 


(B.86d) 


f{R/I) _ 
•^BB — 


G«,(e) G^e) Gf/^(e - c^) G^i,(e)(£^/£^),3(c^) , 


(B.86e) 


where j'}^^^^ 


denotes a first-order contribution from (iS^/S^), with interaction 



vertices that lie on contours a and a' . These expressions agree with those of standard 



diagrammatic Keldysh perturbation theory. Thus, the basic building blocks of the 
influence functional approach, including its treatment of the Pauli principle, have 
survived their first test. 

Next, we shall derive a general expression for the self energy of the Cooperon 
propagator. Usually, the Cooperon self energy is defined, after Fourier transforming 
to momentum space and disorder averaging, by a Dyson equation of the form Cq — 
Cq -I- CgSqCq, whcrc Cg, the free Cooperon in the absence of interactions, is the 
contribution to (G'''^G'"^)dis of time-reversed paths [cf. Eq. (|F.3bp ]. To identify a 
similar structure in position space and before disorder averaging, we need to write 
the first order {N = 1) term of Eq. ([RSS]) in the form U§ ■ t, ■ U§, i.e. a self- 
energy insertion sandwidched by two forward-backward propagators {Ub)YJi'b = 
ljOF,i^jpjjOB^ ^ f^2QR,i^jj.QA^_^ (^^^^Yi of which wiU producc a Cooperon upon 
disorder averaging): 

-If f V\R){rL^lPU'{hM) = (C/| • Sf// . U^)]^X (B.87a) 

-I JF -I ]B 
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= J dx,,dx,^ I dxi^dxs, iuE)YX {^al')'^EtZ ■ (B.87b) 

(For ease of recognition, we here and henceforth in this section write indices asso- 
ciated with the forward (f) or backward (b) paths as superscripts or subscripts, 
respectively). As made exphcit by Eq. (IB.87b[) . the first (or second) dot product on 
the right-hand side of Eq. (jB.87aP indicates integration over the two coordinates 
associated with the two "outgoing" (or the two "incoming" ) vertices at the corners 
of the self-energy box [see Fig. lAlb .]. Now, the left-hand side of Eq. (|B.87ap contains 
two vertices, associated with the indices of [iL^ / U)^^/^ , [Eqs. (|B.84p ]. as insertions 
into a double path integral, and therefore contains four Green's functions G [cf. the 
rule of thumb (jP.lOp of App. D.3]; however, for ■ E • i/^, we formally need six 
Green's functions G and four vertices, one for each corner of the self-energy box. To 
achieve this, we proceed as follows: the two corners to which the interaction lines 
are connected [black dots in Fig. lAlj can be naturally labelled by a and a', which 
take the values F/B, according to the contour that the corner sits on; let a and a! 
similarly label the other two, "free" corners [empty circles in Fig. lAlj . For the free 
corner a (and similarly for a'), we use the identity {tk is an arbitrary time between 
ti and tj) 

Gli^: = (sa^/^) ^ cf/,^ Kk. G^'J, . (B.88) 

taking R/A and Sg ~ ±1 if a = F/B, to write one Green's function as the con- 
volution of two, and regard the 5a function as the "vertex" at the corresponding 
free corner of the self energy box0 In this way, the self-energy contributions S^^/ 
are found to be given by the first lines of the following equations (summarized 
diagrammatically in Fig. lAip : 

V /4i53i5 h 

J (G^/«)3.3-G^^3^(£«/£^)3.4f ^ 89^) 

(w^+Wb^ -uE-h ~5^)[''' hR/if'i, 

^ ^ 4b3b 



' By using Eq. I IB.881 1 twice at the two free corners, an extra overall phase factor of (isa)(is^/) = 
—SaSa' is generated. The latter cancels the overall phase factor (— SaSa') occuring in the first 
lines of Eqs. l|B.83|l for — (iL^/L^), which is why this factor does not occur in the first lines of 
Eqs. llR89l l. 



V /4b3b 

V^^^J4b3b" 

V^^^J4b3b~ 
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= -f G«^^-^-(G^/^)4«3.(^~^/^^)4.3. . (B.89d) 

To obtain the second lines of Eqs. (|B.89[) from the respective first lines, we proceed 
similarly as for Eqs. (|B.84p [but now with 642 + 6*24 — 1, since we use approx- 
imation (ii), as explained in the paragraph after Eqs. (|B.84[) ]. In particular, we 
exploit the fact that the time-integrals in a path integral are time-ordered for the 
upper contour and anti-time-ordered for the lower contour to replace f/"^ by ihG^ 
and by -ihG"^ [cf. Eq. (|R46|], or, if they are pre- or post-contracted with 
(5- 2po), by ihGK [Eqs. (|B.47bp ]. For example, to obtain Eqs. (|B.89a|) and (|B.89d|) 
for flpp^QQ, we used: 

(w^+Sb -UE ■SbW^'Y'^^'' (B.90a) 

V /4B3i3 

= +l^inG^''-''-){-^hGl,^) , (B.90b) 
(s^Wb- ■ UE ■ W^+S^) [^^^ (B.90c) 

=j dx-,^dx4,J dXi.dx-,^ ^3-3- isB(l - 2poh,4s U°^-^'^' UZzs hs3s 

= -i(^;.G«•3-4.)(^ftGf^3J. (B.90d) 

Satisfactorily, the second lines of Eqs. (|B.89|) . summarized diagrammatically in 
Fig. IA1[ are identical to what one obtains from Keldysh perturbation theory, as 
can easily be verified starting from Eq. (|E.24p of App. IE. 31 Moreover, they are ev- 
idently consistent with the first order results listed in Eqs. (jB.86p above. (In fact, 
the latter could have been used to guess Eqs. (|B.89|) : the reason for nevertheless 
going through the above analysis was to check that the signs can be organized in a 
manner that allows for a series to be summed up.) In Sec. IF. 21 we shall calculate 
the Cooperon self-energy explicitly by starting from Eqs. (jB.87bp and (jB.89p and 
performing the disorder averaging diagrammatically. 

B.6.2. Fate of the PauU Factor {5 - 2p°) 

One instructive outcome of the analysis of the previous section is that we have 
learnt quite generally how to deal with the Pauli factors (S — 2p°) occuring in Sr: 
All Keldysh functions in Eqs. (|R86)) and (|R89)) arose from exploiting the identities 
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Ug{S - 2p%j = {d- 2p")aU^j = ihG^ [Eq. (|B.47bll ]. Since its frequency Fourier 
transform obeys [Eq. (IB.47bl) ] (s) = [l-2no{he)][G^-G^]{e), andinEqs. (iRSe)) 
and (|B.89[) all Keldysh functions come in the combination G^(e— tD)£^/'^((Ii), we can 
deduce a rule of thumb: by transforming to the coordinate- frequency representation, 
one generates the replacement 

{S - 2/5° ^ tanh 

Actually, in deriving the general structure of the self-energy above [Eq. (|B.89[) ] , this 
replacement has, in effect, already been deduced directly, and to all orders in the 
interaction, from the general form oiiL^^, in Eqs. (jB.84p . by exploiting the fact that 
in the path integral, each L^^, is sandwidched between propagators f/°. Since this 
point is so important, let us spell it out once more: depending on whether a vertex at 
time sits on the forward (time-ordered) or backward (anti-time-ordered) contour 
(a' = F/B), the factor {5 ~ 2p°)C^/^ occuring in L^^, is sandwidched as follows 
between two G^ . . . G^ or G^ . . . G^ functions [see bottom two diagrams of Fig. lB2| : 



■UJ 



IT 



(B.91) 



Gf,4.(^-2p°)4.4 

(5-2p°)4,4«G, 



Gf^iF - ^) 4\ i^) Gg,, (^")' (B.92a) 



A _ 

B^B "-^4b«£ 



-G 



FiF 

A _ 
Jb4b 



{e)CiJu>)G 



K 

4:BlB 



FJF ' 

{e-Co) 



(B.92b) 



Here the left- and right-hand sides are written in the time and frequency domains, 
respectively, and the replacement rule (jB.91|) follows from Eqs. (|B.92p since G^ {e~ 
oj) contains a factor tanh[ft(e — <D)/2T]. To be very explicit, the arrows in Eqs. (|B.92|) 
are shorthands for the following series of manipulations on the above factors of 
G^4p(^ — 2p°)4^4^ or {S — 2fP)igigG'^^^^ occuring on the forward or backward 
contours [indices are now dropped, for brevity]: 

G^(5 - 2pO) - [G^^ - G^]{5 - 2p") 
{5 - 2pO) G"^ = ((5 - 2pO)[G^ - G^] 



G^ ^ 



[G^ - G^] tanh ^ G^ tanh , 



[G^ -G^] tanh ^ G^ tanh . 

(B.93) 



Beginning in the position-time representation on the left hand side, we exploit the 
fact that the upper or lower contours are time- or anti-time-ordered to add an extra 
_GA/R = inside the square brackets, thereby obtaining a ±G^ . Step (1) indicates 
Fourier transforming to the position-frequency domain, in which the tanh factor be- 
comes explicit. (Step (2) will be discussed later below.) The expressions obtained 
after step (1) are the ones used to produce the right-hand sides of Eqs. (|B.92[) : satis- 
fyingly, the latter are precisely the combinations produced by the Feynman rules of 
diagrammatic Keldysh perturbation theory, illustrated in Fig. IB2I [As Eqs. (|B.89p 
show, the signs work out correctly, too, if the bookkeeping is done sufficiently care- 
fully]. The above argument is indeed completely general, and holds for each vertex 
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separately (but with different fte's at each vertex) , to all orders in perturbation the- 
ory. Thus, we have succeeded in recovering the Feynman rules from the influence 
functional approach. 

In Eqs. (jB.Qip and (|B.92p . the variable he represents the energy of the electron 
line on the upper (or lower) Keldysh contour before it enters (or after it leaves) 
an interaction vertex at which its energy decreases (or increases) by hiu [see lowest 
two figures in Fig. IB2| . The subtraction of uj in the argument of tanh thus reflects 
the physics of recoil: emitting or absorbing a photon causes the electron energy to 
change by HQ, and it is this changed energy h{e — uj) that enters the Fermi functions 
for the accessible final states. (A standard back-of-the-envelope argument for the 
origin of the Pauli factor, based on the availability of initial and final states, is given 
in MDSA-l^'^, Section V.A.) Of course, he will have different values from one vertex 
to the next, reflecting the history of energy changes of an electron line as it proceeds 
through a Feynman diagram. 

The final step (2) in Eqs. (jB.93P [not contained in Eq. (|B.92p ] indicates an 
approximation that occurs if one chooses to evaluate the path integral by including 
only time-reversed paths [as GZ do, see Sec.|l]of main text]: one thereby drops terms 
containing interaction vertices at which changes to on the upper contour, 
or G^ changes to G^ on the lower contour [so-called Hikami box terms], and thus 
drops G"*/^ tanh terms on the upper/lower contour. Of course, this last step (2) 
is optional: the Hikami terms can be retained, if one so chooses, and we do so in 
App. [R2l when diagrammatically deriving a Dyson equation for the Cooperon that 
includes the Hikami box terms. The result of that analysis is used in the main text 
[Sec. [5] to calculate the decoherence rate; remarkably and unexpectedly, it turns 
out that the Hikami-box contribution to the decoherence rate happens to be zero 
for the special form of the interaction propagator used in the main text, namely the 
unitary limit of Eq. (j4a[) . This fact implies that, for the specific purpose of deriving 
the decoherence rate (but not necessarily for other, more general quantities) from 
an influence functional, we may indeed adopt step (2) and drop Hikami-box terms. 
We shall do so henceforth. For the remaining terms, comparison of the very left 
and right-hand sides of Eqs. (IB.93|) clearly shows that one really can simply replace 
{S — 2fP) by tanh, without worrying about signs, etc., as specified in Eq. (jB.9ip . 

Having adopted step (2) of dropping Hikami-box terms, our rule of thumb re- 
placement (|B.9ip can quite easily be implemented "to all orders" in the infiuence 
functional approach: Fourier-transform the kernels {iL^ / L^)^^^ , [Eqs. (IB.84P ] of 
the effective action {iSr + Si) [Eq. (jB.83P ]. and simply make the replacement (IB.9ip 
in the Fourier-transformed version of iL^, now using the same energy he = he as 
that which enters the overall weighting factor [— 7ig(fe)]. The resulting form of the 
effective action is summarized in Eqs. ([2]) to (j4|) of the main text, which serve as 
the starting point of our calculation of the decoherence rate there. 

Diagrammatically speaking, the procedure just proposed amounts to using the 
same e inside each tanh[ft(e - u})/2T]C^/'^{uj). If one intends to consider only 
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self-energy diagrams and to treat infrared divergent frequency integrals with a self- 
consistently-determined lower cutoff {as GZ in fact do themselves in 
and as discussed in detail in Sec. H] and [5] of the main text) then this procedure 
would in fact not introduce any further approximations: the energy entering and 
leaving each self-energy insertion then is indeed the same for all such insertions, so 
they all should have the same tanh[;i(£ - u,)/2T] factors. 

Of course, once one includes vertex diagrams too, as needed if one wants to cure 
infrared problems "properly" (as in GZ99 ) instead of "by hand" (as in GZQii^, 
then the proposed procedure of using the same e everywhere amounts to a fur- 
ther approximation, since it neglects the accumulation of energy changes that are 
generated by vertex terms transferring energy between the forward and backward 
contours [as illustrated by the frequencies wi and u)2 in Fig. IBlj . Nevertheless, the 
mistake incurred by this approximation is insignificant, since the vertex terms are 
not ultraviolet divergent, and the frequency transfers contained therein are limited 
to the range h\u;\ ^ T, just as for self-energy terms. In fact, vertex terms become 
important only in the infrared limit where ~ 1/i (as required, of course, to cure 
infrared problems of the self-energy diagrams), so that we may replace cD by wher- 
ever else it occurs in a diagram. More formally, it suffices to treat the w-dependence 
explicitly only for that part of a diagram where it occurs as energy transfer, while 
Taylor-expanding in w all other factors of the diagram to which this w-dependence 
has propagated; only the zeroth-order terms of this Taylor expansion need to be 
retained, since the others contain higher powers of a) ~ 1/i, and hence produce 
contributions with a subleading time dependence. 

Note also that the accumulation of energy transfers manisfests itself only in 
diagrams of second or higher order in the interaction propagator. However, the in- 
fluence functional approach proposed by GZ and rcderived here features an effective 
action that is linear in the interaction propagator, and hence is equivalent to re- 
exponentiating the first order term in the expansion of the Cooperon in powers of 
the interaction propagator (as shown explicitly in DMSA-li2Ql) . Hence an accurate 
treatment of effects occuring only in second or higher order is beyond the accu- 
racy of the influence functional approach, in both GZ's original formulation and 
the modified version proposed here. The accumulation of energy transfers is such 
an effect. Fortunately, it only produces corrections that are subleading in time, as 
argued above. 

It is shown in the main text that if the replacement Eq. (jB.9ip is used in a 
"nonperturbative calculation" of a la GZ, a result consistent with conventional 
wisdom is obtained. Conversely, the reason why GZ obtained a different result is 
that they, in effect, omitted the —uj in the tanh-function in Eq. (IB. 911) . and hence 
lost the physics of recoil, as first suspected by Erikscn and Hedegard^. 

B.6.3. Alternative Representation of Effective Action 
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To facilitate a comparison of the infl uence functional approach developed in the 
present review with that of MDSA-i^Ql, it is convenient to rewrite the effective action 
derived in Section lB.5.81 and summarized in Eqs. (jA.7|) . (jA.Sp . in the following form 
(to be compared to Eqs. (21) of MDSA-l^O): 

[iSr + Si]{h,to) f'dh r dh ■ (B.94) 

aa'=F/B 

In particular, the integrands are to contain nonzero contributions not only for t^^ > 
[as is the case in Eqs. (|A.8|) ] but also for ^34 < 0. To this end, we follow the second 
of the routes mentioned after Eq. (|B.82p . We start from the latter, but instead 
of exploiting any 3 <-> 4 symmetries and inserting any factors of 2^34, as done in 
Section lB.5.81 ( "route one"), we now write out all terms explicitly, while still making 
approximation (ii) of Sec. IB.5.7[ namely to replace all factors of (^4^,2 + y"^ ^24^,) 
and {64 ,2 + y'^ O24 , ) by 1. [A perhaps quicker way to obtain the same results is to 
start directly from Eqs. (jB.83p . (|A.8|) . but to symmetrize the integrands w.r.t. 3^4 
by replacing J2aa' ^3o.4'a Saa' 1 [^^a-^'a ~^ '^4a3„/] ■] The result can be written in 
the form of Eq. (jB.94p . with £3% , being a shorthand for the following expressions: 



•'-'33,44 


= '^3f3f'^3^4f'^4f4f 


+ ^3.3.4>4.[^- 


2p°]4,4, 




2pO] 


3f3f''-3j?4f"4f4f 


nBF 
-'-'33,44 


= ^3B3B'^g^4F^4F4F 


+ '^3B3B^fB4.[^" 






2p« 


3b3b'^3b4f'^4f4f 


nFB 
•'-'33,44 


= '^3F3F'^g^4B^4B4B 




2P°]4b4b 


+ [5- 


2pO 


3f3f'^31'4b^4b4b 


fiBB 
•'-'33,44 


= <^3b3b'^3^4b'^4b4b 


- hB3B^5BiB t*^ ~ 


2P°]4b4b 


~[6~ 


-2p° 


]3b3b'^^b4b'^4b4b 



(B.95) 

[The double spatial indices, 33 for the forward and 33 for the backward contour, 
are associated with the same time and are both integrated over in the path 
integral (similarly for 44,44 and ti), see point (iii) after Eq. (jB.58[) ]. As explained 
in Sec. IB. 6. 21 upon Fourier transforming, the Pauli factors can be converted via 
Keldysh Green's functions into tanh functions. However, we now need to use a 
more general replacement rule (of which the one discussed in Sec. lB.6.2l was a special 
case), involving either of the expressions thzp = tanh[fi(e =F w)/2T]. The reason is 
that we now have to distinguish two types of vertices: For vertices of "type one" 
[Fig. iBSja)] , the arrows of the L^^^ and correlators that get generated both 
point in the same direction (i.e. both away from or both towards the same vertex), 
in which case we get the combination L^/'^(lj)G^(£ — ui): 

Ll-,J-2p%^4^^th^C^{Cu) , [S-2p%^,^Ll-,^,^th^C^{Cu). (B.96a) 

For vertices of "type two" (the occurence of which was studiously avoided in Sec- 
tion lB.5?8)) . the arrows point in opposite directions (one toward, the other away from 
the same vertex), [Fig. IBSf b)]. which gives the combination L-^/^{u!)G^ {e + lD): 

Ll-^J - 2p%^4B^th+C^'iLu) , [S - 2p\^-,^Ll-,^, ^ th+C^iu) . (B.96b) 
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44 £-C0 £ ^ 33 £+ 0) 





R^3 
CO 



£+C0 

K 33 A K 44 A 



B: 

(a) "type one" (b) "type two" 

Fig. B3. (a) Vertices of "type one" and (b) of "type two" arising in Keldysh perturbation theory; 
the accompanying Keldysh Green's functions are (e^^ui), respectively, producing Pauli factors 
tanh[(£ =f (Z')/2T] that dress the associated interaction propag ators Cg (Oi) and Cg (il>) [Eq. l|R96) l]. 



Using these replacement rules, the effective Fourier representations of Eqs. (|B.95p 
are readily seen to have the following forms: 

£3„4„, = J (dcD)(dq)e^^"[''"'*^")-''°'(*v)]e--(*3.-*v)£f (cD), (B.97a) 
Zfiu) = C^iu,) + Sa' th_,^, C^iio) + Sa th+,„ (u,) . (B.97b) 

Eqs. (IB.94|) and (|B.97|) together constitute an alternative and perhaps more compact 
expression for the effective action of Eqs. (2]) to ([3]). 

Appendix C. Relation between Path Integral and Cooperon 

In this appendix we show how the general path integral expression derived for the 
conductivity in the main text in terms of Ji2',2i' [Eqs. (|B.55a[) and (jB.80[) ]. can 
be rewritten in terms of the Drude conductivity cr^^'^'^ and the familiar Cooperon, 
and thereby clarify how they are related to the standard relations familiar from 
diagrammatic perturbation theory. We begin [Sec. IC.1| by reviewing the noninter- 
acting case before disorder averaging, then [Sec. IC.2| recall how disorder averaging 
produces the standard result for . Next [Sec. IC.3j we discuss the first order in- 
teraction contribution and subsequently [Sec. lC.4] generalize the analysis to include 
interactions to all orders, before disorder averaging. In particular, we elucidate how 
the average energy he of the two counterpropagating trajectories is fixed in this 
formalism. Finally, we perform a disorder average for the general case with interac- 
tions [Sec. lC.5] to establish a connection to the general Cooperon propagator in the 
presence of interactions, and [Sec. IC.6| review its structure in the coordinate space 
representation. 



C.l. Noninteracting Limit before Disorder Averaging 
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Let us check that in the noninteracting hmit but before disorder averaging, 
Eqs. (|B.55p for ctdCi with J'(0) given by Eqs. (|B.54|) and (|B.52|) . reproduce famihar 
expressions for the conductivity CTq^""*. If interactions are neglected, both and 
Ul^ in Eq. (|R52l) reduce to C/^. Using Eqs. ((R45l) and (|B.47aj) in Eqs. (|R52|, one 
then readily obtains 

•^i2',2i'(^i'*2;^o)= ^ J;[2'_'2i/(ii,t2;^o) = ^G'f2'G'^i/ + ^G^j'^^i'- (C.la) 

a=F,B 

Inserting Eq. (jC.laP into Eq. (|B.22a|) . we obtain a standard expression for crj?,™™', 
before disorder averaging. To evaluate its real part (Ji^c^^^ai [^Q- (!B.55a|) ]. we have 
to Fourier transform J according to Eq. (IB.54|) . Writing the result as J^^/ 21' ('^o) — 
J{de)j[f,'lj^,{ujo), we get 



Gf2,(e+)G<i,(£-) + G<2,(e+)G^i,(e_) 



(C.lb) 



with e± = e± ^ojq- Now expand J^°^-%ujo) = J("'^'^(0) + J' (°)'^(0), as needed 



for Eq. (IB.55al) . Using G^{e±) 



-no 



(fe±)[G^ - G^](£±), replacing G«M(e±) 



by G^/'*(£), and dropping terms in J'(''^''^(0) containing deG^/^ie), since they are 
smaller than those kept by a factor T/ep, we obtain 



"^12' '21' ('^) 

j;^?^^(o) 



-n',{hs)h'G^,,{e)G^,,ie) 



GUe)G^v{e) 



(C.2a) 
(C.2b) 



Here no(0 = 5jno(^), hence, in the J'('')'^(0) correlator of Eq. (|C.2bp . the energy 
argument he is constrained to be < T. The desired result for crg^'j."*^) of Eq. (jB.55aP 
thus is: 



nonint 



e) h [-n'oihe]] dx^ jiv •^22' h G^2' (e) G^r (e) 



(C.3) 



This is a standard result; it still has to be averaged over disorder, a step that we 
review in App. IC.2I 

The J(o)'^(0) correlator of Eq. (|C.2ap . in which the energy argument is not con- 
strained, turns out to cancel the (first) diamagnetic term in Eq. (|B.55bp . implying 
that cr^Q^^^g = 0, as expected. This cancellation can be verified, even before dis- 
order averaging, by using an exact identity, 

i / dx.jn' ■ [G^^/,^(£) J22' G^l^ie)] = "^cf/^e) , (C.4) 

proven below, to rewrite the contribution from J^^"^ (0) to Eq. (jB.55bp as follows: 

E A-2in'-i22'^^%^ = / (de)m{ns) [Gf,(s) - GA(£) 



ujom 



J2 ide)G<,ie) = 



(C.5) 
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which indeed cancels the first term of Eq. (|B.55bp . Since the DC conductivity is 
a real quantity, the latter cancellation of the two contributions to crPJ^g, namely 
the diamagnetic term and a term containing an integral J deriQ^he) over the entire 
Fermi sea, must hold order for order, to all orders, in perturbation theory in the 
interaction. Therefore, we shall henceforth not keep track of these terms, and take J 
to represent only those terms that end up containing a factor —nQ{he) that restricts 
e to the vicinity of the Fermi surface, as in Eq. (jC.2bp for J'^^\ 

It remains to prove Eq. (jC.4[) . It follows directly from another exact identity, 

d.uu' [GfJ%)Gf/^{e)] = - '-^Gf/\e) , (C.6) 

which can be derivecP^ before disorder averaging by evoking gauge invariance: let 
— {xi\X) and be exact eigenfunctions and eigenvalues of the single-particle 
Hamiltonian Hq [i.e. Ho\^) = '?a|A), cf. Eq. (jB.lSj) ]. and let A be a spatially uniform 
vector potential. Then the gauge-transformed wave- functions e~'"^^'^^^'^tpx{xi) = 

'ip\{xi) = {xi\X) are eigenfunctions of the gauge-transformed Hamiltonian Hq = 
Ho{P + eA) = Ha{P) + A-j + 4^ , again with eigenvalue ^a, i.e. Hq\X) = ^aIA). 
Consequently, the gauge-transformed version of G^^"^{e) can be written in two 
equivalent ways, as follows: 



Expanding both the left- and right-hand sides to linear order in A, and representing 
the latter in terms of the non-gauge transformed wave functions (x^jA) = ipxixi), 
we obtain 



A' A 



This readily yields Eq. I|C.6|) , since the matrix elements of the current operator are 
given by (A'|j|A) = Jdxijw [■ipy{xi')tlJx{xi)]. 

C.2. Disorder Average of Noninteracting Case 

Evaluating the disorder average {G^G^)dis needed in Eq. (|C.3p is a textbook excer- 
cise: Introducing an extra dummy integration Vol^^ Jdri into Eq. (|B.55ap . using 
Eqs. (jF.ip and (jF.3p from App. |F] and performing the momentum integrals using 
Eq. dElbl, we find: 

pp' 



Jp,p' 



Vol 2TTVT^{^/h 



(C.8a) 
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Fig. CI. Feynman diagrams for the correlators J^J^!^^ of Eqs. lIC.lOl l. 



Drudc 
'DC 



1 



deh[-n'oihe)] (dq) / dr C^ir) 



(C.8b) 



Here crj^™'^° = 2e^vD is the Drude conductivity and D = VpTc\/d is the diffusion 
constant. For the second term of Eq. (jC.Sbp . we introduced the variable q = p + 
p' and set q = everywhere except in Cq{uj = 0), because the latter 's infrared 

singularity as q ^ dominates the / (dq) integral. [Since ^p(O) from Eq. (IF.3bp 
has no singularities, its contribution to Eq. (jC.8b[) was dropped.] The J de integral 
in Eq. (|C.8b|) . which trivially equals one, is displayed here explicitly only for the 
sake of comparison with later results. 

The fact that the weak localization correction is small compared to the Drude 
term is often made explicit by expressing the prefactor of the Cooperon term in 
terms of the dimensionless conductance g^{L) [see Eq. ()B.9|) . and the discussion 
thereafter]: Using j (dq) = a'^^'^ J d'^q/ {2ttY for the momentum integral over the 
diffusive motion, and introducing, e.g., the dimensionless variables u = ti2/th and 
z = qLn (with Lh = ^/Dth) [if more convenient, e.g. in the absence of a magnetic 
field, one could replace th by here) we obtain from Eq. (|C.8b|) (times a' 



d-d\ 



nonint 
'd,T>C 



d'^Z 



9d{LH) TT J {2itY Jt^Jth 



(C.9) 



where we inserted an ultraviolet cutoff at small times, needed for d = 2,3. Appeal- 
ingly, the prefactor of the Cooperon term manifestly displays the smallness of Cj'^, 
via the largeness of g^. 



C.3. First order Calculation of J' 

In this section, we illustrate the structure of the perturbation expansion generated 
when the influence functional is expanded in powers of the effective action (iSR+Si), 
as m Eq. (|R85|: using approximation (ii) of Sec. IB.5.7| we explicitly calculate the 
first order contribution J'i2'\i'i^) ^'^ correlator of Eq. (jB.55a|) . i.e. the ujq = 
value of first derivative d^g of the tjQ-Fourier transform of 6*12 j}2' 21' (^ii ^2; ^o) = 

J2aaj[jia' + Jaa']- Here Jaf'^\= Y.~aJaa'^'^^ dcuotes the first-ordcr contribution 
to J'^^-' that arises from {iS^/S^) and has interaction vertices lying on contours a 
and a', while the index a in J^'^'i^^^'^ indicates which contour the current vertex is 
located on. 
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Our starting point is Eq. (jB.80|) . expanded to first order in ^[iSu/ Si)/h, using 
Eq. (iRSSl) : 

" J JOf J2b JOb J to ^ 3„4„' 

If interaction and current vertices occur on the same part (forward or backward) of 
the Keldysh contour, then, depending on the relative time orderings of the vertices, 
there can be more than one contribution to each of these quantities, which we shall 
denote by J^a''^^^\ with i = 1, 2, 3, etc. 

Consider j'^p'^^^^ [see Fig. IClf a)]. which has two interaction vertices on the 
forward contour at times and ^4 satisfying to < < < ti, and a current vertex 
on the backward contour at time ^2 satisfying to < ^2 < ti [in GZ's approach, who 
take to = t2, cf. Section [B. 5. 7[ these two sets of inequalities are replaced by a single 
one instead, namely ^2 < ^4 < < ^i]- Inserting Eq. (|B.84a|) for {iL'^ / U)^^/^^ into 
the first of the above equations, we obtain: 

JiBlM/I^ ^012 f\ f% f, fjo , fjO P-'^p'')i.iXr"./K 

'^ii 1^3f^3f3fC^3f4f| ^^J^ _ J^3f4f 

^^ipOp POfOb ^0B2'g ^2b1'b 

=-\ih^\ dtJ dt^G^^G^l^G^.G^yt^l^ . (C.lOa) 

Here integration over repeated spatial indices such as Oi? or Os or 3i? is implied; 
those over time are displayed explicitly, to keep track of the integration boundaries.] 
Eq. (jC.lOaP [whose index contractions are illustrated in Fig.|B2ja)] follows from the 
first line by relations such as Eqs. (jB.45|) and (jB.47[) (and dropping the subscripts 
F,B on indices). Moreover, taking the limit to —00 [but keeping ti fixed], the time 
integrals were extended to range over [—00, 00]. This is possible, since G^ contains 
a factor 9ij, G^ a 9ji, and £§4^ a 634, so that the product of Green's functions 
under the time integrals automatically vanishes for time arguments lying outside the 
integration ranges stipulated by the integration boundaries and ^-functions occuring 
in the first line. [However, if to had erroneously been replaced by ^2 in the first line 
above, as GZ do, the second line would have integration limits dt^^ dt^p , 
since G421 contains no ^42.] 

The case of J^'pp^^^^ is similar, but since both the interaction vertices at times ^3, 
^4 and the current vertex at time ^2 ah reside on the forward contour, three separate 
have to be considered [see Fig. IClT b) to (d)], corresponding to the three possible 
time orderings, namely (1): io < ^2 < ^4 < ^3 < ^i, or (2): to < ^4 < ^2 < ^3 < t-i, or 
(iii): to < ^4 < ^3 < ^2 < ti [since GZ implicitely take to = t2, the latter two cases 



J{F,R/I) ^ 
■^aa' — 



j{B,R/I) 
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do not occur in their approach]: 



FF 



t3f 

I dtAp I dXo^fig J7i^3^ ^3^3^ U'^p4p 



= I I dt^ ^13 ^42' ^21' 

J — oc J — oo 



(C.lOb) 



j{F2,R/ III 0l2 f 'j^^^l dti^ fdxop,OB ^1° 

'J to J to J 



to 



^■3j. O-ip^F '^3f2' ^2Fif 



dt- 



rii r<R nR n< rR/K 

ai4 Lri3 Lt32/ <_r24 U^^, 



'34 ^4f4, 



3f4f 



(C.lOc) 







to 

VO -0 



~"2^ ' " I " I "f:Ob '-^1f2'p '-'2f3f '^3f3f ^SpiF 



'34 ^4ir4, 



-3^.4^ 



^^4F0F^0F0Bf^0Bli 

/* 00 /* 00 

= —^ih^ I dt^ I dt^ G^2' ^23 ^ZA ^tv ^34 



00 



(C.lOd) 



00 



Eqs. (jC.lObp corresponds to Figs. IB2r b). The absence of a factor {5 — 2p°) in the 
first hnes of Eqs. (jC.lOcp and (jC.lOdp . and the corresponding absence of a G^- 
function in the respective second hnes, reflects the fact that we took — and 
that t4 < t2 in these integrals, so that the factor {642 + yF^242)2/5'' in Eq. (|B.84a[) 
for {iL^/L^)3p4p vanishes. Eqs. (|C.10c|) and (jC.lOdp are examples of contributions 
for which one or more interaction vertices occur on the same contour as the current 
vertex, but at earlier times. As discussed in approximation (ii) of Sec. IB. 5. 7) such 
terms contribute to "interaction corrections" but not to decoherence, and thus will 
henceforth be be excluded from our considerations. 

Adding the two terms [ (jC.lOap . (jC.lObI) ] that survive under the said approxi- 
mation (ii), we obtain jj/^ ~ J2d j'pp^ ■ The other three correlators, J^gp^\ 
•^FB^^ and J^^\ can be calculated in an entirely analogous manner. The results 
are: 



J 



FF 



jWI) 
"^BF ~ 



j{R/i) 

'FB ~ 



j{R/i)_ 
'bb - 



dt^ dt4 (§(3 G^l^' 



-00 
00 



\ih^ / dh dt4 Gf/^ 



-00 
00 



^42' ^21' + ^42' ^2V {^^ I ^^)34 ' (C.lla) 

G^i,(£«/i£^)3^, (C.llb) 



42' '-^23 + ^-^42' ^-"23 



dt^ dt4 



—00 
00 



r^R r'< A. r<< r<^ ,^k/a / fiA ii pK\ 

•-^32' '-^24 + '-'32' '-'24 ^-^41' \'~- I 2'- 



(C.llc) 



^^3 dt4 



r^R I r<< c'A 

'-'12' '-'24 + '-'12' '-'24 



GltGt^t^lt'')^^- (C.lld) 
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Satisfactorily, these expressions agree completely with those [Eqs. (|E.3ip ] obtained 
in App. IE.4j using diagrammatic Keldysh perturbation theory. 

To obtain J'{0), we have to Fourier transform these equations w.r.t. ti2, and 
then calculate Ja^^^\o) = [dujoJaa'^^^i'^o)]^ _q- For example, j'j^J^\ujQ) is given 
by 

Ji'^fHcoo) = -^^h l{de){dCo)GUe+)G'^i''{e+-Co)L^'r{Co) 

G%{e+)G<,,{e^) + G<^,{e+)Gi,{e.)\ , (C.12) 



and J'p^^^\q) is easily calculated by noting that the factor in the second line of 
Eq. (jC.12p equals J42/2i'('^o) [cf- Eq. (jC.lbp ]. whose first derivative is given by 
Eq. (|C:2bl) . namely J42?,2i'(0) = -n'o{h£)n^ G%{£)G^y(e). Thus, the final resuh 
for j;^^/')(0)is 

j;W^)(0) = -\^h^j {de){dLo) hn'o(fe)] CUe) G^I'^ie-u^) 0^^) GUe)cTi^)- 

Similar expressions for the other contributions j'^^^'^^ (0) can be derived from 
Eqs. (jC.lip in an entirely analogous manner, and are given in Eq. (|B.86p of 
App. IB. 6. 21 In each case the combination hlGf^'G^ + G^/G^] produces a factor 
h^hn',ihe)]Gg,{s)G^^{s). 

Actually, it is clear from the above derivation that in every order of perturbation 
theory in the interaction, such a factor will be produced for all terms that survive 
the abovementioned approximation (ii): in analogy to Eqs. (jC.lOaP and (jC.lObp . it 
will arise from a factor 

~ Jd-XOF.QB ^ipOp POfOb ^^B2'g ^2b]b ~ ^iF2'j, ^^2^0^ POfOb^^bJb ' (CIS) 

where tip and tj^ are the times of the earliest interaction vertex on the upper or 
lower Keldysh contour, respectively. 

To conclude this section, we wish to emphasize once more the significance of the 
fact, illustrated by Eqs. (jC.lip but valid for all contributions to J^^^ (including the 
"interaction corrections" ) , that all time integrals occuring in Keldysh perturbation 
theory can be extended to range over the entire real axis. Importantly, this implies 
that the Fourier transforms that are needed to obtain J^^\ujq) (and from there the 
conductivity) are always given by simple convolution integrals, such as Eq. (|C.12|) . 
In contrast, in GZ's calculations all time integrals J dt^ dt^ have t2 as lower limit, 
see e.g. (GZ-III.A20) and (GZ-III.A23) in GZofl, whose t' corresponds to our t2. 
This means that instead of obtaining simple convolution integrals, they erroneously 
end up with sin and cos functions, see (GZ-III.58) and (GZ-III.61). This leads to 
numerous incorrect complications and conclusions, such as the claimed existence of 
an "oscillating cos-term" in (GZ-III.70). Thus, GZ's perturbative analysis in Sec. IV 
of Gzod^l, in particular their discussion of the "breakdown of the Fermi golden rule 
approximation" in Sec. IV. B, is invalid, since its starting point is based on the 
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replacement to — s- t2, which is incorrect (and unnecessary, since the correct limit 
to — > — oo can be incorporated into GZ's approach, as emphasized in Sec. IB.5?7l and 
illustrated explicitly above). 

C.4. Thermal Weighting and Path Integral, before Disorder Averaging 




Fig. C2. (a) Diagrammatic depiction of Eqs. l|C.18| l for Jj|, or Eq. ||C.29|| for cr^c ■ Be- 
fore disorder averaging the black box represents P^^, (e — io); — a)), thereafter it represents 

Cq ^ {—1^, i^) I (27r!^Tei'^//i). (b) Real-spacc depiction of a typical pair of Drude (dashed) and time- 
reversed (solid) trajectories contributing to pi^.Drudo pviy^lj ^ corresponding to Eqs. l lC.31al l 
and IIC.31bl l. respectively, (c) Definition of variables used for Fourier-transforming the double 
path integral P^g (£", f22) of Eq. 1 IC.I6I 1. In (c), frequency and momentum variables are cho- 
sen such that Oi and O2 Bxe, respectively, the outgoing and incoming "Cooperon frequencies" 
(i.e. frequency differences between upper and lower lines); q ± P3 are the outgoing and incoming 
"Cooperon momenta" (i.e. sum of momenta of upper and lower lines); E it ^3/2 are the average 
(between upper and lower) frequencies flowing out of our into the Cooperon, respectively. The time 
variables ri 2 and fi 2 and coordinate variables pi_2 and pi 2 are purposefully defined in such a 
way [Eqs. IIC.17al l. ||C.32| |] that the Fourier exponents in Eqs. l lC.16bt. IIC.16cl l and IIC.351 1 are free 
of factors of 2. (Our labelling convention differs from that of AAl4ifi( which has typos involving 
factors of 2.) 



The presence of interactions will, in general, modify the result (jC.Sbp for a\ 
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in two ways: firstly, it can renormalize the value of ctdc"^"' this effect is not in- 
teresting for present purposes and will be ignored here. Secondly, it can reduce the 
life-time of the Cooperon propagator, thereby contributing to decoherence, which is 
the effect we are interested in. Our goal in this section is to express the conductivity 
of Eq. (|B.55aP in terms of double path integral expressions for J[2i 2i'i^)^ obtained 
from Eq. (|B.80p . in a way that is generally valid in the presence of interactions, 
before disorder averaging, and properly accounts for thermal weighting via a fac- 
tor [—nQ{he)], as in Eq. (jC.8b|l . Hence, we will have to find appropriate Fourier 
transforms of our path integral expressions that relate them to the energy he. 

An important first clue comes from the first order relations (jB.86|l for J^j/ 2i'(0)- 
each term contains a factor J {de)h'^[—n'Q{He)]G^,{e)G2j{£), thus the current vertex 
j22' is always sandwiched between a retarded and advanced function with energy e, 
^i2' (s)j22'G'^ , and thermal weighting is always governed by a factor [— nQ(?ie)]. As 
explained after in App. IC.3I [just before Eq. (jC.13|) ] these properties actually hold 
in every order of perturbation theory in the interaction, for all terms that survive 
approximation (ii). Of course, the other current vertex j'n/ is similarly sandwiched, 
too, but in general with a different energy argument, G^,(e — cj)jii/Gf|(e — uj). 
The general expression for that part of the conductivity containing the Cooperon 
propagator, relevant for weak localization, is by definition the sum to all orders of all 
such terms containing [— »^o]G'||,j22'G'^- . . . Gj^/jn'C^j. In path integral language, 
it will thus have the following form. 



written in analogy to Eq. (jB.55a|l for crDc.roai, where the integral equals J^j/ 2i'(0); 
and J(|/ equals [—nQ{hs)] times some suitable frequency Fourier transform 
(needed to set the energy to e) of a double path integral whose forward path con- 
nects the points r2' and ri, while the backward path connects r2 and ry. To find 
the appropriate expression, we begin by considering the general double path integral 



depicted schematically in Figs. IC2F b) and IC2r c). It ranges from rf^ at time to 
rf at time tf (> on the forward contour and from rf at time tf to at time 

(> tf ) on the backward contour. These times are understood to be the limits of 
the J dta time integrals in and {iS^ + Sj), and tf , tf are in general not equal to 

tf , since they will have to be Fourier transformed independently [as required, 
e.g., to properly define the variable e in Eq. (|C.14[) ]. For general time arguments, 
we adopt the following conventions, depicted in Fig.[C2jc), for Fourier transforming 
from the time to frequency domain and back: 




(C.14) 




(C.15) 




{d£){dni){dn2){dQ3) 2TTS{n3)rl3i£;ni,n2) expi \ -t( £ + 
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£ - 



+ ti 



E- 



{d£) {du) {du') Vli{£; UJ + \J, UJ - \uj') 



Pli{8;ni,n2) = J dT,dT2dfi2 e*["^''^-"^"^+^^^^lPil(n2,n2,Ti2) . 



(C.16a) 
(C.16b) 
(C.16c) 



(For P, the indices 43 stand for both coordinate and time variables, for its frequency 
Fourier transform V, distinguished from the former by using calligraphic script, they 
stand for coordinate variables only; we use a similar convention for the Cooperon, 
(7 or C, defined below.) For Eq. (|C.16b[) we changed frequency variables to w = 
i(rii + 1^2) and uj' — — fl2, and introduced various sum and difference times [see 
Fig.Elc)]: 



Tl = 
T12 



Tl - T2 



t 



t 



B 



T12 



2 

Tl - f2 



T2 —^2 + I3 , 



T2 



_ '•2 



T12 



Tl + T2 



(C.17a) 
(C.17b) 



On the right-hand side of the back transformation ()C.16c|) . ^43 (ti2, T12, fi2) by 
definition is given by P4I of Eq. (|C.15|) . with the understanding that the indices 
12,43 now only specify the path end points rf,r|^, rf,r^, but that the time 
arguments tf , if, if, in Eq. (|C.15|) are chosen such that Eqs. (IC.17|) hold. 

The frequencies introduced in Eq. ()C.16P have evident physical interpretations 
[see Fig. IC2f c')]. The "Cooperon frequencies" fli and Q2 are the outgoing and 
incoming frequency differences between upper and lower lines, respectively, while 
£ ± ^fia are the average (between upper and lower) frequencies flowing out of or 
into the Cooperon. In general, the presence of external time-dependent fields would 
require fl^, the total frequency difference between outgoing and incoming lines, to 
be nonzero. However, for the present purpose of calculating the conductivity in lin- 
ear response, such external fields can be set to zero; hence in Eq. (jC.16a[) we use 
a delta-function to set ^3 equal to zero, thus recovering translational invariance in 
time for P^^. 

Having identified the meaning of the frequency arguments fli and fl2 
[Fig. IC2f c')]. and inspecting the frequency labels of the standard diagrammatic de- 
piction [Fig. IC2l ^a). where an integral over the "internal" frequency uj is implied] of 
the current-current correlator needed for the conductivity, it becomes evident that 
the average frequency is £ — e ~ ^uj, while the outgoing and incoming Cooperon 
frequencies are fii = —uj and CI2 = uj, respectively (i.e. uj = and uj' = —2u). 
Moreover, the upper line runs from r2' to ri, while the lower line runs backwards in 
time from ry to r2. Thus, the particular Fourier transformed version of P needed 
for J[^, 2v in Eq. (|C.14|) is 



n'(fe)] / [d2Zj)V\l 



12' , 



-UJ,UJ) 



(CIS) 
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To check that that the normahzation factors and frequency assignments are cor- 
rect, let us expand P21, in Eq. (jClSp to zeroth and first order in the interaction in 
order to calculate J(|/ 21' = [J'^^-*'^ + J'^^-* '^]i2'.2i' , and compare the results to our 
previously-obtained expressions for these [Eqs. (|C.2bp and (|B.86p ]. To this end, we 
begin with Pg^, as given in Eq. (|C.15|1 and with general indices, expand it to first 
order in — [iSji + Si]/h, and express the resulting terms in terms of G^^^^^ func- 
tions. The details are analogous to those presented in Sections IC . 1 1 and IC .31 to derive 
jj^) 21' and j'^, from Ji2'%v of Eq. (jB.SOP (except that the latter's first line is 
not needed for P^^, and the limits of the path integrals are different). The result 
can be written as P^^ = P[%, + Y^aa' \Pia' + PaJ] , . where P^^^e^ and P^^/'^ 

L J 12,65 

are given by Eqs. (|CTa| and (|ClT|l . respectively [with (12', 21') (12, 65)], except 
that all occurences of the combination h[G^,G2j + G^/G^] have to be replaced 
by h'^G'^G^y Fourier transforming the result for Pg^| [via Eq. (jC.16cP ] to obtain 
■PgKf ; Oi, ^2), specifying the spatial indices as (12,65) —t (12', 21') and then inte- 
grating as stipulated in Eq. (jC.lSj) . one recovers J(|, 21' — [J'^^^^^ + J'^^^'^]i2',2i', 
with the first and second terms given by Eqs. (|C.2bp and (IB.86j) . respectively, 
as expected. Thus, our check worked. [Also, the reason for the 2 in J{d2uj) in 
Eq. (jC.18[) becomes c\eai:VQ^{£;il 1,^12) turns out to contain factors of 27r(5(rii — f22) 
or 2ttS{^Ii — + ■ ■ ■) for self-energy and vertex terms, respectively, which under 
the integral J {d2uj)V2ii {£ — \^\—^^^) of Eq. (jC.18|) have to collapse to unity, 
/(d25)27r(5(...-2(i) = 1.] 

Finally, let us rewrite Eq. (|C.18p in a more suggestive form. Transforming back 
to the time domain using Eq. (|C.16cp and writing the result in terms of the time 
variables of Eq. (jC.17bp , we find 

J[l,.2v = [-n'{ns)] / dT,2ply''{Ti2), (C.19a) 
P^^','%Ti2)^ dTi2 {d2^) e-'^^^'^ dfi2e*^i^(^-5=^) (C.19b) 



^P21' (t12,Ti2,Ti2) , 

/oo 
dfi2e'^i^^F2r'(n2,-iTi2,fi2) . (C.19c) 
-00 

We need to consider P21, only in the limit r2 ^ ri, since the Cooperon con- 
tribution to it is negligible for |ri — r2| > Ap, where Ap is the Fermi wavevector 
(assumed to be much smaller than the mean free path, Ap «C lei)- The purpose of 
the time integrals in Eq. (|C.19bp is to project out from the general path integral 
P21, of Eq. (|C.15p . defined in the position-time domain, an object depending in 
an appropriate way on both the average propagation time T12 of the forward and 
backward paths and the energy e occuring in the thermal weighting factor. (The 
simultaneous specification of both a time and an energy does not violate the time- 



energy uncertainty relation, as incorrectly argued by G2^, because p2i,''^{Ti2) is 
constructed from two electron propagators, not one). To see how this projection 
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works in detail, we use Eqs. (jC.17bp to write the time differences ri2, ti2 and fi2 
as follows: 

ri2 = i [(tf - tf) + {ti tf )] , n2 - i [(tf + 1^) (if + tf)] , 

ri2-(tf-tr)-(if-tf). (C.20) 

The / c?fi2 integral in Eq. (jC.lQbp fixes the average energy of the upper and lower 
electron lines (in diagrammatic language) to be e — lo/2 [where fi2 is the length 
difference between the forward and backward pieces of the contour]. The J{duj) 
integral averages over all possible frequency differences uj between the upper and 
lower electron lines, as is necessary when vertex terms are present that transfer en- 
ergy between them. And finally, the / dTi2 integral projects out the ri2-dependence 
of P211 [where T12 is half the difference between the midpoints of the forward and 
backward pieces of the contour]. The only remaining time variable, T12, is the av- 
erage of the lengths of the forward and backward pieces, and can be viewed as the 
"observation time" as a function of which P21, {T12) will decay. P21, {T12) will con- 
tain a contribution resulting from time-reversed paths that corresponds to the full 
Cooperon in the position-time representation, Cp=o(''"i2)- The time scale on which 
it decays is the desired decoherence time t^. 

Now, the / (duj) integral in Eq. (jC.lQbp yields 6{t + jf) [here and henceforth we 
drop the subscripts on r, r and t], leaving us to consider a path integral with time 
arguments iji' (''"i^i^i^); ^ indicated in Eq. (jC.lQcp . These time arguments can 
be obtained by choosing, e.g., ti = = \t and t2.4 — — ^(r ± f), resulting in: 

Pllr'{T)= dfe'^'W i V'{R)e-^'^"+^'^'^ . (C.21) 

Eq. ([CHll . together with (|C.19aP and JOlI]), are the central results of this section, 
because they express the conductivity in terms of a general path integral influence 
functional, with thermal weighting taken properly into account. The main difference 
to the path integral (|lb|l used in the main text (and by GZ) is that in Eq. (jC.2ip 
the duration of the forward and backward paths differs by a time f that is being 
integrated over in Jdf e*"^*^. The remainder of this section is devoted to justifying 
the replacement of Eq. (jC.21|) by the simpler Eq. ([Tb|) . 

The combination J de J df of integrals from Eqs. (jC.14[l and (jC.2ip have the 
effect of fixing the average energy of the forward and backward trajectories to be 
close to the Fermi energy, with energy spread of roughly ±T (in a way reminiscent 
of App. B of the review^ by Chakravarty and Schmid). To see this, consider first 
the noninteracting limit (i.e. ignore iSu + Si) in the semiclassical approximation, 
where the path integrals in Eq. (|C.21[) are restricted to all possible classical for- 
ward and backward paths f^/^its) having the specified boundary conditions, with 
corresponding classical actions S^^^f ~i "F §)■ Since these paths follow diffusive 
trajectories through a disordered potential landscape, for any given r and f the 
path integral still includes many such classical paths, with a range of different clas- 
sical energies (and correspondingly different diffusion constants). Now, the energy 
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integral in Eq. (ITa|) restricts the / df integral in Eq. (IC.21|) to the range |f | < h/T, 
since 

dehn'(fe)]e-^ = . if^^-rr., ■ (C.22) 
sinh(7r/irJ ) 

The relevant values of f are thus much smaller than the typical propagation times 
T relevant for determining the decoherence time [r ~ Ty ^ hg{L^)/T 3> h/T, see 
Eq. P^ ]. so that the classical actions can be expanded^ to first order in f, 

qF/B/T T r ^ ^ qF/B,r t ^ -r l^c-f^Z-^ /p oox 

■^O.cl V2' 2 + 2/ — "^0,01 ll' 2 + 2^^cl ' [^.ZOj 

where £^/^ is the classical energy at the endpoint of the corresponding classical 
path ^"^^^(^3). Using this in Eq. (jC.21[) . the / df integral is seen to fix the average 
classical energy of the forward and backward classical paths to be close to the Fermi 
energy e^^ = 0, with an energy spread of order T: 

nhfT 
sinh(7rfi,fr) 

(The right-hand side follows from using the integral representation (jC.22p for the 
sinh-function.) Note that the energy spread is consistent with the time-energy un- 
certainty relation in the limit of present interest, tT 3> h. 

Now, in the absence of interactions, the only effect of fixing this average energy 
e to be roughly ej? is that the velocity appearing in the diffusion constant is the 
Fermi velocity, D = VpTei/d. However, in the presence of interactions, the energy 
e also plays a role in determining the phase space available for electrons to get 
scattered upon absorbing or emitting a noise quantum. In particular, in perturbative 
calculations it shows up in the tanh[fi(e =F tD)/2T]-factors of the Keldysh electron 
Green's functions G^(e T '^). In our influence functional approach this can be 
kept track of by replacing Eq. (|C.21|1 by Eq. (jlbp . which mimicks the effect of 
the former's integral J dfe'^"^ by using (i) forward and backward paths of equal 
duration r and (ii) an effective action whose time integration boundaries are fixed 
at ±r/2, but which depends explicitly on the average propagation energy e. Note 
that GZ's approach in effect employs the same simplification, since they likewise 
have no / dfe^^'^ integral and use forward and backward paths of equal duration r. 

The e-dependence of the effective action enters through the Pauli factor {5 — 2p) 
occuring in Su [Eqs. (jA.8[l or (jB.QSP J. which we treat differently from GZ. In our 
approach, it produces factors of tanh[?i(e =p uj)/2T] in the frequency representation 
of Sr [ci. Eqs. (|4e| or (|B.97P ]. chosen in such a way as to be consistent with Keldysh 
perturbation theory, as discussed in Sec [3] and (more extensively) IB.6.2[ IB. 6. 31 In 
GZ's approach, the tanh-argumcnts contain e instead of e =F tD (i.e. their effective 
action depends on the average energy too). However, lacking the =p(D recoil shifts, 
the tanh-terms turn out to yield zero after averaging over random walks, so that 

The strategy just described for arriving at forward and backward paths of equal 
duration is of course not exact; but it is sufficiently accurate for our purposes: the 
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errors incurred by it are of order h/iTr) (<C 1 for r ~ t^), as can be shown by a 
detailed comparison with Keldysh diagrammatic perturbation theory (App. B.6.2 



of this review, and App. A. 3 of DMSA-Il 

C.5. General Cooperon, after Disorder Averaging 

Let us now disorder average Eqs. (|C.18[) . in order to arrive at an expression for 

in terms of a general Cooperon propagator, in the presence of interactions. To this 

end, we have to Fourier transform from position to wave-number variables. 



Vol^ ^ 
pi,P2,P3,q 







r 2 J 





-Pl 



q + P3 



- r. 



-P2 



q P3 



(C.25) 



as depicted in Fig. IC2r c). (Again, the dp.^^o guarantees translational invariance.) 
According to the standard diagrammatic approach for disorder averaging [cf. Fig. [FT 
in App. |F], the disorder average of 'Pq^pllp^ can be separated into a "Drude" and 
a "weak-localization" contribution, 



\ q,pi,P2 / 



dis 



+ir!i)^fg_p^(£-ir!i) |2^(f]i - n2)6p,,p. 



Vol 2lTVT^{^ /h 



(C.26) 



where in the second term the contributions from the four external electron lines 
were separated and a conventional prefactor {2-KVTe\^ / h)~^ was split off. The nor- 
malization of the general Cooperon in the presence of interactions, Cq(f7i,572), is 
fixed by requiring that when interactions are switched off, it reduces to its free 
version, Cq(rii), according to 

1 



no int 



(C.27) 



Just as Cg(ili), the full Cooperon C^{yii,Q,2) does not depend on the external 
momenta Pi,2, because, in diagrammtic terms, it is separated from external lines 
by impurity lines. 

In a purely diagrammatic approach, where one typically works exclusively in the 
wavenumber-frequency domain, Eq. (jC.26P would be the standard starting point 
for further calculations. Since the dominant contribution to C^{Vli,Vl2) typically 
comes from small q (with ql^x ^ 1) and small f2i_2 (with Oi^2''ei ^ 1), while £ is 
likewise small (< T), it is customary to neglect the terms i^q and £ ± \^i,2 in 

the arguments of the external Q functions, which simplifies the J dpi^2 integrals. 
To explore the effects of interactions, one would proceed to expand Cq(rii,il2) in 
powers of the interaction propagator, etc. 
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Instead, here we shall use the general Eqs. (|C.25P and (|C.26|) for Vl^{£; Qi,il2) 
to analyse the general structure of the disorder- averaged version of Eq. (|C.14[) . as 
needed for {cr-Dc}^^^- As intermediate result, we obtain 



dx2jiv ■j22'{r^l{£;n,,n2))^,^ (c.28a) 



4m2Vol 

piP2q 



■ ~}2q 

27T{n,-n2)^^ f{dq)cl{ni,n2)\ , (c.28b) 

where Eq. (jF.6bp was used (under neglect of £ ± 5^^1/2 in the frequency arguments 
of all electron Green's functions) to perform the momentum integrals, i.e. the / (dq) 
integral for the Drude contribution to 'Pq',pl'p2 1 and the J{dpi) integral for the 
Cooperon contribution (for the latter, the ^q arguments in the external electron leg 
Green's functions were neglected). Inserting Eqs. (|C.18|) and (jC.28bp into Eq. (jC.14[) . 
we readily find: 



^ ^Drudc 
O^DC — O-DC 



l--i^y deh[~n'„{he)] J (dq) y(d25)Cq ^"(-w,w) 



(C.29) 



Eq. (|C.29|) is the desired generalization of Eq. (jC.8bp [and in the absence of inter- 
actions, duly reduces to the latter, via Eq. (|C.27|) ]. 

C.6. Cooperon in Position- Time Domain 

For our present purpose of relating the diagrammatic and path integral aproaches 
to each other, it is instructive to understand the consequences of Eq. (jC.26p also 
in path integral language. To this end, let us transcribe Eq. (|C.26|) back into the 
position-time domain, in which the Cooperon is defined as: 

C^(ri,T2) = J{dq){dn,){dn2)e'^P'i-'''^'+''-^-^Cl{nu^2). (C.30) 

Inserting Eqs. Kim and (IU:26l) into Eq. (IC.16al) yields (P)^.^ = pDrudo ^ pwL^ 
with 

pi2,WL f dfi df2 dii di[ dk di'2{d£) n.r p^u: 7m\ 

xh^g^^_,^ ih - h) Qt-r, iu - i\) gf^_,^ {t2 - 12) gi^_,.^ {P, - 1,) . (c.aib) 

Fig. lC2r b) offers an intuitive interpretation of these expressions: p^^^^'"^'^° gives the 
amplitude for propagation from (r2,i2) — * (''1,^1) (forward in time) times that for 
('"37 is) ~* {fij ^4) (backward in time). And P43 ' gives the amplitude for forward 
propagation from (r2,i2) — * (^^2,^2) ^ (ri,ti), times that for backward 

propagation from (r3,<3) (7^2, ^2) ^ {''^iii'i) ^ ('"4:i4)- The middle part of the 
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forward and backward paths have the same beginning and end points in space, 
albeit not in time, and hence can interfere constructively if the paths connecting 
them are time-reversed partners. 

The approximation mentioned above of neglecting ±i<7 and £ ± ^i^i,2 in the 

arguments of external Q functions has a counterpart in the position-time domain: 
when performing the integrals in Eq. (|C.31bp . it corresponds to exploiting the fact 
that Gr{i) has a short range in space (|r| < Tei) and time (|t| < lei) [cf. Eq. (jF.2dp ]. 
To be explicit, the latter fact means that the disordered Green's functions occuring 
in the second line of Eq. (|C.31bp act effectively as delta-functions in time as far 
as the factor e"*^'- 'Cf^^^,^ ( ) is concerned. Thus, in the latter we may make the 
replacements ti — > ii, t'l — > ^4, t2 — > ^2, ^2 ~* ^3' softer which the four time- integrals 
each yield a zero-frequency Green's function, J dtQr{t) = Qr{£ = 0). Introducing 
the sum and difference coordinates 

_^^r^^^ p,^rf-rf, ^ !1±^ , P2 ^ r[ ^ , (C.32) 

recalling similar definitions (|G.17p for the time variables, and shifting the space 
integrations according to fi^fi + pi for i — 1,2, Eq. (|C.31b|l gives: 

312,WL _ L^; f I ^li -P2+ri-.^2 (^1 ' ^2) 



xh^g? - (0)^-^1 - (0)^^ 1 (o)^^^i (0). (C.33) 

Since the zero-frequency Green's functions ^^(0) decay with distance as e^l'''/^'"' , 
we note that ~ ~ — f^, which implies that \fi\ < l^i and \pi\ < l^i- Thus, we 
may drop the terms f i — f2 from the argument of in Eq. (IC.33P , whereupon the 
spatial integrations can be performed explicitly, using 

dU r di, f dri gf^'^{U)gtJ''it^) = Mil , (c.34a) 



(C.34b) 



where {r^ ) is a "smeared-out delta function" of normalization J dr^ Si^^ {vy ) — 1 
and width ~ l/fep, the Fermi wavelength (since l/fcp ^ ^oi, the width is set by 
the oscillating factor sin(fcFr)/r, not by the exponential e^''/^'°'). Thus, Eq. (|C.33p 
becomes: 

^4"'^' = ^kMKM J {d£) e-^^(^-^^) ^tt^^n'^ ■ ^^-^'^ 

This useful result clarifies the relation between the coordinates 1,2,3,4 of ^43'^'^, 
and the times and spatial coordinates relevant for the Cooperon. In particular, we 
see that P43' is nonzero only for \pi\ = \ri — r^l < l/fcp and |p2| — |t'2 — 
rsl < 1/kp. Moreover, if we want to describe a Cooperon with a specified average 
energy £, we need to Fourier-transform P^^'^^ with Jdfi2e'^^'^^'^ . Note that for 
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P2i''^^i as needed in Eqs. (|C.19|) . the Cooperon position argument is identically 
zero, pi — P2 = 0, while |pi| — \ri — ensures that ri and r2 lie close together. 

Appendix D. Time-Slicing of Path Integral for IJ^ 

In this appendix, we give an exlicit time-slicing definition for the path integral rep- 
resentation (jB.56P of the propagators used the main text, and derive various 
properties thereof. Our discussion is very (perhaps overly) detailed, since the ob- 
ject of interest is somewhat unconventional, namely a path integral for a non-local 
hamiltonian. We begin [Sees. D.l to D.3] by defining it in terms of a path integral 
/ DR J DP over paths in both coordinate and momentum space, which is the form 
used by GZ; then [Sec. ID.4j we explicitly perform the J VP integral to arrive at a 
"coordinate-space-only" path integral J VR, which is the form used in App. B.5 to 
B.8. Finally [Sec. ID. 5] we present explicit expressions for the effective Hamiltonian 
in the position-momentum representation used by GZ, and [Sec. ID.6| recover 
from this GZ's expressions for the effective action {iSn + Si)[R°', P"]. 



D.l. Time-Slicing Definition 

The propagators U^j are defined by the requirement that they have to satisfy both 
the conditions Eqs. (jB.39|l . This fact can be used to give meaning to the formal path 
integral of Eq. (jB.56[) . by using the standard time-slicing procedure to construct 
an object that satisfies this requirement. To this end, we divide the interval [t',t\ 
into M = {t — t')/e time intervals, with tn = t' + ne for n = 0, . . .M, and write 
j,a _ ^a(i^) [^.g = rj, rlj — ri] and p'^ — p°'{tn)- Then the following construction, 
illustrated in the first row of Fig. IDll has the desired properties: 



I?i?/l?Pe(*''"/'')^°[-^°-P"l . (D.lb) 



The second line, with action S*" = eX^n -^n; is a formal shorthand for the detailed 
time-slicing construction of the first line. Here and below, t > t' , the index value 
a = F or B should be used for the upper or lower term in the curly bracket, 
and Sa stands for sp/B = ±. The multiple products in Eq. (ID. lap contain one 
momentum integral (M in total) for each interval, and one position integral (M — 1 
in total) for each boundary between intervals (see Fig. lDip . The Lagrangian Z° and 
Hamiltonian = H"^ (tn, Rn, Pn) associated with the n-th interval are given by 
(here = hp-): 

H^P^-^-^- , (D.2a) 
K = J d(50e-*^"P"^<iJ'^(i„, RI + Sail - ba)6r'^, R^ - Sahjr^^) . (D.2b) 
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(a) U^(t.f) (b) U^(t',t) 

Fig. Dl. Three representations of the propagators (a) U^{t,t') and (b) U^{t',t), with t > t' . 
Arrows point from the second to the first index of propagators. The first row illustrates the position- 
momentum (dots-squares) time-sliced path integral representation of Eq. I ID.lal l (with the choice 
i>a = &aB in cf. Eq. ||D.3| |. so that = r^_j); the wavy line indicates which end of the n-th 
time slice the interaction field V_(4„,rJJ_j) is attached to. The second row depicts the first order 
perturbation expansion of Eq. IID.Sbl l. obtained after performing the momentum path integral, 
using Eq. (ID.9I I to convert hy to hy. The third row shows the A^-th order perturbation term of 
Eq. l ID.llbl l. The double dots remind us that the vertices hy^- and hy-^ are nonlocal (since they 
contain factors of pnfi or pan)'- they arise from "pulling together" the two local vertices at times 
tn and t„—i of the second row of this figure into a single nonlocal vertex at time tn, with which 
we hence associate a double integration J dx^f^ or J dx^„ . The dot carrying a bar indicates which 
of these two integration variables occurs in the argument of V— (Xn), namely the one drawn on the 
side of earlier times. 



Here we introduced relative and "asymmetric center-of-mass" coordinates for the 
n-th interval. 



5< = <- <_i , K = { 5^1+ sabjrl = ( 5 - 5,(1 - ba)5rl , (D.3) 



where the "asymmetry parameter" ba is a real number with < ba < 1, which 
in general can be different for a = F or _B. The actual values chosen for ba do 
not affect any of the final results, hence they can be chosen according to taste or 
convenience, or left unspecified, as we shall do for now. It is to be understood that 
under the path integral, the notation R°-{tn) and P°'{tn) [e.g. as arguments of the 
fields V±(i„, -R°(t„)], should be interpreted as i?" and P°, respectively. 

The arguments of H"^ in Eq. (jD.2bp were purposefully constructed such that the 
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inverse Fourier transform of Eq. (jD.2bp yields 

This equation can be regarded as the defining relation for H'^ (and Eq. (jD.2bp as its 
consequence): -ff^ is the (generally asymmetric) Fourier transform, with respect to 
the relative coordinate 8r'^, oiH°'{tn), in which the position arguments and r'^_i 
occur in a time-ordered or anti-time-ordered fashion for a = F or B, respectively 
(i.e. the coordinate associated with the later time, t„, appears to the left or right of 
the earlier time, in-i, respectively). This, of course, is required to ensure that the 
path integral representation for U^{t,t') and U^^{t',t) produces time-ordered and 
anti-time-ordered expressions, respectively, as illustrated in Fig. IDll The reason for 
using a factor Sa in the Fourier transform exponentials e"*'*"'' ^ in the definition 
(|D.2bp of and its inverse, Eq. l|D.4p . is simply that the factor e*'*"'" ^ occuring 
in the latter is generated by the combination iSaL'^ in the action of Eq. (jP.lap . 
Finally, note also that H°-{R°-,P°') is independent of P" if and only if ^"(rf , r^) 
is proportional to (5(rf ~ r^)- 



D.2. Verifying the Defining Equations and Composition Rule 



It is straightfoward to verify that Eq. (ID.laP satisfies all the requirements expected 
of a propagator. We shall now first show that it fulfills the defining conditions for 
Ufj, namely Eqs. ()B.39p . and then check that it satisfies the usual composition 
rule. Since the manipulations for a ~ F and a ~ B are very similar, but differ in 
numerous minor details, we shall mostly consider the former case only. Hence, a 
will be understood to stand for F below, except when explicitely noted otherwise. 

Normalization: To recover the normalization condition Eq. (jB.39a|) . take the 
limit t ^ t' hy taking M — 1 and e 0. Then the entire path integral reduces 
simply to 

dpi 



limf/^(t,i') 



(D.5) 



Equation of motion: To recover the equation of motions for and U^^, namely 
Eqs. (jB.39bp and (|B.39c|) . add one time slice in Eq. (|ma|) {M ~> M + 1, so that 
now ~ '"m+i)' expand the corresponding exponential e^*'''''^/''^^*f+i to first 
order in e: 



U^(t + e,t'). 
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drli 
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(D.6a) 
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(D.6b) 



Here Eqs. (|D.4[) and r^^^^ — rf were used to obtain Eqs. (jD.6a[) and (jD.6b[) . which, 
in the limit e ^ 0, reproduce Eqs. (|B.39b|) and (|B.39cp . 

Composition rule: Next we check that Eq. (jP.laP also satisfies the usual com- 
position rules for propagators, namely 



dx, C/fi {t, h) {h,t') = C/| {t, t') , j dxi U^, {t\ h) U^, {tut) = (t', t). 

To this end, let Mi be the number of intervals between ti and i.e. write ti =t' + 
eMi and ri = tmi- Then, by concatenating two expressions of the form Eq. (jP.lap 
for and U^, , we find that the left-hand side of the above equation can be 



written, up to a factor So- 



[dr[ lim TT 

—Ml 
Ml 

n( 



Mi-l 

X n 

n=l 



n=Mi + l 
Ml 

n=l 



dri 



(27r)° 



M 

n 



n=Mi + l 



(D.7) 



This is equal to U^{t, t') as given by Eq. (jD.la|) . since J drf = J dr^j^ . The deriva- 
tion for Uf^ is entirely analogous. 



D.3. Power Series Expansion in h"" : 

The power series expansions of U^{t,t') and Ui^{t',t) in powers of hy and hy 
are given by Eq. (jB.48p . To illustrate how they come about from the time-slicing 
definition (jP.laP of the path integral, we begin by considering only the first order 
terms (the higher order terms will be discussed subsequently). To this end, we 
expand each factor e^'*"'^^''^^"! in Eq. (jP.lap to linear order in hy^^, to obtain 
^{^s^e/h)Ll'^^ _^g»'*aPS, <5<, (-s<,ie/?i)/i^„^ . Here L^^^l is the ^-independent part of Z^^, 
and for the second term, all contributions of order or higher were dropped (in 
particular, we replaced e~^**"'^^'''''*°"i by 1). Then, to leading order in e, Eq. (jP.laP 
readily yields the following expression: 



m{t,t')~mt,t') 



umt',t)~mt',t) 



{~Saie) 




(P.Sa) 
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(D.8b) 



7S /■* f f^^l^f/llf^T 



dh dx,A , (D.Sc) 



in agreement with the iV = 1 terms of Eq. (|B.48|) . For Eq. (jP.Sbp , which is iUustrated 
in the second row of Fig. IDll we have evoked Eq. (|D.4p to make the identification 

From the above excercise, we extract the following rule of thumb: when a function 
f°-(ti) = f°-(ti,R°-{ti),P°-{ti)) [e.g. hy above] occurs at time ti along the forward 
or backward parts of the Keldysh path integral / VR^VR^ , the / dp°^e*^''^"i ''^''"i 
momentum integral at that the corresponding time slice tm — ti converts it into 
/?fini-i(^ni) /ifi-irii (^"i)' Combining this with the propagators imphcit in e*'*°'^o , 
generates terms of the form U^.^ f^^U-^ or U^^ f^^U^f , respectively [where /° and 

/° are Fourier transform pairs, in analogy to H"" and H'^ of Eqs. (jD.2bp and (|D.4p ]. 
To be explicit, we have 

5„„. Pil"^ /pP'^e^""^" =/ dxii <^ . (D.IO) 

JB,^it,)=r, J J ' [UlJfJj^, 

Having found the rule (jP.lOp , it is straightforward to go beyond the first order and 
to recover the full perturbation expansion from the path integral: 



= S^^^^ J VR^ih) y"pp"(ti)e±[*'5o (*.*')-//' (D.Ua) 

X^^^ / dh dh... dtNhl{ti)h''y{t2) . . .h'^vitN) 

pii nil piN-i r 

/ dti / dt2 . . . / dtjsi I dxi idx2^2 ■ ■ • 



(D.llb) 
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oo 

E 

N=0 



dti . . .dt 



N 



dxi I 



. dx 



il Vll 12 



N.N 



^jN VNN ■ 



If f^R 

■ ■ WNN^Nj 



(D.llc) 



''21'Vll li 

Eq. (jD.llbp . which is iUustrated in the third row of Fig lDll was obtained from 
the Une preceding it by multiple applications of the rule of thumb (jD.10|) . and 
reproduces the expansions of Eqs. (jB.48p . For Eq. (|D.llc|l we recalled Eq. (|B.45P 
to set JJ^^/^ = ±i hG^/^ along the forward or backward contours, respectively. 



D.4. Coordinate- Space- Only Path Integral 



Since the power series expansions (jP.llbp for Uf^ do not contain any explicit mo- 
mentum integrals, they may be used as starting points for deriving coordinates-only 
path integral expressions containing no / VP"- integrations at all, so that only the 
coordinate integrations J Pi?" remain. To this end, we simply perform the J DP"" 
integrals in the definition of the free propagators Uf"" explicitly, with the well-known 
result: 

M 



lim 

Af^oo 



M-l 

n 



dpi 



X exp 



M 

h ^ 

n=l 



f^Pn 



5rl h^p"^ 



2m 



(2n) 



Vhnp (Rn) 



(D.12a) 



Md/2 M-l 

n 

n=l 



dK exp 




Sr"- 



25_Rag(is„/ft)So (t,t') 



SS{t,t')[R^t,)] 



dt. 



iTOi?'^2(t3)-14mp(i2"(i3)) 



(D.12b) 



(D.12c) 



Here Sq is the standard action for a noninteracting electron in a disorder potential, 
and the tilde indicates that [in contrast to iJ" of Eq. (jD.2bp ] it is a functional of 
R°-{t3) only, not of P^-its) too. The tilde on JVR in Eq. (ID.12bj) indicates that 
the measure includes the prefactor in the line above it. Now, if we take the power 
series expansion (jP.llbp for U^j and insert Eq. (jD.12bp for each occurrence of J7|°, 
we obtain for a coordinate-only path integral expression with a precise (though 
cumbersome) time-slicing definition. In the main text, we have used for the path 
integral so obtained the formal path integral notation (|B.56p . with actions defined 
by Eqs. (|B.57p and (|B.58[) . and measure / VR, where the prime reminds us of the 
double position integrals J dxi^i occuring in Eq. (jP.llbp . The points discussed after 
Eq. (jB.58|) in the main text all follow directly from the explicit construction given 
above. 



D.5. Explicit Expressions for H° 
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The material presented up to now in this appendix was general, applicable to any 
nonlocal Hamiltonian of the form H°j = Sy hoj + hy^^ . Let us now be more concrete 
and specialize to the Hamiltonian defined by Eqs. (|B.36p . in order to verify GZ's 
expression for the effective action derived for their / VR J VP path integral. 
Inserting Eqs. (jB.36P into Eq. (jD.2|l . we readily find that 

= ho{Rn,P:) + Kitn, to; K,Pn) = hon + , (D.lSa) 

Kn = E ^""(*«' ^n' ^n) K) , (D-13b) 

a=± 

ho{R\P'') = ^ + V,^p{R^)-fM, (D.lSc) 

!f{t,t„;W,P') '■p<'">(R° + ,,„(l-t„)r»,fi°-»„6.r°). (D.lSe) 

Here ho{R°',P'^) and p'^{t,to; R°- , P"-) are, respectively, the free Hamiltonian and 
the single-particle density matrix (in the presence of interactions but without source 
terms) in the mixed representation. In the definition (jP.lSdl) of w°'~ , it is to 
be understood that i?'" should be equated to il" after evaluating the action of 
the exponential differential operator on the function V-^{tn, R^) to the right of 

{t,to; R"" , P"") in Eq. (jD.13bp . and all equations derived therefrom. 
For general choices of ba, the shift operator e~*^^''~''"^^^pS'^«° in Eq. (jP.lSdp 
is needed for the following reason: In the defining Eqs. (jB.36a[) for and i/^, 
the arguments of the field V_j are evaluated at and rf , respectively. When 
considering the n-th interval (for which = r!^, r?" = these arguments of 

become rf = r^_i = Rn - bpSr^ and rf = rf_i = R^ - (1 - bB)Sr^ [cf. 
Eq. (|D.3p ]. which are evidently shifted relative to the argument at which the field 
V- {tn, i2°) is evaluated in Eq. (|D.13bp . namely by an amount —Sa{ba — SaB)Sr'^- 
The exponential shift operator implements this shift [as can be verified by inserting 
Eqs. ((DT3|l into Eq. (1^4]) to recover and iJ^f]. Evidently, though, one can 
achieve R^ ~ fn-i{~ '"'j) R^ = r^{= rf) and hence avoid the need for 
shifts, by making the special, "maximally asymmetric" choice ba — 6aB- Indeed, 
for this choice, which we shall adopt henceforth, the exponential shift operators 
g-4(6a-'5oB)Vp- Vr reduce to unity. Moreover, since then depends on only one of 
the position coordinates r'^ and r°_]^ associated with the n-th time interval, namely 
the second, so does = H°-{tn,r'^_i, Pn ), which greatly simplifies subsequent 
manipulations. The "price" to be paid for this simplification is not high - one merely 
has to remember that the definitions (jD.2bp of in terms of the Fourier transforms 
of i?" and with respect to the relative coordinate become fully asymmetric: 

H^^jd{5r^)e-^p'- '<H^{t^,r^_,+5r^,r^_,) ^ho{rf,_„p^) + h^^ , (D.14a) 
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ijf^y'rf(<5rf)e'P"-^'-"i7^(i„,r,li,rf_i-<5r,f) ^ hoirf^_„p^) + h^^, (D.14b) 

where ft^„ = Eq=± ^n" ^a(^nj<-i)> with = e and = e^Sa[l~2p^], and 
p° is defined in terms of p^^^Htmta) by Fourier transform relations [Eq. (|D.13eP ] 
that are analogous to those [Eqs. (|D.14p ] for H'^ in terms of H°j{tn)- 

D.6. GZ's Effective Action in Position- Momentum Representation 

Having specified the time-sliced versions of -ff" in the position-momentum repre- 
sentation, it is straightforward to also derive the effective action (15*7^ + 5/)[-R°, P"] 
for this representation: simply repeat the strategy followed in App. B.5.5 to 5.8, 
but use the position-momentum representation (denoted by bars instead of tildes) 
throughout. Since the details are very analogous, we shall be very brief, and indicate 
only the main differences. 

The starting point is again Eq. (jB.59p for (^^/ 22' 21' (^i' '•2' ^o))y ns' with 

the coordinates-only path integral measure ^ 1^ 'D'{R) replaced by a position- 
momentum path integral measure, ^ ^ ^{RP), which is a shorthand for 

m V{RP)...= / VR^'it^) VP'"{t3)e''^oit^.W/r^ (D.15) 

/ I?H^(t3) / PP^(t3)e-<(*"*^^/''... . (D.16) 

SQ[R°'{t3), P°-{t3)] in the weighting factor is the action for a single, free electron, 

SSit^,tJ) = /*' dt3 [P'^it^i) ■ dt.R'^its) - /io(H'^(i3),P"(i3))] , (D.17) 

and the bar on 5*0 (and B, Sn/j below) indicates that [in contrast to Sq, B, 
Sfi/i of App. B] they are functionals of R°'{t3) and P"'{t3), not of i?°(i3) only. 
In Eqs. (|B.6ip . Bas is replaced by 

F„(i3,r-3) = ^s^T^aXra - R^its)) , W3+ = e , (D.18a) 

a 

Wl~ =e^Sa[l~{e32+lf023)p''{t3,to;R''{t3),P''{t3))] . (D.18b) 

Now use precisely the same set of approximations and arguments as in App. B.5.6 
to B.5.8 to derive the effective action iSji + Sj. One readily arrives at an equation 
just like (lR83|, but with {iL^/U) of Eqs. fRM)) replaced bji 

(zZ«/L^)3,4„, = \saSa'Wt[wi^{2Rli)3^i^, , (D.19) 



X 



JThe 934(iR/ 1)34 occuring in Eqs. IIB.841 1 was written as ^(2iR/ 1)34 here, exploiting the symmetry 
I34 = l43- 
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where the density matrix occuring in W now is the free one, p^. Muhiplying out 
the terms in Eq. (|D.19p explicitly (and setting {O4 ,2 + 024 ,) — ^ for reasons 
explained in footnote [h] on page [55)1 , we find 




SR{ti,to) = — dta dtA (D.20a) 

^ Jta Jto [ 

1 - 2pi (R^iU), P^iU))] (i?[i34, R^its) - R^iU)] - R[h4, R'^ih) - R'^iU)]] 
1 - 2p^{R^{U),P^{U))] (7?[t34,fi^(t3) - R'^iU)] - R[h4,R''{h) - R'^iU)]'^ 

Sl{tl,to) f dt3 f dtJi[t34,R^{t3) ~ R'^iU)] - /[t34,i?^(i3) - ^''(^4)] 

-i[t34,R^{t3)-R''(M)]+i[t34,R''{t3)-R^{t4)]} ■ (D.20b) 

This reproduces GZ's expressions for the effective action, since Eqs. (|D.20p are the 
analogues of (GZ-II.54) and (GZ-II.55) [our 1st, 2nd, 3rd and 4th terms, having 
aa' = FF, BF, FB, BB, correspond to GZ's 1st, 4th, 3rd and 2nd terms, respec- 
tively]. The only difference is that in their Pauli factor, GZ have evidently replaced 
our by n{R°-{t), P^it)) , which they define as the Fermi function 

n(ho), evaluated at energy /io(-R°(i), -P°(i)) ■ 

GZ offered no justification for the latter replacement in but have de- 

fended it in subsequent papers^ by arguing that it amounts to a quasiclassi- 
cal approximation that neglects terms of order h. We have argued in a previous 
publicatioiiJ^ that the "small parameter" that would protects this approximation 
is actually rdfi/T, which evidently is not small in the T ^ limit of present in- 
terest. Much more alarming, though, is that when averaging over all self-returning 
random walk paths, GZ proceed to make the assumption that "no depends only 
on the energy and not on time (our emphasis), because the energy is conserved 
along the clasical path" [see discussion after Eq. (GZ-II.68)]. As argued in Sec. 4 of 
the main text, however, this neglects recoil, and produces incorrect results. A more 
accurate way of treating the Pauli factor, that properly includes recoil, is discussed 
in Sec. 3 of the main text. 

Appendix E. Diagrammatic Keldysh Approach 

In order to facilitate comparison between GZ's notation and our's, this appendix 
collects some standard definitions (following Rammer and SmitljH) and results for 
electron and field correlators used in the Keldysh approach. [Where relevant, GZ's 
notation is given in brackets.] Below, subscripts i are abbreviations for (ti, Xi) when 
used for fermion fields or for {ti,ri) when used for interaction fields. Gij is a short- 



J. Rammer and H. Smith, Rev. Mod. Phys. 58, 323 (1986). 
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hand for Gij{tij) = G{tij;Xi,Xj) [and similarly for Cy], i.e. the time-argument, when 
not displayed explicitly, will be understood to be tij = U — tj. As elsewhere, the 
tilde signifies the matrix structure in coordinate space, while bold symbols are used 
for matrices in Keldysh space, e.g. Gy. 

E.l. Electron Correlators 

We begin with the electronic Green's functions Gy, and consider for the moment 
only those for free, noninterating electrons (i.e. evaluated for = 0): the basic 
correlators 

Gt= {{'^\{tj,xj)i'i{tuxi))a [=G^}(ij% (E.la) 



G>j^-'-{i,,{ti,xS\{tj,Xj))o [=G'i^{jo% (E.lb) 



are used as follows to construct the time-ordered, anti-time-ordered, retarded, ad- 
vanced, Keldysh and contour-ordered Green's functions, respectively: 

Gl^ e{ty)G>^+e{tji)G^ [=G'{^m, (E.2a) 

Gl^ e{tji)G^+6{tij)G^ [=G'i}{ij)l (E.2b) 

e{ty){G^-G^) NG^'^'Czj-)], (E.2c) 

G^ ^ -e{tji){G^ - G|) [= (E.2d) 

Gij = G"^ + G^ = G^ + G'^ 5 (E.2e) 

G^ = G^ — G^ = Gfj — G^ , (E.2f) 

G^ = G^ + G^ = G^ ~ G^ + '^G^ ) (E.2g) 
£fSsi(G«-G^)=i(G?-G5), (E.2h) 

«^-{| l:: tzt: '^-^^ 

where ti >c tj means that t,; is further along the Keldysh contour than tj, and %■ 
denotes contour-ordering along this contour. (The Keldysh contour runs from the 
initial time to to -|-oo and back.) Under complex conjugation, the following relations 
hold: 

{Gf^r = Gf/'', (Gfr = -Gj|, (G</>r = -G'</>. (E.3) 

It is customary to represent the contour-ordered Green's function G| by a 2 x 2 ma- 
trix Gij in Keldysh space, whose components are the quantum-statistical averages 

of contour-ordered operator products, 

{Gir' ^ {d;f)o , C!f ^ [ri{ti,xi)^i\tj,xj)\ , (E.4) 

and arc labeled by indices a, a' that take the values F and B, with the convention 
that if a = F (or B), then ti resides on the forward (or backward) part of the 
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Keldysh contour, and similarly for a' and tj . In matrix notation, we have 

where we used a boldface notation for the fermion fields to indicate that it has two 
components in Keldysh space, ipi = (^b|("!^'|)- [Note that (G.y)"'"' corresponds to 
GZ's G'^^,{ij), with F ^ 1 and B ^ 2]. A more convenient, since tridiagonal, form 
is obtained using the representatiorQ 

EE Lt^, , ;0j = , 4- = ' (E-6a) 



G| ^ (G,)o = Lr^GlL^ = (^^ ^ j , (E.6b) 

where t^''^^^ denote the Pauli matrices acting in Keldysh space, L = and 
Eq. (|E.6bj) follows from the definitions (|E.2p . 

For future reference, note also that density operators nf-j{ti) located on the 
forward or backward branches of the Keldysh contour have the following represen- 
tations (suppressing the time argument), for a = F, B: 



%i = V'f C = ^jJ'aV', = ±^jPa±^j , (E.7a) 



Pp^B - -(1 ± t3) , = Lt^Pf/bL^ = -(ri ± 1) . (E.7b) 



E.2. Field Correlators 

Next we consider the "interaction propagators" Cy, i.e. correlators involving the 
real, bosonic fields Vi that were introduced via the Hubbard-Stratonovich transfor- 
mation (|B.28ap . Below we shall use as a shorthand for Va{ti, ri), taking it to be 
understood that if a = (or B), then ti resides on the forward (or backward) parts 
of the Keldysh contour. The basic correlators 

C< ^ '-^(V,V,)v ^ C> , (E.8) 

are averaged over all field configurations according to Eq. (|B.29dp . The definitions 
of the correlators and £|- in terms of and £^ are identical 

to those of the corresponding electronic Gij's in terms of G^ and G^ in Eqs. (|E.2[) . 
The matrix representation Cy of the contour-ordered interaction propagator 
with matrix elements 

= —{%Va^Va'j)v = , (E.9a) 



' A. I. Larkin and Yu. N. Ovchinnikov, Zh. Eksp. Theor. Fiz. 68, 1915 (1975) [Sov. Phys. - JETP 
41, 960 (1975)]. This is also the form used by Rammer and Smith (see footnote Ikt . 
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takes a form analogous to Eq. (jE.Sp . namely: 

!£!f , ^ , ~. , . (E.gb) 





n \{VB^VFJ)v {VB^VBJ)v 

Following AAqI^, we shall use the transformation i = J) to obtain a tridi- 

agonal representation, reminiscent of Eq. (jE.6bp , 

with matrix elements to be denoted by , where a, a' take the values ±. The 
last equality of Eq. (|E.10a|) was obtained by using 

(E.ll) 

[cf. Eq. (jB.37[) ] to rewrite LCyL^ in terms of the correlators e'^{VaiVa'j)v- The 
relations (jE.lOap are general. The explicit expressions for these correlators given 
by Eq. (jE.lObp . which are specific for the present model, follow from Eq. (jB.TSap . 
[Incidentally, comparing Eqs. (jE.lOap and (jE.lObp proves Eq. (|B.75b| . Using the 
explicit forms for Rij and lij of Eqs. (|B.76p . it can easily be checked that 

and that their Fourier transforms w.r.t. ty satisfy the relations 

£f (u.) = Cf,{c.) = £f (-^) - >Cf (c^) , C^ic) = £j^M = -Af (c^) , (E.13) 
C^iuj) = coth(nc./2r) [£f (a;) - C^{u;)] . (E.14) 

Eq. (|E.14p [cf. Eq. (|B.72p ] has the form required by the fluctuation-dissipation 
theorem. 

Explicit expressions for the interaction propagators are most readily written 
down in the Fourier representation. For disordered metals, where small frequencies 
and wave numbers dominate, we obtain from = —e^R and (jB.77p the following 
relations (in ag reement with Eq. (5.8) of AAcfl^: 

C^iu;) , >Cf (c.) = 2zcoth(nc^/2r)In<(c.) . (E.15) 



E.3. Keldysh Perturbation Theory 

In this section, we recall how the Feynman rules for Keldysh perturbation theory are 
derived, and use them to obtain an expression for the self energy S of the Keldysh 
electron Green's function [Eq. (jE.24p ]. 
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In the Keldysh approach, expectation values of the form occuring in Eq. (jB.27bp 
are written as follows (following Rammer and Smith, see footnote |k| ■ 

l^O{t))K = ^^^B{t,to)di{t)UjF{t,to))n ^ {%S„Scv di{t))o ^g^^ 

Sci = Tc e^T^^o Hu{t^) ^ = 7-^ g-i rft3 fi/(t3) ^ (E.16b) 

where dti and Tc indicate integration and time ordering along the familiar 
Keldysh contour [H^I and v are defined in Eqs. (IB.lip and (IB.27dp ]. In Eq. (IE.16ap . 
the operator Oi[t) can be written as either Of (i) or Of{t), where the superscripts 
indicate that the operator resides on the upper or lower branch of the Keldysh 
contour, since the contribution from the portion of the Keldysh contour from t to 
00 cancels that from 00 back to t. Consequently, we can also represent Oi{t) as 
\\0^ + Of\{t), which turns out to be most convenient and will be used henceforth. 
For examples, the reduced single-particle matrix p-^j^,(t,to) of Eq. (jB.27bp can be 
written as^ 

~p,y{tiM) ^W^vi + nfvi)K = (ii, hr^±,)K - -iKlTK [\T^&^i^] , (E.17) 

Here Ttk denotes a trace over Keldysh indices, = {Gij)K (and likewise — 

{Q-iv) K,\isi which will occur below, too), G ij has thejsame matrix structure as in 
Eq. (jE.6ap . and the superscript "full" (or "ns") indicates that the average is to 
be evaluated in the presence of the full interaction and including (or excluding) 
all external perturbations, i.e. with ( )k (or ( )K.ns) instead of ( )o. As a check, 
we note that in the absence of interactions, Eq. (jE.17|) reduces to —ih^G^y = 
^(V'J/V'i ~ V'iV'i')o ^ which is equal to the desired result of (V'i''0i)o (recall that 
V'J/ and V'l anticommute, since is equated to xi only at the very end of the 
calculation) . 

By writing J^dt^viits) = dt^lv^j — v^j], and switching to the Keldysh rep- 
resentation of Eq. (jE.6ap . Scv takes the form 

*33(i3) = n3ih)[n33j - n^3j] = ^'33(^3) V^glV^3 ■ (E.18b) 

As a special case of Eq. (|E.18bp . we note that the external perturbation, iJext of 
Eq. (jB.17|l . generates vertices of the form 

{-t/n)fQ' = {-z/h)h''^.il,l\xl,^ . (E.19) 



"^An alternative but equivalent form to Eq. IIE.lTI l is often used (e.g. by AAG-'^^, Eq. (5.1), where 
the factor 2 in front of fi is a typo), namely p^^, = ^(x-"- — 1)^^)^ = —ih\\{T^ — 1)G^"*,*] , 
where it is to be understood that t-^i = ti -|- 0+. 
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To linear order in h'^^^, where each fermion line is simply decorated by the insertion 
of a single external vertex, we thus have the Feynman rule that each full G^j"'' is to 
be replaced by 

iE{ujo) ■ 322' 



^22' '^G2j)K,i\s 



{G_i2' 'i-G2j) K.ns 



(E.20) 



U!0 

where the subscript "ns" denotes "no (external) sources" , and the term in brackets 
indicates the form which /i^f^ assumes under Fourier transformation, if we use the 
gauge of Eq. (|B.21bp . 

For any expectation value of the form {0{t))K, the interaction term Hn in Sd 
can be decoupled using the Hubbard-Stratonovitch transformation of Eqs. (|B.28p . 
just as in Sec. lB.3[ using the fields Vp and Vb for the forward and backward branches 
of the Keldysh contour, respectively. One then readily finds that {0{t)) k can be 
expressed as follows as a functional average over all fields Vp/s'- 

{di{t))K ^ {0^{tM))v (E.21a) 



iV/. . \ _ {%ScV Scv Ol{t))Q 



{t, to) ^ ' Y{t to) ' ^^-^^^^ 

Z(t,io) = (r,4v'5,D)o , (E.21c) 
S,v = rce"^^*r''*-'^'^"^^3 ^ (E.21d) 

= e[n^^,VF{r^) - ng^VB{r^)\ - e[V^3,(l^+3 + \T^V^:i)±,] ■ (E.21e) 

Here the functional average ( )y over all field configurations is defined, as before, 
by Eqs. (|B.29dp . where the functional Z occuring in Eq. (jB.29e|) is now given by 
Eq. (|E.21c|) . 

To obtain an perturbation expansion within the Keldysh approach, one expands 
Scv in powers of (— i/^)il3i which thus serves as basic interaction vertex, and 
then applies Wick's theorem to the fermion fields. In the n-th order term, there 
are n! equivalent ways to connect the n vertices with n fermion lines of the type 
{'ip.ip\)Q = ?fi.G°, yielding a combinatorical factor of (ifi,)"n! which cancels the 

{—i/h)"/{n\) from the expansion of the exponent of Scv- Next, the average ( )v 
over all field configurations is to be performed, which yields contractions between 
the interaction fields pairs of vertices. These contractions have the form 

{V,V,)v = -5*^EiI'^"i.^r'i-'7"'i,- , (E.22) 
where we introduced the "vertex matrices" 7"'' = 1 and 7" — t^, the field propa- 

— act' 

gator in the Keldysh representation is given by Eqs. (jE.10[) . and the Feynman 
diagram corresponding to Eq. (|E.22[) is the leftmost graph in Fig. lEll Eq. (jE.22p 
implies the following Dyson equation (cf. Eq. (5.6) of AAcfl^, 

G^f = GO + / dU fdx3 dxi G% S34 G^f , (E.23) 
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Fig. El. Feynman diagrams for the interaction propagator [Eq. | |E.22| |] and the correlators 
Jff, Jbb and Jfb+bf [Eqs. l)E.28| l] that give the leading correction to the conductivity due to 
electron-electron interaction. Solid lines denote matrix Green's functions , wavy lines interaction 

propagators , and the symbols a and a' the vertices •7°' and 'y°' . Arrows point from the second 

to the first index of propagators. 

where, to lowest order in the interaction, the self-energy is given by 

S34 = -^^hJ'^G!i^'r"'l34 ■ (E.24) 



E.4. Conductivity 



In this section we derive a general expression for conductivity ctdc hi the Keldysli 
approach and expand it to leading order in the interaction propagator. This will 
allow us to check the perturbative expansion (jC.lip of our influence functional 

^12', 21' of Sec. EH 

We start by using Eq. (|E.17p to express the quantum-statistical average of the 
current density operator J7ff(ii,ri) of Eq. (IB.16P as follows, 



jiv A{ti,ri) 

m 



(E.25) 



Next we expand Eq. (|E.25p to first order in /loxt [using Eq. (jE.20|) ]. and then 
use Eq. (|B.20p to calculate ctdc; the result has the form of Eq. (jB.22bp . where 
Ji2',2i'(^o) therein is given by the Fourier transform w.r.t. ti2 of the following 
expression: 



■^12', 21' — "-LrA- 



I'T {G_12'G21')ka 

In the absence of electron-electron interactions, this readily reduces to 



(E.26) 



x(0),Kcldysh 
•^12', 21' 



\h{Gi2'G2l' + G'^2'G21'^ — ^^{Gi2'G21' + Gf2'G2\'^ ■ (E.27) 



The second equality follows from Eq. ( |E.2gP (with G^^'^G^/^ — 0) and confirms 

Eq. (icia . 

Let us now obtain the leading correction to ctdc due to the electron-electron 
interaction. To this end, we have to expand Eq. (jE.26p for J^2'^'^f^ to second order 
in V^g. One readily arrives at the following result [which can also be obtained by 
starting directly from Eq. (|E.25P , expanding G^"'/' therein to first order m using 
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Eq. (jE.23|) , and then expanding each G| in the latter equation to first order in /lext 
using Eq. jESOl]: 



T(2),Koldysh 
•^12', 21' 




dt^ dt4 J dxz dx4 (^Jff + Jbb + 


jBf\ 
■^FB J J 




Jff 


aa' 


1^1 r'O .,Qr^O „,a'/=.0 /=.0 /^""' 
2^ LlisI ti347 ^£42,0:21' i:i34 




(E.28a) 


Jbb 


aa' 


5'''^^12'^237"^§47" ^41' :^34 




(E.28b) 


JBF 

'Jfb 




2^ iil37 ^3132/0:247 Llu' L.3i 




(E.28c) 



aa' 



The correlators Ji? ^ , Jbb and Jj?^ are illustrated in Fig. lEll and correspond to self- 
energy insertions in the upper and lower Keldysh contours, and a vertex correction, 
respectively. Multiplying out the Keldysh matrices explicitly, taking the trace and 
omitting all terms involving the combinations G^^Lf^ or Gf^Cf^, which vanish (since 
^34^43 = 0), we obtain: 

Jff — ^Gf^ [G^^C^^ + Gf^Cf^] (G'^2'G'2i' + G^/G^^/) (E.29a) 
I 1 iriR (r<K pK , /^_R pR , p^A r-^\ ^ (n^ . ^^k pA w /^A 

"•"2 L'-'13l<^34-'-34 + '-^34-'-34 + <-^34-'-34j + ^131'-^ 34^34 + '-^34-'-34jJ W2''-^21' J 

Jbb ^ ^{Gi2'G23 + G^2'G23)[G34:^34: + G34C^4\Gfy (E.29b) 
I 1 r<R \(nR rK . r<K pR \ riK , //^k pK . r<R pR , riA pA \ riA n 

+ 2*-'12'^23 Ll'-'34''-34 + '-^34''-34 j '-^41' + l'-'34-'-34 + *-^34''-34 + ^34-'-34 j Wl' J ; 

JBF _ 1 (f^R /^K I /^K I^A \ r^A I^R pK , /^iR pR . ^A j^A ] ^t;, r,Q„N 

•JfB — 2 v'-^32'<-^24 + '-^32'<-^24j LWl'<-^13''-34 + '-^41''-'13^34 + '-^41''-^13^34j 

, 1 riR riR \nK nR pK . ^r. pR , f^K < aa \ pA ] 

+ 2'-'32'<-^24L<J'41''-^13''-'34 + "^41' "-^13 ''-34 + l^41''-'13 + Wl'<-^13j''-34j 
I 1 AA AA rAA AiiT pK , AA A_R pA , fr<K AK , A_R AA \ pR ] 
+ 2^32''^24 L'^41'*J'13-'-'34 + "-^41' "^13^34 + 1,^41'^13 + "-^41' "-^13^^34] ■ 

Now, terms that involve the combination G^/G|^ or G^G^ contribute to the so- 
called interaction corrections, and do not contribute to "decoherence". Hence, we 
retain only the first lines of Eqs. (|E.29p henceforth. For these, we use the identity 
[cf. 

\{G^G^ + G^G^) = G'^G< + G<G^ + i(G«G« - G^G^) (E.30) 

and drop the last term, for the same reason. The remaining terms then take the 
following form: 

J_FF=Gf3 [Gf4£^ -I- G^£f4] (Gf2'G^i, + Gf2'G2v) (E.31a) 

JbB = {G^2'G23 + ^12' ^23) [^34-^34 + ^34-^^4] G^y (E.31b) 

JFB={G32'G2i + G32/G24) [Gfi,Gf3£|^ -I- G|i/Gf3£f4 -I- Gfi,G^£34] . (E.31c) 

These expressions agree with the expansion (IC.11|) we obtained from the influence 
functional approach, as can be seen by relabelling 3 <-> 4 in some terms. [Jp^ here 
accounts for both Jbf o-nd Jfb there.] 
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Appendix F. Diagrammatic Disorder Averaging 

In this appendix we summarize, for reference purposes, some standard and well- 
known conventions and results used for diagrammatically performing disorder av- 
erages, using notations summarized at the beginning of App. [B) 

F.l. Definitions, Standard Results and Useful Tricks 

To perform the disorder averages, we take the impurity potential to be short-ranged, 
Ri)^ (^imp has units of energy times volume), represent the 
fermion fields as ipa-{t, r) = VoP^^^ e'P ''Cp(t), and Fourier transform as follows: 

Gf ^ = G^/^{t„x.,x,) ^ 5..., I ge-*-^^ 5: ^iP.r.-Prr.)G^/A^^^ . (F.l) 



PiPj 



Using standard diagrammatic techniques, the disorder-averaged single-particle 
propagator is found to have the form [Fig. IFlf a)]: 

oris) = /dr,e--^e--(Gr>.. = ns-^,Ln/2r.^ ' ^^'''^ 

Here = h / {2Tn' Ci„-^p vf^^) is the elastic scattering time, Cjmp the impurity con- 
centration, = p^It? /2m — ep, and calligraphic symbols will be used throughout 
for disorder-averaged quantities. The corresponding position-time expression, found 
by inverse Fourier transforming, is: 

Q^'^it) = J (dp)e^P--^ J {de)e-^^%'^{e) , (F.2c) 

^ie-^t/h ^-\t\/2T^, . (F.2d) 



d/2 



2nt 

The disorder-averaged products (G^G'^)dis have the form [Fig. IFlT b)]. 
(G^,p(e)G^,p(e-))^.^ = V,P V,P^p(£)^i(^") 

+ V+P'.P+P Qp' (e) (e) Qp' (e) Qp (e) Vo\2iTVT,^/h ' ^ ^ ^ 

Cq(a;) and ■P^(w) being the bare (i.e. without interactions) Cooperon and diffuson, 
respectively. Fig. IF II summarizes the standard calculations of Cq(w) and I?q(tj), 
and of the diffuson-dressed interaction vertex Tqiui) and polarization bubble Xq{'^)^ 
which is defined as the Fourier transform of Eq. (|B.66a[) : 

{xqH)^^ = -i2e^h I (de)(dp)(^p+q(£ + c.)g<(e) + g<+g(e + c^)^p(£)) jF.4a) 
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(a) 



(b) 



R/A, 



XRJA 



R/A,z,p 



P 

P' 



kk = 



1 -i^{-Rjil2iIji) 



1 



dis 



R,£.,p 
A,e,p 



7 AIK tt>,? , A/K 

R,£+(i),k 



(dky[ A,t,q-k 



r^(co) 



"p +p,p+p 

Vol 



R,£,p 



A,£,p' 



e-E,P +P 
n 



R,£,p' 



A,^,p 



R,^.p 
A,e,,p 



e-£,P +P 
ij i 



R.E,P' 



A,e,p' 



(C) 



c;(co) 

2jivi^i//i 



£+a,p e+cop e+co,/7 



co,9 III 

iiiUliiiiilll 



2jivi^i/A 



e+cop e+cop' 



/;,e+co,A; 



e+co,p' 



^£.q-p' 



E+O.p' 



\P'-P + ik-p A,z,k-, 



iiiiniiiiiiiii 
itiiilli 



'^Z,p'-q 



(d) 



r,(co) 



CO, ^ 




e+coj' 



(e) X,(co) = 



e+co, 



£+C0, 



e+co, p+^r 
A 



Fig. Fl. The building blocks of diagrammatic perturbation theory: (a) Basic definitions for the 
electron lines G^^^ and Sp^^is), impurity lines, the function Ilq{iL>) of Eq. l)F.5a| l. and the in- 
teraction lines Cij or Cq{u>) of Eq. l)B.75b| l. For all correlators, arrows point from the second to 
the first indix. creation to annihilation operators] Internal impurity momenta are to be integrated 
over with J (dk), as in nq(aj). (b) Eq. IIF.3bl l. (c) The bare Cooperon Cq{uj) [Eq. l IF.Sbl l] and bare 
diffuson (lu) [Eq. IIF.ScI l] ; (d) the diffuson-dressed vertex Fq (uj) [Eq. l IF.SdI l] and (e) polarization 

bubble Xq{^) [Eq. l IF.Sel l] . For each of Ilq(uj), C^{u)) and rq(cj), the frequency argument u) is de- 
fined as the frequency of the corresponding retarded Green's function minus that of the advanced 
one. 
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-i2e'h / {de){dp){[-u:n'^{e)]g^+^{e + oj)g^{e) 



-no{e) 



yp+q 



p+q{e + uj)gp (e) 



dis 



The results are: 



-J {dk)gl{s + u)g%{s) = ^ [1 - re,(i?g^ + 



Tcl 



1 



l-nq(a;)/(27rt/Toi/n) 

Tol 

l-liq{uj)/{2TIVT,y/h) 

^ ^ lLq{u;f^Ju) 



1 



i-KVTcY^ /h Tc\{Dq'^ — iuj) 



Dq^ 



^2 _.Drudc 
'f "DC 

Dq^ — iu) 



(F.4b) 

(F.5a) 
(F.5b) 
(F.5c) 

(F.5d) 
(F.5e) 



Here D = VpT^i/d is the diffusion constant in d = 3 or 2 dimensions, is a 
magnetic-field cutoff and the dots indicate subleading terms that are small in lutc\ <^ 
1 and qlc\ <C 1. 

For convenience, we also summarize here some results that are useful for evalu- 
ating momentum integrals that arise in diagrammatic perturbation theory. Usually, 
the energy parameter he of the disorder-averaged Greens' functions Cy^ (e) is con- 
fined to the vicinity of e-p, typically by the presence of a factor —di^riQ^he) in an 
J de integration, so that terms of order he /e-p can be neglected. [The second term 
of Eq. (|F.4bp does not contain a factor —deUo, but one can be generated by inte- 
grating by parts.] The explicit form (jF.2p for tJp^^(e) then implies the following 
"identities" : 



Gp (e) Gp (e) = 1 



2TrUTel/h P 

(dp) 



{dp) R/A -R./A 

Gr, (e) G„ (e) = 



2TTVTe\/h 



Gp (e) 



2'KVTc\/h ^ 



P 

ri-1 



171 + 11 — 2 

m — 1 



(F.6a) 
(F.6b) 



Furthermore, in the limit of small frequencies {i-u,lu <C 1/tci) and wavenumbers 
(q^, q^ ^ 1/Drei), integrals of the following kind can be evaluated by a systematic 
expansion in the small paramters, combined with repeated use of Eqs. (jF.6[) : 



(dp) 



Gp^^is)G;i;is + Lo) = ^-r4Dq^±^u;] 



-AjR 



2'KVTe\lh 



^ f (dp) -^R/A s-^AjR -.rA/R _ 

^ J 2^^i ^^^'^ + + 

^±i{l-T,i[D{q + q)^±i{u; + Lj)]} + ... , 



(F.7a) 
(F.7b) 
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(dp) 



-A/R 



A/R 



Qp (e) Qp+q{£+^)Qp+q{£ + ^) 



Tol^ J 2TlVTc\/h 

2 - T^i[AD{q + qf ±ii{uj + u))] + ... , 

2 - Tel [D[A{q + qf + A{q'f - 6q' ■ {q + q)] ± 3i{uj + oj - w')] + ■ ■ 



(F.7c) 



(F.7d) 



F.2. Cooperon Self Energy 

In this section, we provide some details for how the Cooperon self energy can be 
calculated by performing the disorder average diagrammatically. As starting point 
we use Eqs. jEMI, which we derived in Sec. IB. 6. II from the influence functional 
approach, but which are equvalent to the standard Keldysh expressions following 
from Eq. (jE.24p . According to Eqs. (|B.89|1 . there are four self-energy contributions 
to the Cooperon self energy, which we write as: 



—self _ 1 



dcu 
2^ 



(dq) 



' 

^FF + ^BB 1 ^FF ' ^BB 



The diagrams for are depicted in Fig. EDjb), those for in Fig. |F2tc) to 
IF2r f) (which correspond one-to-one to Fig. 2(b) to 2(f) of AAV"*^^). Starting from 



Eq. (|B.89|) . the corresponding algebraic expressions can be written as 

-iiZf(cs)cV(c.-<i)ri^) 



^FF- 



(F.9) 



^BB- 



-•FF- 



^BB- 



4 

iz£j((i) tanh[;.(£ + c - i^)/2T]{cl_,{cj - Co) Y^'^ - r,,Tl{u,) ^ Y^"^}, 

4 

^{C^ioj) Unh[h{s - cS)/2T]{-C;_,(a; + Co) Y^^ + T,{r%{~uj) ^ Y^^^]. 



n=2 



In the expressions for E^^^, the minus signs before Yp^\ Yp'^'' and Yj^'' arise 

from the minus sign in — tanh( ) [Q^ — and the F's represent integrals 



over internal momenta, that can be performed using variations of Eqs. (jF.7P : 



Y, 



(1) 



(dp) j^A -^R -^R _ 

Qq-p{£) Gp + ^) Qp-qi.£ + - Cj) 



1 J 2TlVTc\/'h 



(F.IO) 



[dp') -,A -^R -\-r^ I 

6q-p'{£)Gp'-q{£ + cj - tj) gp,{e 



1-T, 



D{2q - qf - i{2uj - U)) 
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p-q I p'-q 

P'^P 
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R e+co, p e+co, p' /? 

- k - 
E e-co e-co 



q-p 



q-p-q q-p'-q 



e 



-^'-5^-^Jl^V^r_,(-co) 



e, q-p A 
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e 
9-P 



e+co, p R 

e-co 
q-p-q 



e 



r_5(-co)t^ „co^ Ipjv' r_5(-co) 



r^(co)£' R CO, 9 

R R m. K 



e+co \ e+co 
P V P' 



e+co-co 
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e+co' e+co 
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P-p ^ 


R 


e e e-M 


£'\ e 
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Fig. F2. (a) Dyson equation [Eq. | |12| |] for the Cooperon, for the case that the Cooperon self 
energy contains only self-energy contributions. The latter arc shown in diagrams (b) to (f), which 
depict how to perform the disorder average of the various contributions S^^f' to the Cooperon self 
energy, leading to Eqs. l|F.9|l and ||F.10| |. Diagrams (b) depict Tjpp/gg; the diagrams (c) + (d) + 



(e) + (f) depict Y^pp^gg, the four contributions corresponding to the terms in Eqs. (IF.IOI I that 

contain Ya^^ , Ya"^^ , vi^' and vi*' , respectively. [To avoid cluttering the figure with factors of i , 
the energy and momentum labels e, u! and q used here were assigned in a less symmetrical way 
between upper and lower lines than in Fig. lC2f c): to transcribe the expressions used in this section 



into the notation used there, make the replacements (e + 



Phc 



(P2 + |q), and identify {£ + |fii)t 



e 



2^J^ Pherc 



(Pl + \<l). 



(2) f (dp) r-R _ -^A 
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2-rn 



4:D{q + qf - 3i{LU + Q) 



{dp') 



UJ) Gp-q{e + (JJ-U}) Gq-p{e) 



1 J 2TlVTc\/h 



(4) 



1 - Tel 



D{q + q)' - iito + Lo) 



{dp){dp') -A , .rp^R, , _ R ,n2 



l-T, 



i:)(2q)2 - i2uj 



}{ 



1-r, 



D{2q)'-i20j +...y 



Performing a similar set of integrals for the Y^"'''s, we readily find that Yg"''{uj) = 
Y^\—uj). Note that the sums J2n=2 ^i-zs' which are associated with the so-called 
"Hikami-box" diagrams of Fig. |F2jd) to|F2jf), add up to zero in leading order, 

/ self 

which is why the next order had to be included. Finally, the results for E ' and 
E , given by Eqs. (jl3bp in the main text, are obtained by inserting Eqs. (jF.lOp 
into Eqs. (jF.9[) and (jF.8[) . and making the replacement ehore ^ ^ ^ 5^ (cf- caption 
of FiglF2|). 
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